INVESTIGATION OF CERTAIN DIFFERENTIA] 
AND INTEGRO-DIFFERENTIAL EQUATIONS 
IN ABSTRACT SPACES 


A Thesis submitted 
in Partial Fulfilment of the Requirements 
for 'the Degree of 

DOCTOR'br PHILOSOPHY 

By 

REETA SHUKLA 
(9920861) 



to the 

DEPARTMENT OF MATHEMATICS 
INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

September, 2002 




1 6 OCT ?004 


aY^t «f MV •*'<1^ 

miftT*oA 1498^4- 


Jt-i 

M " .13 J*. 

fi , ,, -v. j 




CERTIFICATE 


It is certified that the work contained in 'the thesis entitled “Investigation of Certain 
Differential and Integro-differential Equations in Abstract Spaces” by Ms. 
Reeta Shukla (Roll No: 9920861) has been carried out under my supervision. In 
my opinion, the thesis has reached the standard fulfilling the requirement of regulation 
of the Ph.D. degree. The results embodied in this thesis have not been submitted to 
any other University or Institute for the award of any degree or diploma. 



Dhirendra Bahuguna 
Associate Professor 
Dept, of mathematics 
LIT. Kanpur 


September, 2002 



Synopsis 


Name of the Student : 

Reeta Shukla 

Roll Number ; 

9920861 

Degree for which submitted : 

Ph.D. 

Department : 

Mathematics 

Thesis Title 

Investigation of Certain Differential 
and Integro-differential Equations 
in Abstract Spaces 

Thesis Supervisor 

; Professor Dhirendra Bahuguna 

Month and year of submission 

: September, 2002 


Many physical problems can be modelled as evolution equations in Hilbert or more 
generally, in Banach spaces. The evolution equations may be viewed as ordinary differ- 
ential equations in infinite dimensional spaces of functions and are generally associated 
with the partial differential equations modelling many physical phenomena such as the 
heat conduction in materials, the wave propagation in materials, the population dynam- 
ics, chemical reactions, etc. The work presented in this thesis aims to demonstrate the 
applicability of rich theory of functional analysis to analyze the existence, uniqueness, 
continuous dependence on the initial data and asymptotic behavior of the solutions of 
such problems. 

The main advantage in analyzing such problems in abstract spaces is that not only 
we establish the results for the problems under consideration, but also for the whole 
class of problems to which these prototype of problems belong. We concentrate on basic 
features those remain invariant when we switch from one problem to another belonging 




to the same class. 

After analyzing the existence, uniqueness and continuous dependence of a solution 
of an evolution equation, we aim to find the solution. Often it is a difficult task to 
find the exact solution. In such a case, the approximate methods of analysis provide 
alternate ways of finding the approximate solutions and analyzing their behaviors. We 
consider certain approximate equations in finite dimensional spaces associated with 
an evolution equation under consideration and establish the existence and uniqueness 
of solutions to these approximate equations. The convergence of the solutions of the 
approximate equations to the solution of the evolution equation is then established. 
Our tools for this are the theory of analytic semigroups and the contraction mapping 
theorem. We consider analytic semigroups generated by the operators appearing in 
evolution equations under consideration. 

Although, our study concerns some special type of evolution equations, namely the 
Sobolev type equations and certain classes of second order semi-linear differential and 
semi-linear integro-differential equations but the techniques are quite general extendable 
to other type of problems. The results obtained may be applied to evolution equations 
arising in the modelling of many physical problems, for instance, the propagation of 
waves with strong damping, the vibrations in a beam, certain viscoelastic materials 
with memory, the flow of a fluid through fissured rocks and the shear in a second-order 
fluid and other related problems. 

The thesis has seven chapters. 

In Chapter 1, first we give the introduction to the problems discussed in the 
subsequent chapters providing a motivation to the study carried out in this thesis. 
Then, we give a review of related works in the literature and the significance of the 
problems considered. Next, we mention some of the results required for the analysis in 
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the ensuing chapters and finally, we give a brief description of the analysis carried out 
in the ensuing chapters. 

In Chapter 2, we consider the Sobolev type evolution equations in a separable 
Hilbert space and study the approximations of the solutions. We extend the technique 
of Milleta (1994) to these type of problems. The Sobolev type evolution equations 
are implicit evolution equations in which time derivatives of the unknown function 
appear implicitly and these type of equations arise in many diffusion and fluid flow 
models. First of all, we consider prehminaries and assumptions required later. We 
then consider an associated integral equation and a sequence of approximate integral 
equations using projection operators. We establish the existence of a unique solution 
to every approximate integral equation using the contraction mapping theorem. After 
proving some estimates for the solutions of approximate integral equations, we prove 
the convergence of the solutions of the approximate integral equations to the solution 
of the associated integral equation. Further, we show that the solution of associated 
integral equation can be extended to the maximal interval of existence and it is unique. 
Finally, we consider the Faedo-Galerkin approximations of solutions and prove some 
convergence results. 

In Chapter 3, we consider a strongly damped semi-linear wave equation and 
reformulate it as a second order semi-linear evolution equation. We study the approxi- 
mation of the solution to a second order semi-linear evolution equation in a separable 
Hilbert space by using similar techniques which are used in Chapter 2. First, with the 
help of a pair of associated integral equations and projection operators, we consider a 
pair of approximate integral equations and prove the existence and uniqueness of the 
solution to this pair. We then establish the convergence of the pair of approximate 
integral equations to the pair of associated integral equations and hmit of the solution 
of approximate integral equations being the solution of associated integral equation 
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after proving some estimates for the solutions of the approximate integral equations. 
Further, we show that the solution of the pair of associated integral equations can be 
extended to the maximal interval of existence and it is unique. Finally, we consider the 
Faedo-Galerkin approximations of the solutions and prove some convergence results. 

The behavior of many physical systems at an instant depends not only on the 
instantaneous values of the functions appearing in the governing equations but on their 
values at previous instants also. Such a system is called a system with a memory. For 
instance, in nuclear dynamics and thermo-elasticity, there is a need to reflect the effects 
of the “memory” of the system. This results in the inclusion of an integral term into 
the governing partial differential equation yielding a partial integro-differential equation 
(FIDE). This provides us a motivation to study the partial integro-differential equations. 
In the next two Chapters, we study a class of second order semi-linear integro-differential 
equations arising in the study of viscoelastic material with memory. 

In Chapter 4, we consider an abstract second order semi-linear integro-differential 
equation in a Banach space and first prove the existence and uniqueness of a local 
classical solution with the help of the semigroup theory and the contraction mapping 
theorem. Further, under some additional growth conditions on the nonlinear maps we 
establish the continuation of this solution, the maximum interval of existence and the 
global existence. 

In Chapter 5, we consider a strongly damped semi-linear integro-differential equa- 
tion and reformulate it as a second order semi-linear integro-differential equation in a 
separable Hilbert space and study the convergence of the approximation of the solution 
using similar techniques which are used in Chapters 2 and 3. 

In Chapter 6, we consider a first order quasi-hnear implicit integro-differential 
equation in a Banach space. First, we prove the existence, uniqueness and continu- 
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ous dependence on the initial data of a strcn, solution by using the ntethod of seuti- 
d.scret.zat.on m t.me in which we discretize the time sods and replace the time deriva- 
tives by the correspond difference quotients. Then, we show that these discrete points 
lies in a ball whose radius rs independent of the d.cretization parameters and prove 
some estimates. After defining the approxrmate solution in terms of these discrete 
pornts, we prove its convergence. Next, we establish the existence of a unique local 
mild solution and its regularity under some additional condition with the help of the 
semigroup theory and the contraction mapping theorem. 

In Chapter 7, we conclude the thesis with some remarks and provide some imtight 

into further research work in this and related areas. 


Relevant references are appended at the end. 
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Chapter 1 


Introduction and Preliminaries 

In so far as the theorems of the mathematics relate 
to the reality they are not certain, and in so far 
as they are certain they do not relate to reality. 

— Albert Einstein 

1.1 Introduction 

Many physical problems can be modelled as evolution equations in Hilbert spaces or 
more generally, in Banach spaces. In this dissertation we have studied some classes of 
evolution equations which arise in several physical and reahstic phenomena, for example, 
the propagation of waves with strong damping, the vibrations in a beam, the viscoelastic 
materials with memory, the flow of a fluid through Assured rocks and the shear in a 
second-order fluid. We have established some results for the existence, uniqueness and 
convergence of the approximation of solutions to certain nonlinear first and second order 
evolution equations. By the order of an evolution equation we mean the highest order 
of the time derivative appearing in the evolution equation. 
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After proving the existence of a unique solution to an evolution equation in an 
abstract space, the next aim is to find the solution. Often it is a difficult task to find 
the exact solution. In the such cases, the approximation methods of analysis pro^'ide 
alternate ways of finding the approximate solutions and analyzing their behaviors. We 
consider certain approximate equations in finite dimensional spaces associated with 
the evolution equations under consideration and establish the existence and unique- 
ness of the solutions to these approximate equations. The convergence of the solutions 
of the approximate equations to the solution of the evolution equation is then estab- 
lished. Similar procedure is followed to study the nonlinear Sobolev type evolution 
equations and the second order semi-linear evolution equations. Further, we prove the 
convergence of the Faedo-Galerkin approximations to the solutions whose existence and 
uniqueness have already been established. The Sobolev type evolution equations arise in 
the study of the partial neutral functional differential equations with unbounded delay. 
The strongly damped semi-linear wave equation, the Klein-Gordon equation, the dy- 
namical von-Karman equation and the vibrating beam equation can be modelled as the 
second order semi-linear evolution equations. Next, we have considered the existence 
and uniqueness of the solutions to second order semi-linear integro-differential equations 
which arise in the study of the viscoelastic materials with memory and have proved the 
convergence of the approximate solutions. Finally, we have considered the first order 
quasi-lineax implicit integro-diflterential equations and established some existence and 
uniqueness results. 

Certain partial diflferential equations are of widespread interest because of their 
connection with phenomena in the physical world. The heat conduction in materials, 
the vibrations in the wires and beams, the motion of the elastic bodies, the propagation 
of a small disturbance in a gas and many physical problems are studied by using their 
governing partial differential equations. These partial differential equations are classi- 
fied according to certain common features prevailing in each class. We concentrate on 
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certain prototype of problems, which represent each of the members of these classes. 
For instance, the heat equation, the wave equation and the Laplace equation are the 
prototype of the parabolic, hyperbolic and elliptic equations, respectively. Many prop- 
erties remain same for each member of the same class. These properties are called the 
invariant properties of the class. 

If we are interested in the invariant properties of certain class of problems then 
the best way to study such problems is to consider the abstract formulations of such 
problems. By an abstract formulation we mean a functional analytic representation 
of the problems. The abstract formulations of the prototype of problems are known 
as the evolution equations and their study comes under the area which has the AMS 
classification as Differential Equation in Abstract Spaces (34Gxx). These are basically 
ordinary differential equations in the Hilbert spaces or more generally in the Banach 
spaces of functions. In such formulations only invariant properties come into picture 
and many unnecessary details of an individual problem get suppressed. This allows us 
to concentrate on the invariant properties of a class. We may then use the rich theory 
of the functional analysis to establish the existence, uniqueness, continuous dependence 
on the initial data and asymptotic behavior of the solutions. 

To successfully exploit the tools of the functional analysis, one has to choose a 
suitable space and properly define the operators in that space. Same problem may be 
reformulated differently in different spaces. Therefore it is very important to make a 
proper choice of the space and the operators required for the abstract formulations. 
The main advantage in analyzing such problems in the abstract spaces is that not 
only we establish the results for the problem under consideration, but for the whole 
class of problems also to which the prototype of problem belongs. For example, the 
results established for the abstract formulation of the strongly damped semi-linear wave 
equation can be applied to the Klein-Gordon Equation, the dynamical von-Karman 
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equation and the vibrating beam equation. We concentrate only on the basic features 
those remain invariant when we switch from one problem to another belonging to the 
same class. 

The evolution of a physical system in time is usually described by an initial value 
problem for a differential equation. The differential equation may be an ordinary or a 
partial and mixed initial-boundary value problems are included. Let u(t) describe the 
state of some physical system at the time t. Suppose that the time rate of change of 
u(t) is given by some function A of the state of the system u{t) and the initial data 
m( 0) = / is also given. Then we have the initial value problem 

dfii 

— {t) = Au{t), t>Q, u{0) = f. (1.1.1) 

The function u takes value in a Banach space X and A be an operator from its domain 
D{A) in X to X. 

Abstract Formulation 

Now, we give a brief idea of converting a partial differential equation to an abstract 
form where the space variable is suppressed and the equation looks like an ordinary 
differential equation in time variable in an infinite dimensional space. 

Example 1.1.1 Consider the Laplacian A = ^ in a bounded domain C R” 

OX^ 

with smooth boundary dVl. Consider the initial boundary value problem for the heat 
equation 


^XD 

— = Ary, e (0,oo) x f2, 
w{0,x) = f{x), X E n, 
w(t, x) = 0,t, X E [0, oo) X dQ. 


( 1 . 1 . 2 ) 
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We set 


u{t) = 


regarded as a function of x and take X to be the space of functions on Q e.g., 
for p > 1 or C(Q). The derivatives 

du , dw 

* aF 

are both limits of the difference quotient 

vj{t + h,x) — w{t, x) 
h ’ 

first limit being in the sense of the norm of X and the second limit being a pointwise 
one. We may formally identify ^ with 

To define A we take X = L^Q) and let D(A) = H^{Q) n Define 

Av = Av for V € D{A). 

Equation (1.1.2) are thus written in the abstract form (1.1.1). 

Remark 1.1.2 The boundary conditions of (1.1.2) is absorbed into the domain of 
definition of operator A and into the requirement that u{t) G D{A) for all t>0. 


Example 1.1.3 Consider the initial value problem for the wave equation in R'*' 


= Aw, X G R , t > 0, 
w(Q,x) = fi{x), r G R^, 
^(0,a:) = /2(x), xeR^. 


(1.1.3) 


Let X = X L^(R^) and the domain of A, D{A) = iJ^(R^) x We set 


^ ih) (a 0 
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i.e. 



Then the equation (1.1.3) is converted into (1.1.1). 

Thus, we see that the evolution equations in some sense are the initial value problems for 
the ordinary differential equations in the infinite dimensional spaces and are associated 
with the partial differential equations governing certain physical phenomena. 

In this dissertation we have considered the following type of evolution equations. 

• First Order Evolution Equations. 

(a) Implicit Evolution Equations. 

(b) Quasi-linear integro-differential Equations. 

• Second Order Evolution Equations. 

(a) semi-linear Evolution Equations. 

(b) semi-linear integro-differential Equations. 

The Sobolev type evolution equations are impffcit evolution equations in which the time 
derivatives of the unknown function are not given explicitly. These type of equations 
arise in the study of various diffusion and fluid flow models. Consider an initial boundary 
value problem on a cylinder of the form 

Dt{u + b{x)Lu) + Lu = (1-1-4) 

where Lu is a nonlinear elliptic operator in the divergence form. The coefficient b{.) is 
assumed to be bounded, measurable and nonnegative. The above equation is parabolic 
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where b{x) = 0 and of the Sobolev type where b{x) > 0. The coefficient 6(.) represents a 
quantity with the dimensions of the viscosity. Applications of these and related results 
are considered in [86]. We may write (1.1.4) in an abstract form by choosing a proper 
space and an operator as 

j^{u{t)-^g{t,u{t))) + Au{t) = f{t,u{t)). (1.1.5) 


A partial neutral functional differential equation with unbounded delay is a par- 
ticular case of the nonlinear Sobolev type evolution equation of the form 


— (u(t) + G{t, Ut)) = Au{t) -t- F{t, ut), t > 0, 


( 1 . 1 . 6 ) 


in a Banach space X where A is the infinitesimal generator of an analytic semigroup 
in A”, F and G are appropriate nonlineax functions from [0,r] x W into X and for any 
function u G C((— oo, oo), X), the history function Ut G C((— oo, 0], A) of u is given by 
Ut{d) = u{t+9). The above equation is called the abstract neutral functional differential 
equation (ANFDE) with unbounded delay. 

As a motivational example for this class of equations we consider the boundary 
value problem 



(t,a:)-i- / / b{s — t,ri,x)u{s,r])drids 

J — OO J 0 


+ / a{s — t)u{s,x)ds + ai{t,x), 

J — OO 


u{t, 0) = u{t, tt) = 0, t >0, 
u{6,x) = (l){d,x), d<0, 0 < a: < TT, 


= x) + ao{x)u{t, x), 

t>0, 0 <X < TT, 


where the functions uq, a, oi, bi and (f satisfy appropriate conditions. 
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The partial integro-difFerential equations with infinite delays arising, for example, 
in the study of heat conduction in the materials with memory or population dynamics 
for spatially distributed population, may be described in the form 



(t) — Au{t) + F{t, Ut), 


t > 0, 


( 1 . 1 . 7 ) 


where A is the infinitesimal generator of a strongly continuous semigroup of linear oper- 
ators on a Banach space X . The above equation is called an abstract retarded functional 
differential equation (ARFDE). Some ARFDEs can be conveniently transformed into 
ANFDEs. For example, consider equation (1.1.7) with 

F{ip)= f C{-0)ip{e)d9, 

J — oo 


where C{s) is a strongly continuous map of continuous operator from X into X. As- 
suming that we can decompose C(s) = L{s) -f N{s) where L{s) and N{s) are also 
strongly continuous maps of continuous operators and further L(s) is linear. We define 
the operator V (t) by 


V{t)u 


L{s)uds. 


Jo 

Proceeding formally, the equation (1.1.7) may be transformed into an ANFDE as 


£ 

dt 



V{t — s)u{s)ds 


Au(t) + f N{t — s)u{s)ds, 

J — OO 


which has the form (1.1.6) and in some cases, depending on V and N, it is easier to 
treat this equation than the original one. 


The equation (1.1.5) with a linear g arises in, for example, in thermodynamics, in 
the flow of a fluid through Assured rocks, in the shear in second-order fluids and soil 
mechanics. Several applications of the results for the Sobolev type evolution equations 
to physical problems provide us a motivation to the establish the convergence results 
for the approximate solutions of these equations. 
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The behavior of many physical systems at an instant depends not only the instan- 
taneous values of the functions appearing in the governing equations but their values at 
previous instants also. Such a system is called a system with a memory. For instance, in 
nuclear dynamics and thermoelasticity there is a need to reflect the effects of the “mem- 
ory” of the system. This results in the inclusion of an integral term into the governing 
partial differential equation yielding a partial integro-differential equation (PIDE). 

As an example, Consider the heat conduction equation for a one dimensional rod 
of unit length. In this model the temperature, u{x,t), satisfies 

DtU{x,t) = Da:xU{x,t), 

where Dt is the first partial derivative with respect to time and D^x is the second 
derivative with respect to x, the space variable. In [76] Nunziato observes that the above 
system has two major shortcomings. First, this equation fails to take into account the 
memory effects present in some materials, particularly at low temperatures. Secondly 
this equation implies that a thermal disturbance at any point in the body will be felt 
instantaneously throughout the body. Several alternate formulations which include an 
integral term with respect to time have been considered in [18, 29, 71, 75, 76]. For 
example, [29] considers the following heat flow equation 

Dtu{x, t)= h{t — s)Dx{cr{DxU{x, s))ds + k{x, t), 

Jo 

where the fimctions h, a and k satisfy some appropriate conditions. This provides us 
a motivation to study the partial integro-differential equations. In this work we have 
studied certain classes of first and second order integro-differential equations. 

In many practical applications, the coefficients of the partial differential operators 
depend also on the time and the unknown solution of the problem. Therefore associated 
operator in the abstract Cauchy problem involve time and the unknown function. We 
have considered a class of first order quasi-linear integro-differential equations. 
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For proving the existence and uniqueness results, we have used the theory of 
monotone operators, the theory of semigroups of operators and the method of semi- 
discretization in time which are very efficient tools in the study of nonlinear evolution 
equations. The method of semi-discretization in time is also know as the Rothe method 
or the method of lines. In this method using the time discretization, an evolution equa- 
tion is approximated by the corresponding system of elliptic equations by means of 
which an approximate solution for the original evolution equation is constructed. After 
proving some apriori estimates for the approximate solution, the convergence of the 
approximate solution to the unique solution of the evolution problem is established. In 
these studies either global Lipschitz conditions or local Lipschitz conditions with some 
growth conditions on nonlinear forcing terms have been assumed. 

The Rothe method has its significance both as a numerical method as well as a 
theoretical tool. The existence theorems are proved in a constructive way. Thus, no 
other methods are needed to give preliminary information on existence or regularity of 
the solutions as required in many other numerical methods when question on conver- 
gence or order of convergence etc. are to be answered. The Rothe method is a stable 
method. To the solution of elliptic problem generated by this method, the variational 
method can be applied. As the theoretical results are concerned, they are obtained in 
a relatively simple way. Moreover, the Rothe method, being a very natural one, makes 
it possible to get a particularly good insight into the structure of the solutions. Often 
a brief inspection of the corresponding elhptic problems gives an information what can 
be expected concerning the properties of the solution. 

For the convergence results we are mainly concerned with the approximation of the 
solution to an evolution equation in a separable Hilbert space H by the Faedo-Galerkin 
method and our approach is that of numerical and analytical approximations. Let {uj} 
be orthonormal basis for Hilbert space. Then it is natural to consider the subspace 
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Hn = span{ui,- • • , Un} and the projection Pn : H defined by 

for X = XiUi. In our approach we consider the integral equations associated to the 
evolution equation. With the help of an associated integral equation and the projection 
operators we get an approximate integral equation which has a unique solution. Then 
we establish the convergence of a sequence of the approximate integral equations to 
the associated integral equation and show the limit of the solutions of the approximate 
integral equations be the solution of the associated integral equation. Then we consider 
the projections onto the finite dimensional subspaces Hn of the approximate integral 
equations. Finally, we consider the Faedo-Galerkin approximations and show that it 
converges to the solution of the evolution equation. 


1.2 Literature Survey 


If in an evolution equation, the coefficient of the highest order derivative cannot be so 
isolated, the equation is said to be of Sobolev type. For example the equation 

{A{.)u)tt + {B{.)u)t + C{.)u = g 

is of the Sobolev type. This type of problem called singular if at least one of the 
operator coefiicient tends to infinity in some sense as t — »■ 0 and degenerate if some 
operator coefficient tends to zero as t 0 in such a way as to change the type of 
problem. In 1976 Carroll and Showalter has considered a subclass of singular and 
degenerate Sobolev equations in their book [26] and considered the well-posedness of 
these problems. They have considered the abstract Sobolev type equations 

|(M(«)) + L(u) = /, 


( 1 . 2 . 1 ) 

( 1 . 2 . 2 ) 
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and have used the spectral and the energy methods to analyze (1.2.2) when C is the 
identity operator and B and A are (possibly degenerate) operator polynomials involving 
a closed densely defined self adjoint operator. The equation (1.2.1) is called strongly 
regular when M is invertible and is continuous. Similarly, (1.2.2) is strongly 

regular when C~^B and C~^A are continuous. The strongly regular case is studied in 
[26] and the operators are permitted to be time dependent and nonlinear. They also 
have given well-posedness results for the linear weakly regular equations in which only 
the leading operators are invertible, by means of the classical generation theory of linear 
semigroups in the Hilbert spaces. 

Brill [25] has studied the abstract Cauchy problem for the semi-linear Sobolev type 
evolution equation 

d 

+ Au{t) = g{u{t),t), (1.2.3) 

in which B and A are the linear operators with their domains contained in a real Banach 
space X and the ranges contained in a real Banach space Y. Function g maps X x [0, T], 
T > 0, into Y. He has developed an existence theory for (1.2.3) under the fundamental 
requirement that the operator B dominates the operator A : D{B) C D{A). The 
abstract problem, he has considered arises as a realization of an initial boundary value 
problem for the pseudo-paraboUc partial differential equation 

dMu{x, t) ^ (D°‘u{x, t))|a|<2m-i)- (1.2.4) 

Here, M and £. are the linear elliptic differential operators of orders 2m and 2Z, re- 
spectively, with m>l. The equation (1.2.4) appears in a variety of physical problems, 
for example, in the thermodynamics [27], in the flow of fluid through fissured rocks 
[19], in the shear in second-order fluids [94] and in the soil mechanics [93]. The case of 
a function g which satisfies a local Lipschitz condition has been treated by Showalter 
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[85, 86]. His results, which are proven without using compactness properties, apply only 
to a small class of nonlinearities. 

Lightbourne and Rankin [67] have established an existence result for the Cauchy 
problem for the partial functional differential equation of the Sobolev type 

{Bu{t)y + Au{t) = g{t,ut), t>0, (1.2.5) 

u{t) = ^(t), -r <t <0. (1.2.6) 

The function g is continuous and A and B are closed linear operators with their domains 
contained in a Banach space X and ranges in a Banach spaces Y. The approach taken 
in [67] is to consider the related integral equation 

v{t) = T(t)B(}){0)+ [ T{t - s)g{s,B~^Vs)ds, t>0, 

Jo 

v(t) = B4){t), —r<t<0, 

where T(t), t > 0 is the semigroup of bounded linear operators generated by —AB~^. 

Hernandez and Henfiquez have established in [40] the existence of a mild and a 
strong solution while in [41] the result of the existence of a periodic solution for the 
quasi-linear partial neutral functional differential equation 

+ F{t, xt)) = Ax(t) + G(t, xt), 
at 

where A is the infinitesimal generator of a strongly continuous semigroup of the linear 
operators on a Banach space and F and G are the appropriate functions defined on a 
phase space. 

In [39] Hernandez have proved the existence of a regular solution for the quasi- 
linear evolution equation of the Sobolev type 

iixit) + 9{i, ^(t)) = + /(^> 

at 
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where A is the infinitesimal generator of an analytic semigroup of bounded linear oper- 
ators defined on a Banach space and the functions / and g are continuous. 

To solve an infinite dimensional equation, we adopt the most natural and ele- 
mentary approach of replacing it by an appropriate finite dimensional problem. We 
then solve that approximate problem and prove the convergence to the solution of the 
original problem. For this, we consider a procedure by using the projections onto a 
finite dimensional space. This is referred as Galerkin approximations although every 
slight variant has its own name in literature e.g. ‘Bubnov-Galerkin’, ‘Petrov-Galerkin’, 
‘Faedo-Galerkin’ etc. The Galerkin approximations used for the nonlinear equations 
are called the Faedo-Galerkin approximations. 

Despite the widespread use of the Faedo-Galerkin method (in many applications 
it is referred to as the method of harmonic balance), the convergence behavior in many 
cases is not known. Bazley [20] has considered the IVP for a nonlinear wave equation 
in a separable Hilbert space H 

(Pu 

-|- Au -f M (u) = 0, 

u(0) = (j>, u'(0) = (1.2.7) 

where A : D{A) —^Hhe strictly positive, hnear, self-adjoint operator and (j) € D{A), 
Ip G D{A^) and M ; D{A^) H he & single valued nonlinear map. He has proved the 
convergence of an approximate solution to the solution of (1.2.7). He has first introduced 
the nonlinearities M which are “reproducing” relative to a complete orthonormal se- 
quence {uj}?® and obtained a finite system of the Faedo-Galerkin approximating ODE. 
Then he has shown that if {ui}'^ are eigen functions of A, then the solution of the 
Faedo-Galerkin approximating ODE converges to that of (1.2.7) under some assump- 
tions. His procedure may be considered as generahzed separation of variables for a 
class of nonlinear wave equations. It can be extended to the inhomogeneous equation 
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u" + Au + M{u) — /, as well as the parabolic equation u' + Au + M{u) = f. In [21] , he 
has proved the global convergence of the Faedo-Galerkin approximations to nonlinear 
wave equations. 

Miletta [72] have proved the convergence of the Faedo-Galerkin approximations to 
the solution of an abstract semi-linear evolution of the form 

^ + Au(t) = M(u(t)), 

u{0) = 0, 

in a separable Hilbert space H where A is closed, positive definite, self-adjoint, linear 
operator with domain D(A) dense in H. Also assumed is that A has a pure point 
spectrum 

0 < Ao < Ai < .... 

and a corresponding complete orthonormal system {uj} so that 

Axii and 

where = 1 if i = j and zero otherwise. 

These assumptions on A guarantee that —A generates an analytic semigroup 
M is a continuous nonlinear operator from D{A°‘) — >■ H for some a, 0 < a < 1, which 
satisfies the following Lipschitz condition on the balls in D{A°‘): For each c > 0 there 
exists a constant K{c) s.t., 

(a) ||M(u)|| < K{c) for u € D{A°‘) with ||A.“'u|| < c, 

(b) \\M{ui) - M(u 2 )\\ < K{c)\\A°‘{ui -'U 2 )|| for Ui e D{A°‘) with ||A“Ui|| < c for 
i = l,2. 

Extending the technique of Miletta [72], Bahuguna, Srivastava and Singh [16] have 
proved the global convergence of the Faedo-Galerkin approximation to the solution of 
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an integro-differential equation of the form 

—^ + Au{t) = f{t,u{t)) + J a{t- s)g{s,u{s))ds, 
u{0) = (f), 

where A satisfy the same assumption but nonlinear maps f and g are defined from 
[0,T] X Xo:{T) into H, Xa{T) = C{[0,T],D{A°')) for 0 < a < 1 endowed with the 
supremum norm and the kernel a is in IA{0,T) for some 1 < p < oo. 

Bahuguna [4, 5] have established the existence, uniqueness, continuation of solu- 
tions to the maximal interval of existence, and the global existence of a strong and 
classical solution of the abstract second order semi-linear equation, 

(^) ^ ^ 

u{0) = xq, u'{0) = rci, (1.2.8) 

under the suitable conditions. 

Problem (1.2.8) for the particular case L = — A, the n-dimensional Laplacian, is 
dealt with by Duvaut and Lions [31] and Glowinski, Lions and Tremolieres [34] in the 
context of some viscoelastic materials using the method of variational inequalities. 

Sandefur [84] has considered (1.2.8) for a special case where the operators A and 
B are such that the equation can be written in the following factorized form, 

— {Ai + A 2 ) + A 1 A 2 U = t e (0,T], 

u(0) = a:o, ^(0) = ^1- (1.2.9) 

Using the theory of semigroups together with the method of successive approximations, 
Sendefur [84] applied the results to various physical problems, for instance, the damped 



1.2 Literature Survey 


17 


and strongly damped semi-linear equations, the telegraph equation and the equation of 
motion for a thin panel. 

Aviles and Sandefur [3] have modified the techniques of Sandefur [84] by changing 
the order of integration in the definition of the mild solution to (1.2.9) and applied 
the results to the Klein-Gordon equation, the von Karman equation and the vibrating 
beam equation. 

Webb [95] and Ang and Dinh [2] have considered the following initial boundary 
value problem. 


utt - XAut + f{u) = 0, (x, t) G f) X (0, T), A > 0, 

u = 0, (x,£) G dn X [O,!"), u{x,0) = wo{x), ut{x,0) = wi{x), (1.2.10) 

where is a bounded domain in R'^ with sufficiently smooth boundary dO,. For 1 < 
N < 2, (1.2.10) governs the motion of a Hnear Kelvin solid (a bar if A = 1, a plate 
if A = 2) subject to certain nonlinear elastic constraints. Webb [95] has studied the 
existence and the asymptotic behavior of a solution to (1.2.10) for A = 1,2,3. Ang 
and Dinh [2] have generalized the results of Webb [95] for (1.2.10) using the method of 
successive linearizations. 

Engler, Neubrander and Sandefur [32] have used the theory of semigroups to prove 
the existence and uniqueness of a mild solution of (1.2.8). Balachandran and Park 
[17] has also proved the existence of a mild solution by using the theory of strongly 
continuous cosine families and the Schaefer fixed point theorem. 

The Rothe Method, introduced by E. Rothe [83] in 1930, is an efficient theoretical 
tool for solving a broad scale of evolution equations. Unlike some abstract method for 
the analysis of the existence and uniqueness of a solution for an evolution equation, 
the Rothe method has a strong numerical aspects. The problem solved originally by 
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E. Rothe was very simple one. However, his method turned out to be a very efficient 
tool for obtaining the solution of substantially more complicated evolution problems, for 
example, linear and quasi-lineax parabolic problems of the second order in n dimensions, 
a parabolic problem of arbitrary order, nonlinear problems, hyperbolic problems, the 
Stephen problem, integro-differential problems, mixed parabolic-hyperbolic problems, 
etc. 


The Rothe method has been used and developed by many authors, e.g. O.A. 
Ladyzenskaja [66], T.D. Ventcel, A. M. Iljin, A. S. Kalasnikov, 0. A. Olejnik, S.I. 
Ibragimov, P.S. Mosolov and K. Rektorys [80, 82] in linear and quasi-linear parabolic 
problems. Nonlinear and abstract parabolic problems has been studied by J. Kacur 
[44, 45, 47], J. Necaa [74], A. G. Kartsatos, M. E. Parrott [51, 52] etc. linear and 
quasi-linear hyperbolic equations has been considered by J. Jerome, E. Martensen, M. 
Pultar [79], J. Kacur [46] etc. Many of the obtained results by the Rothe method were 
obtained as well by the other methods like the method of compactness, the theory of 
semigroups, the method of monotone operators, the Fourier transform, etc. 

Two powerful theoretical tools, among many methods, are the theory of semigroups 
and theory of monotone operators for proving the existence and uniqueness results for a 
differential equation in an abstract space. The theory of semigroups of linear operators 
developed quite rapidly after the establishment of the theorem of Hille and Yosida in 
1948 and by now it is an interesting mathematical subject with substantial applications 
to many fields, in particular, to partial differential equations. The books by Hille 
and Phillips [43], Ladas and Lakshmikantam [65], Kato [56], Pazy [77], Martin [70], 
Showalter [90] among several others are extensive treatises on the subject of semigroups 
of linear operators. One of the important results in the theory of semigroups of bounded 
linear operators is due to R.A. Phillips [78] for the case of a Hilbert space and its 
extension in a general Banach space is due to Lumer and Phillips [68]. This result 
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characterizes the infinitesimal generator of a contraction semigroup of bounded linear 
operators as an m-dissipative operator (see Pazy [77], p. 13-14). 

The theory of monotone operators with the Rothe method has been incorporated 
for the first time by Necas [74]. Necas has studied the following nonlinear equation in 
a real Hilbert space H: 

d/ii 

— {t)+Au{t) = f{t), 0<t<T, u{0) = uq, (1.2.11) 

where u is defined on [0, T] into H, A is & nonlinear m-monotone operator defined on 
the domain D{A) C H of A into H and /, defined on [0,T] into H, is a continuous 
function of bounded variation. 

Kartsatos and Zigler [53] have extended the results of Necas for the semi-linear 
problem 

dv 

— (t) + Au{t) = G{t, u{t)), 0 <t <T, 

at 

u(0) = uo, (1.2.12) 

considered in a Banach space X whose dual space is uniformly convex where A defined 
on the domain D{A) C X into X is m-monotone and G defined on [0, T] x X into X 
satisfies a Lipschitz-like condition of the type 

(Al) |lG(t,a:)-G(s,y)|| < ||(^(t) - (/)(s)|| + L||x - 1/|| 

with 0 : [0, T] — > X a continuous function of bounded variation and L > 0 a constant. 

In [47] Kacur has applied Rothe ’s Method to the nonlinear “implicit” integro- 
differential equation 

^{t)+Au{t) = G{t,K{u){t)), a.e. t€{0,T), 
u = (j) on [—q,T], q>0. 


(1.2.13) 
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where u is defined on [0, T] into a Hilbert space if, A is a nonlinear operator defined on 
a reflexive Banach space V into its dual V*, the nonlinear map G is defined on [0, T]xH 
into H and if is a Volterra operator defined on Loo{[-q,T],H) into Loo{[-q,T], H). 
Kacur has proved the existence of a unique strong solution under the following hypothe- 
ses for a given Lipschitz continuous function (j) : [—q, T] H. 

(A2) A is a coercive maximal monotone operator. Coercivity is assumed in the form 

<Au,u> > \\u\\P{\\u\\) - C^\u\^ - C 2 , 

for all u in V where |1.1| and |.| are the respective norm of V and H, < represents 

the usual duality product between V* and V and P{s) — >■ oo as s oo. 

(A3) G satisfies a Lipschitz condition of the type 

lC?(t,u) - G(s,u) I < - s| -h |t- s|lu| + |u - u|], 

for all t, s G [0, T] and all u,v ^ H where C > 0 is a constant. 

(A4) The nonlinear Volterra operator K maps Lip(5'r, if) into Lip(S'r, H) where St’, = 
[—q, T] and satisfies 

(i) \\K{u) - K{v)\\ciST,H) < C\\u-v\\c(Sr,H), 

for all u,v & C{St, H) where G > 0 is a constant; and 

(ii) \K{u){t) - K{u){s)\ < \t- slL(||ti||c(Sr,if)) (l + ll^LoolST.iy)) ’ 
for all s,t e St, s <t, u e Lip(S'r, if) and L : R+ — >• R+ is continuous. 

In the earlier studies on Rothe’s method, we notice that one proves an apriori 
bound for difference quotients without proving any bounds for the discrete 

points Uj, hence required global Lipschitz conditions of the types (Al), (A2) and (A4) 
which has been considered by Kartsatos and Zigler [53] and Kacur [47]. Bahuguna 
modified Rothe’s method which allows to weaken the global Lipschitz condition to 
a local Lipschitz condition plus a growth condition. Using the growth condition, he 
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has proved that the discrete points lie in a ball whose radius is independent of the 
discretization parameters and then using the local Lipschitz condition he has proved 
the bounds for the difference quotients. 

Bahuguna and Raghavendra [8] have considered the abstract integro-differential 
equation 


a{t - s)k{s,u{s))ds f{t), t>0, 

ti(0) = no, (1.2.14) 

in a Banach space X whose dual is uniformly convex where A is a nonlinear single 
valued operator from D{A) C X into X, the mapping k is defined on [0, t] x X into 
X, the X- valued function / and the real- valued function a are defined on [0,T]. He 
has proved, using Rothe’s method, the existence of a strong solution of (1.2.14) under 
the assumption that the operator A is m-monotone such that (/ + A)~^ is compact, 
mapping k is continuous in both the variables and satisfies the growth condition 

||/c(f,a;)l| < Ci||x||+C2, 

for all t E [0,r] and a: E X, where Ci, C 2 > 0 are the constants, the functions / 
and a are the continuous functions with bounded variations on [0,T]. The uniqueness 
is proved under an additional condition of Lipschitz continuity on k{t, x) in x, a.e. 
te[o,T]. 

Bahuguna [9] has also used Rothe’s method to establish the existence of a unique 
strong solution of the nonlinear implicit abstract integro-differential equation 

rill 

^{t) + Auit) = f{t,uit),F{u){t)), fora.e. tE{0,T), 
at 

u = (f) on [-q, 0], </> E Lip{[-q, 0],H), 




(1.2.15) 
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in a Hilbert space H. 

The initial work using the theory of m-accretive operators to the quasi-linear evo- 
lution equations is due to Kartsatos [48]. Kartsatos [49] and Kartsatos and Parrot [50] 
have established the existence and uniqueness and used the Galerkin method for the 
approximation of solutions to the quasi-linear functional problem 

vi!{t) A{t,ut)u{t) = 0, te[0, T),, ■u(O) = 0, 

in a Banach space X, where A{t,(j))v is m-accretive in v E D for every pair E 
[0, T) X Co, D is a subset of X, Co is a closed subset of the space of all continuous 
functions (j) on (— oo,0] into D with ll<^||oo < f for some fixed r > 0 and Ut E Co, given 
by Ut{s) = u{t -b s), s E (— oo, 0]. Murphy [73] has constructed a family of approximate 
solutions to the homogeneous quasi-hnear evolution equation 

u'{t) + A{t,u{t))u{t) =0, tG(0,T], 7jl{0) = (p. (1.2.16) 

He has shown that the approximate solutions converge to a limit solution and this limit 
solution becomes a unique solution. Kartsatos [48] established a theorem concerning 
the existence of a unique strong solution to (1.2.16) under the assumption that A(t, u)v 
is Lipschitzian in t, u and m-accretive in v. In this paper, conditions on A{t, u)v is 
motivated by Kartsatos [48]. 

For the study of the abstract quasi-linear evolution equation 

u'{t) + A{t, u{t))u{t) = f{t, u{t)), t E (0, T], u(0) = 0, (1.2.17) 

we refer to T. Kato [61, 58, 60, 62], S. Kato [54], Amann [1] and references cited therein. 
The crucial assumption in these papers is that the linear operator A(t,u), depending 
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on t and the unknown u, has the dual property of being the negative generator of a 
Co-semigroup but not necessarily of an analytic semigroup and at the same time, a 
bounded operator on Y into X. Moreover, in [61, 58, 60, 62] there is an isomorphism 
C of y onto X with the property that 

SA{t, w)S~^ = A{t, w) + B{t, w), 

where B{t,w) is in the space B{X) of all bounded linear operators from X into X 
whereas in [54] it is assumed that there is a family {S{t,w)} of the isomorphisms of Y 
onto X such that 


S{t,w)A(t,w)S{t,w) ^ = A{t,w) + B{t,w), 
where B{t, w) € B{X) and S{t,w) satisfies the Lipschitz-like condition 
|15(t,u;i) - 5(s,W2)lly.x < f^i\t - s\ + \\wi - w^Wx) , 

where /x is a constant. Kato [61] has proved the existence, uniqueness and continuous 
dependence on the initial data of the solution of an abstract quasi-linear and has shown 
that these results are applicable to the different kinds of quasi-linear equations such as 
the symmetric hyperbolic systems of the first order, the wave equation, the Korteweg- 
de Vries equation, the Navier-Stokes and the Euler equation, a magnetohydrodynamics 
equation, the coupled Maxwell and the Dirac equations, etc. Kato [54] has proved the 
existence of a strong solution to (1.2.17) under the conditions on A{t, u) and /(t, u) for 
(t, u) e JxW similar to that of Crandal and Sougandis [28]. Amann [1] has treated the 
various cases of (1.2.17) in interpolation spaces using the theory of analytic semigroups. 

Bahuguna [4] has also proved the similar results by using Rothe’s method for the 
abstract quasi-linear explicit integro-differential equation 

^(t) -f A{u(t))u{t) = [ a{t- s)k{s, u{s))ds + f{t), 
dt Jq 
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for 0 < t < T with the initial value Uq in real reflexive Banach space X whose dual is 
uniformly convex, A(u) is a linear operator in X for each u in an open subset W of 
Y where Y is also a real reflexive Banach space which is continuously and compactly 
embedded in X, a is a real- valued and / is a F- valued functions defined on [0,T] 
and k is a. Y- valued map defined on [0, T] x W. Using the semigroup theory together 
with the contraction mapping theorem Bahuguna [5] has established the existence of 
a unique classical solution to quasi-linear explicit integro-differential equation in the 
Banach space X, 

du 

— {t) + A(t, u{t))u{t) = J a{t- s)k{s, u{s))ds + f{t), 
u(0) = Uq, 

where A{t, u) is a linear operator in X for each u in an open subset W of Y, maps a, 
/, k and spaces X, Y have same properties as in [4]. Also using similar techniques of 
paper [4] , we [10] have established the existence, uniqueness and continuous dependence 
of a strong solution to the quasi-linear “implicit” integro-differential equations 

du 

— (t) -i- A{u{t))u{t) = fit, u{t),Giu)it)), 0 <t<T, u(0) = uo, 

at 

in a reflexive Banach space whose dual is uniformly convex. 

1.3 Preliminaries 

In this section, we mention some of the results required for the analysis in the ensuing 
chapters. We mostly mention these results without proofs with proper references. 
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1.3.1 Semigroups 

Definition 1.3.1 Let .X" be a Banach space. A one parameter family T{t), 0 <t < oo, 
of bounded linear operators from X into X is a semigroup of bounded linear operator 
on X if 

(i) T(0) = J, where / is the identity operator on X. 

(ii) T{t + s) = T{t)T{s) for every t,s > 0 (the semigroup property). 


The linear operator A defined by 

D{A)-- 

and 


T(t)x — X 

xeX:]im-^ exists 


tio 


t 


} 


T{t)x — X d'^T{t)x 
Ax = hm ^ ' 


tiO 


dt 


for x 6 D{A) 


t=o 


is the infinitesimal generator of the semigroup T{t), D{A) is the domain of A. 


Example 1.3.2 Let T{t) : R — > R be given by 

T{t) = e^^Xo, t >0, xo G R. 

Then T{t) is a semigroup of bounded linear operators on R and its infinitesimal gener- 
ator is A = XI. 


Example 1.3.3 Let T{t) : R"' — > R'^, be such that 

T(t) = e^^xo, t > 0, xo G R'^, 

where B = [6ij] is an n x n matrix, bij are the constants. Then T{f) is a semigroup of 
bounded linear operators on X = R” and its generator is A = B. 
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Example 1.3.4 Let T(t) : X — > X be given by 

T(£) = e^^xo, t>0, xo € R, 

where B : X — X is a bounded linear operator. Then T{t) is a semigroup of bounded 
linear operators on X. and its generator is A = B. 

Definition 1.3.5 The semigroup of bounded linear operators, T{t) is uniformly con- 
tinuous if 

lim||r(()-/|l= 0 . 

We note that the semigroup T(£) in each of the above examples is uniformly con- 
tinuous. In fact, we have the following result. 

Theorem 1.3.6 A linear operator A is the infinitesimal generator of a uniformly con- 
tinuous semigroup if and only if A is a bounded linear operator. Furthermore, the 
semigroup T{t) generated by a bounded linear operator A : X — >• X is unique and 

(a) There exists w > 0 such that llT(t)ll < e“h 

(b) The map t Tit) is differentiable in norm and 

= AT{t) = Tit) A. 

For applications to partial differential equations, we can not expect the infinites- 
imal generator of a semigroup to be a bounded linear operator. Therefore, the class 
of uniformly continuous semigroups is of little use for this purpose. Thus, we need to 
consider other type of semigroups. 

Definition 1.3.7 A semigroup Tit), 0 < £ < oo, of bounded linear operators on X is 
a strongly continuous semigroup of bounded hnear operator if 

lim Tit)x = X for every x e X. 

tio 
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A strongly continuous semigroup of bounded linear operators on X will be called a 
semigroup of class Co or simply a Co-semigroup. 

Example 1.3.8 In the Banach space X = (7(0,1] of continuous functions with sup- 
norm, define the family of operators {T(t)}, t > 0 by the formula 

r(t)r(0 = x[^/(i + t0], x 6A, ee[o,i]. 

Then {T{f)}, t> 0, is a (7o-semigroup of operators in X. 


If T{t) is a (7o-semigroup, then there are constants w > 0 and M > 1 such that 

||T(t)|| < for 0 < t < oo. 

Also, for every x € X, the map t T{t)x is a continuous function from Rg (the 
nonnegative real line) into X (cf. §1.2 in Pazy [77]). 

We also have the following main result for (7o-semigroups. 


Theorem 1.3.9 Let T{t) be a Co-semigroup and let A be its infinitesimal generator. 
Then 

(a) For X e X, 

lim — / T(s)xds = T{t)x. 

h — ►O fx j 

(b) For X £ X, f^T(s)xds € L)(A) and 

A T{s)xd^ = T{t)x — X. 

(c) For X € D{A), T(t)x € D{A) and 


(d) For X € D{A), 


^T{t)x = AT{t)x = T{t)Ax. 
at 


T(t)x - T{s)x = J T{T)AxdT — AT{T)dT. 
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(cf. Theorem 1.2.4 in Pazy [77].) 

Let T{t) be a Co-semigroup, then there are constants w > 0 and M > 1 such that 
||r(t)ll < for t > 0. If a; = 0, then T{t) is called uniformly bounded semigroup 

and if moreover M = 1, then T(t) is called Cq- semigroup of contractions. 

Example 1.3.10 In the Banach space X = C[0, oo) of continuous, bounded functions 
on [0,oo) with sup-norm, define the family of translation operators {T{t)}, t > 0, by 
the formula 

T{t)x{^) = x{t + xeX, ^>0. 

Then {T(t)}, t > 0, is a Co-semigroup of contraction of operators in X. 

For a linear operator A not necessarily bounded, in X, the resolvent set p{A) of A 
is the set of all complex numbers A for which XI — A is invertible, i.e. {XI — A)“^ is a 
bounded linear operator in X. The family R{X : A) = {XI — A)“^, A G p{A) of bounded 
linear operators is called the resolvent of A. 

One of main results in the theory of semigroups of bounded linear operators is the 
Hille-Yosida theorem which was established in 1948 by Hille [42] and Yosida [96]. This 
result establishes a characterization of the infinitesimal generator of a Co-semigroup of 
contraction. 

Theorem 1.3.11 (Hille-Yosida). A linear (unbounded) operator A is the infinitesimal 
generator of a Cq- semigroup of contractions T{t), t>0 if and only if 

(i) A is closed and D(A) is dense in X, 

(ii) the resolvent set p{A) of A contains [0, oo) and for every A > 0 


||iJO:/l)|l < i. 
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We note that if T{t) is a Co-semigroup with A its infinitesimal generator such that 
for some uj >0, ||T(t)l| < then S(t) = is a Co -semigroup of contractions 

and its infinitesimal generator is A — ojI. 

The Hille-Yosida theorem gives us a characterization of the infinitesimal generator 
of a Co-semigroup of contractions. Another characterization of such operators is given 
by the burner- Phillips theorem which we state below. 

Now, we consider the notion of the dissipative operators. To define the dissipative 
operators we require duality map defined from X into the dual space X* of X, given 

by 


F(x) = {x’€X' ■.<x\x>= IWI“=I|i’K}. 

where < ., . > denotes the duality paring between the dual space X* of X and X. 
Clearly, duality map is a multi-valued or set valued map in general. Prom the Hahn- 
Banach theorem it follows that F{x) ^ 0 for every x E X. In particular, when X is a 
Hilbert space then the duality map is single- valued. 

Definition 1.3.12 A linear operator A is called dissipative if for every x E D{A) there 
exists ax* E F{x) such that 

Re < x*,Ax > < 0. 

A dissipative operator A for which R{I — a A) = X, for every a: > 0 is called 
m-dissipative. If operator A is dissipative the —A is accretive or monotone operator. 


The following theorem gives the characterization of the dissipative operators. 



1.3 Preliminaries 


30 


Theorem 1.3.13 A linear operator A is dissipative if and only if 

A||x|| < ll(A/-/l)x|| 

for all X G D{A) and A > 0. 

Theorem 1.3.14 (Lumer and Phillips) Let A be a linear operator with dense domain 
D{A) in X. Then we have the following. 

(a) if A is dissipative and there exists a Xq > 0 such that the range, R{XoI — A) of 
{XqI — A) is X , then A is the infinitesimal generator of a Co-semigroup of contractions 
on X . 

(b) if A is the infinitesimal generator of a Co-semigroup of contractions on X then 
R{XI — A) = X for all X> 0 and A is dissipative. Moreover, for every x G D{A) and 
every x* G F(X), Re < x*,Ax > < 0. 

(cf. Theorems 1.4.2 and 1.4.3, resepectively, in Pazy [77].) 

As pointed out earlier that what we know is the operator A associated with some 
partial differential operator and the problem is to find the semigroup generated by it. 
For a bounded linear operator B, we already know that the semigroup is uniformly 
continuous and given by The main task is to find the representation of the semi- 
group generated by an unbounded linear operator, if there is any. For this purpose we 
use the complex operator calculus. For a bounded linear operator B we have another 
representation given by the following theorem. 

Theorem 1.3.15 Let B be a bounded linear operator. If 'y > ||5||, then 

1 /•7+ioo 

= e^‘(AJ - 5)-^dA. (1.3.1) 

27rz J y — joQ 

The convergence in (1.3.1) is in the uniform operator topology and uniformly in t on 
bounded intervals. 
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(cf. Theorem 1.7.1 in Pazy [77].) 

The following theorem gives the sufficient conditions for an operator A to be the 
infinitesimal generator of a Co-semigroup and a similar representation as in (1.3.1). 

Theorem 1.3.16 Let A be densely defined operator in X satisfying the following con- 
ditions. 

(i) For some 0 <5 < 7r/2, p{A) D = {A : \arg\\ < (7r/2) -h U {0}. 

(a) There exists a constant M such that 

ll(A/-4)-|| < M 

for A G Ei with A > 0. 

Then, A is the infinitesimal generator of a Co-semigroup T{t) satisfying ||T(t)|| < C 
for some positive constant C and 

T{t) = ~ [ e^\XI - A)-'^dX 

2tti jp 

where T is a smooth curve in running from ooe“^^ to coe^^ for n/2 < 6 < 7r/2 -I- 5. 
The integral converges for t > 0 in the uniform operator topology. 

(cf. Theorem 1.7.7 in Pazy [77].) 

We are more interested in the special kind of semigroups known as analytic semi- 
group from the application point of view. Therefore we now consider the possibility of 
extending the domain of the parameter to a region of the complex plane that includes 
the negative real axis. 

Definition 1.3.17 Let A = {2 G C : ipi < axg z < (f 2, (pi < 0 < 1P2} and for 
2 G A let T{z) be a bounded linear operator. The family T{z) , 2 G A is an analytic 


semigroup in A if 
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(i) z -^T{z) is analytic in A; 

(ii) T’(O) = I and lim^_^o T{z)x = x for every x e X\ 

(iii) T{zi + Z 2 ) = T{zi)T{z 2 ) for zi,Z 2 € A. 

A semigroup T{t) will be called analytic if it is analytic in some sector A containing 
the nonnegative real axis. 

If T{t) is a Co-semigroup which can be extended to an analytic semigroup in some 
sector then the semigroup e“*T(t) can also be extended to an analytic semigroup in 
some sector. Therefore we may restrict ourselves in most of the cases to uniformly 
bounded Co-semigroups. For convenience, we suppose that 0 € p{A) where A is the 
infinitesimal generator of the semigroup T{t). This may also be achieved by multiplying 
the uniformly bounded semigroup T(t) by for e > 0. 

We assume that the densely defined operator Ain X satisfies 

^ X] ^ ^1 < I + < 5 } U {0} (1-3.2) 

for some 0 < 5 < |, and 

l|R(A:.4)|| < A for AeJ], A#0. (1.3,3) 

It follows that if A is densely defined operator, A satisfies (1.3.2) and (1.3.3), then 
the semigroup generated by A can be extended to an anal 5 h;ic semigroup in the sector 
A^ = { 2 : ; [arg z\ < 5} and every closed sub-sector Ay = {z : |arg z\ < 5' < <5}, 
1|T(2:)|| is uniformly bounded. In fact, we have the following theorem. 

Theorem 1.3.18 Let T{t) be a uniformly bounded Co-semigroup. Let A be the in- 
finitesimal generator ofT{t) and assume that 0 6 p{A). Then the following are equiva- 


lent. 
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(a) T{t) can be extended to an analytic semigroup in a sector Aj = {^ : \arg z\ < <5} 
and ||T'(t)|| is uniformly hounded in every closed sub-sector A;/ of As, 5' < 5 of As- 

(b) There exists a constant C such that for every a > 0, t ^ 0, 

||B(£r + iT:/!)|| < 

rl 

(c) There exist 0 < 5 < xjl and M > 0 such that 

p(^) D £ = |a : \arg A| < — + (jj- U {0} 

and 

l|fi(A:.4)|| < ^ for AeE, A # 0. 

(d) T{t) is differentiable for £ > 0 and there exists a constant C such that 

prwii < ^. i>o. 

(cf. Theorem 2.5.2 in Pazy [77]). 

1.3.2 Fractional Powers of Operators 

For an operator for which —A is the infinitesimal generator of an analytic semigroup 
T{t), it is possible to define fractional power of A. More precisely, we suppose that A 
is a densely defined closed linear operator satisfying 

P(^) T) = {A : 0 < a; < |arg A| < tt} U V, (1.3.4) 

where V is a neighborhood of zero, and 


l|B(A:A)l| < for AsE^- 


(1.3.5) 
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If M — 1 and uj — tx/ 2 then —A is the infinitesimal generator of a Co-semigroup. 
If w < 7r/2 then, —A is the infinitesimal generator of an analytic semigroup. 

For an operator satisfying (1.3.4) and (1.3.5) and a > 0 we define 

where C is a path in the resolvent set A from ooe"*^ to ooe'®, to <9 <tt, avoiding the 
negative real axis and the origin and z~°‘ is taken to be positive for real positive values 
of 2 . The integral (1.3.6) converges in the uniform operator topology for every a > 0 
and thus defines a bounded linear operator A~°^. 

We observe that for a, /? > 0 


and there exists a constant C such that 


\\A “II < C for 0 < a < 1. 

For A satisfying (1.3.4) and (1.3.5) with oj < 7 r/ 2 , for a > 0, we have 

A“ = (A-“)-^ 


For a = 0, A“ = J. 

For the main application of the analytic semigroups, we require the following the- 
orem. 

Theorem 1.3.19 Let —A be infinitesimal generator of analytic semigroup T{t). If 
0 e p(A), then we have the following. 

(a) T{t) : X — > D(A“) for every t>0 and a > 0. 
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(b) For every X G D{A°‘), we have T{t)A°‘x = A°‘T{t)x. 

(c) For every t > 0 the operator A°‘T{t) is bounded and 

for some constants M^, depending on a, and ^ > 0. 

(d) For 0 < a < 1 and x E D{A°^), 

\\T{t)x-x\\ < CaniA^xW. 


(cf. Theorem 2.6.13 in Pazy [77].) 

Corollary 1.3.20 Let A be the infinitesimal generator of an analytic semigroup T(t). 
If f E Lf{t),T\X) is locally Holder continuous on (0,T], then for every x E X the 
initial value problem 

= Au{t) + f{t), 
ii(0) = X, 


has a unique solution u E (^((OjT] : X) given by 

uft) = T{t)x + f T{t — s)f{s)ds. 

Jo 


(cf. Corollary 4.3.3 in Pazy [77].) 


1.3.3 Evolution Systems 

Definition 1.3.21 A two parameter family A two parameter family of bounded linear 
operators U{t,s), 0 < s < t < T, on X is called an evolution system if the following 
two conditions are satisfied: 
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(i) U(s,s) = I and U{t,r)U(r,s) = U{t,s), 0<s<r<t<T. 

(ii) The map {t,s) U{t,s) is strongly continuous for 0 < s <t <T. 

Definition 1.3.22 Let X be a Banach space. A family {A(i)}tg[o,T'] of infinitesimal 
generators of Co-semigroups St on X is called stable if there are constants M > 1 and 
u) (called the stability constants) such that 

p(A(t)) D]a;, oo[ for te[0,T] 


and 


nfi(A : Afe)) 

J = 1 


< M( A — oj) ® for \> u! 


and every finite sequence 0 < < ^ 2 ) • • • ^T, k = 1,2, 


Let X and Y be two Banach Spaces such that Y is densely and continuously embedded 
in X, i.e., y is a dense subset of X and there is a constant G such that 

for VJ ^Y. 

The subspace T of X is an invariant subspace of a linear operator S : D{S) C X X 
as : D{S)nY ^Y. 

Let St{s), s > 0, be Co-semigroup generated by A(t), t G [0,T]. A subset Y of 
X called Aft) -admissible if Y is an invariant subspace of operator St{s), s > 0 and the 
restriction St(s) of St{s) to T is a Co-semigroup on Y. Moreover, A{f), the part of A(t) 
is, in this case, the infinitesimal generator of the semigroup St{s) on Y. We make the 
following assumptions. 


(HI) {A{t)}tej is a stable family with stability constant M, u. 
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(H2) Y is yl(t)-admissible for t e J and the family {i(f)}te 7 of parts A(t) of A(t) in 
Y, is a stable family in Y with stability constant M, u. 

(H3) For t e J, D{A(t)) D Y , A{t) is a bounded operator from Y into X and t — ^ A{t) 
is continuous in the B{Y,X) norm ||.||y_>x. 

Theorem 1.3.23 Let A{t), 0 <t <T, be the infinitesimal generator of a Co-semigroup 
‘S'f(s), s>0, on X. If the family {A(T')}tg[o_ 7 ’] satisfies the conditions (H1)-(H2), then 
there exists a unique evolution system U{t, s) , 0 < s < t < T , in X satisfying 

(El) ||C/(t,s)|| < for 0 < s < t < T. 


(El) ^U(t,s)w 


= A(s)w for w eY, 0 < s <T. 


t~S 


(E^) -^U{t,s)w = —U(t,s)A{s)w for luEY^ 0 < s <t <T. 


Where the derivative from the right in (E 2 ) and the derivative in (E^) are in the strong 
sense in X . 


(cf. Theorem 5.3.1 in Pazy [77].) 


Theorem 1.3.24 Let satisfy the condition of Theorem (1.3.23) and let 

U{t,s), 0 < s <t <T be the evolution system given in Theorem (1.3.23). If 


(Ei) U{t,s)YcY for 0<s<t<T, 
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(E^) forv G Y, U(t,s)v is continuous in Y forO ^ s ^ t < T 
and f G C([s,T] ; Y), then for every v eY the initial value problem 
^^ = A{t)u{t) + f{t) for 0<s<t<T, 

u{s)^v, (1.3.7) 

possesses a unique Y -valued solution u given by 

u{t) = U{t,s)v j U{t,r)f{r)dr. 

(cf. Theorem 5.5.2 in Pazy [77].) 


Theorem 1.3.25 (Contraction Mapping Theorem) Let (X^d) be a (nonempty) com- 
plete metric space and let F : X ^ X be a function. If there exists 0 < r < 1 such that 
d{F{x), F{y)) < rd{x,y) for all x,y E X, then F has a unique fixed point in X, i.e., 
there exists a unique x EX such that Fx = x. 


Theorem 1.3.26 (Gronwall’s Inequality) Let C be a non-negative constant and Letu 
and V be nonnegative continuous function on some interval to<t < to a satisfying 

u{t) < C+ u{s)v{s)ds, tG[ta, to + a], 

Jto 


then the inequality 


u(t) < C exp 



t E [to, to + o 


holds. 
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Lemma 1.3.27 (cf. Theorem 5.6.7 in Pazy [77]) Let ip{t,s) > 0 be continuous on 
0 < s < t <T. If there are positive constants A, B, a such that 

ip{t, s) < A + B j^{t- s)da for 0<s<t<T, (1.3.8) 

then there is a constant K such that (p{t, s) < AK for 0 < s < t <T. 


Proof. We have a well-known identity 


J {t — t)°‘ ^(r — s)^ = (t — 


r(c. + p) • 


(1.3.9) 


Iterating (1.3.8) n — 1 times using the identity (1.3.9) and estimating t — s by T we find 


Choosing n sufficiently large so that no: > 1 and estimating {t — ^ by we 

<p{t,s) < Ci-fC' 2 ^ (p{a,s)da, 


get 


where 


j=o ^ 

Gronwall’s inequality implies that 


C,=.4E(^T and 

V n/ / r(na) 


rj^na—l 


(p(t,s) < Cie^^^<AK. 


Since Ci and C 2 do not depend on s this estimate hold for 0 < s < t < T. 


Let X and Y are reflexive Banach spaces and Y is continuously and densely, em- 
bedded in X. Let N{X, /?) be set of densely defined linear operators C in X such that 
such that if A > 0 and A/3 < 1, then (/ + AC) is one to one with a bounded inverse 
defined everywhere on X and 


l|(J-hAC)-'|U < (l-A/3)-' 
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where / is identity operator. 


Lemma 1.3.28 Let S : Y X be a linear isometric isomorphism, C € 

Y C D{C), domain of C, P E B{X), the space of all bounded linear operators on X 
and SC = CS + PS. Set 6 = j3 + llPllx- Then for every y E X and A > 0 such that 
Xd < I, the problems 

X + XCx = y 


and 

X -1- X{Cx + Px) = y 

have unique solutions x and x in X. Moreover 


iiiiu < 


i|U < (1-^^) 'Ikiu, 


and if y eY , then x eY and 


\x\\y < {1-X9) lt/||y. 


(cf. Lemma 2 [28].) 

The operator A is called m-monotone if A satisfies the condition 

||ri — u + q;(Au — Ai/)|| > ||m — t>|| for every u,v E D[A), a > 0 (1.3.10) 

and the range R{I + aA) = X ior a > 0. For and m-monotone operator A, let us 
introduce the following sequences operators 

Jn = {I + n~^A)~'^, An = AJn = n{I - Jn) for n = 1, 2, • • • , 

where AJn denotes the composition of the two maps A and Jn- The and An axe 
defined everywhere on X. 
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Lemma 1.3.29 Let A be m-monotonic. JnCLudA^ are uniformly Lipschitz continuous, 
with 

— < ||a; — y||, — A„y|| < 2n||a: — y||, 

where 2n may be replaced by n if X is a Hilbert space. 

Lemma 1.3.30 Let A be m-monotonic. The An are also monotonic. Furthermore, we 
have 

||^n«|| < for u e D{A). 


Lemma 1.3.31 Let X be a Banach space with uniformly convex dual X*. Let A : 
D{A) C X X be m-monotone. Then we have the following. 

(a) Ifun € D{A), n = 1, 2 , • • • , Un—^ u in X and if ll^lwrill bounded then u € D{A) 
and Aun Au. (Here means weak convergence in X.) 

(b) If Xn E. X, n = 1,2, - Xn u in X and if ||AnXn|| are bounded, then u e D{A) 
and AnXn — ^ Au. 

(c) AnU — ^ Au if u E D(A). 


(cf. Lemma 2.2, 2.3 and 2.5, respectively, in Kato [56].) 


1.3.4 Approximate Solutions 


Now consider the abstract Cauchy problem 

^(t) + Au(t) =0, 0<t<T, u{0) = uo, 

dt 


(1.3.11) 


where A : D{A) G X X is m-monotone. Consider the approximate equations 
^^{t) -t- AnUn{t) = 0, Un{0) = Uq. 


dt 


(1.3.12) 
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It is proved by Kato [56] that if uq G D{A) then 

ll“«(t)li < c, ||«'(t)|| < c, 

for all n = 1, 2, • • • where C is a constant independent of n and the strong limit u{t) = 
limn_,oo Wn(i) exists uniformly on [0,T] with u being Lipschitz continuous on [0,T] and 
u(0) = uo- Furthermore, Kato [56] has proved that u{t) G D{A) for t G [0,T] and Au{t) 
is bounded and weakly continuous. 

The main results of Kato [56], we are interested in are the following lemmas. 

Lemma 1.3.32 For each f G X*, < u{t),f > is continuously differentiable on [0,T] 
with 

d 

— <u{t),f>= -<Au{t),f>. 

Proof: Since u„(t) satisfies (1.3.12), we have 

< Un{t),f > = <Uo,f > - [ < AnUn{s), f > ds. 

Jo 

Since Un{t) -4 u{t) and Aun{t) — ^ Au{t) and | < AnUn{t), / > | < C'[|/||, we get 

<u{t),f> = <uo,f>- [ <Au{s),f>ds. ( 1 . 3 . 13 ) 

Jo 

Continuity of the integrand gives the required result. 

1.3.5 Bochner Integral 

Definition 1.3.33 Let (S, B, m) be a measure space and let u be a mapping defined on 
S with values in a Banach space X. u is called weakly B measurable if for any f € X* 
the scalar function s i->-< u{s),f > on 5 is B-measurable. The function u is called 
finitely-valued if it is constant 0 on each of a finite number of disjoint B-measurable 
subsets Bi with m{Bi) < oo and u = 0 on S' - UiBi. A function u is said to be strongly 
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B-measurable if there exists a sequence of finitely-valued functions strongly convergent 
to u m-a.e. on S. 


Definition 1.3.34 The function u is said to be separably-valued if its range {u(s) : 
s E S} is separable, i.e., its range is a closure of a countable set. The function u is 
m-almost separably- valued if there exists a ;B-measurable set Bq of m-measure zero such 
that {u{s) : s E S — Bq} is separable. 


From a result due to Pettis (cf. Yosida [96], p. 131), it follows that u is strongly 
B-measurable if and only if it is weakly B-measurable and m-almost separably- valued. 

Let It be a finitely-valued function defined on a measure space (S, B, m) with values 
in a Banach space X. Let u be equal to Uj 0 on Bj E B, i = 1,2, - • ■ ,n, where Bj’s 
are disjoint and m{Bi) < 00 for i = 1,2, • • • ,n and a; = 0 on (S — Then we 

can define the m-integral J u{s)m{ds) of u over S by 

n 

Y^Uim{Bi). 

i=l 

With the help of a limiting process, we can define the m-integral of more general func- 
tions. More precisely we have the following definition. 


Definition 1.3.35 A function u defined on a measure space (5',B,m) with values in 
a Banach space X is said to be Bochner m-integrable, if there exists a sequence of 
finitely- valued functions {w„} which strongly converges to u m-a.e. in such a way that 

lim / ||u(s) — u„(s)||m(ds) = 0. 

n-*oo Jg 

For any set B E B, the Bochner m-integral of v. over B is defined by 

/ u(s)m(ds) = lim / XB(s)it„(s)m(ds), 

Jb 

where xb is the characteristic function of B, i.e., Xs(s) = 1 if s 6 B and it is equal to 


zero otherwise. 
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Prom a result due to Bochner (cf. Yosida [96], p. 133), we have that a strongly B- 
measurable function u is Bochner m-integrable if and only if ||ii|| is m-integrable. 

Lemma 1.3.36 Au{t) is Bochner integrable andu{t) is an indefinite integral of —Au{t). 
The strong derivative du{t)/dt exists almost everywhere and equals -Au(t). 

Proof: Let Xq be the smallest closed linear subspace of X containing all the values of 
the AnUn{t) for t € [0, T] and n = 1, 2, • • • . Since An.Un{t) are continuous, Xq is separa- 
ble. Since AnUn{t) Au{t) and since Xq is weakly closed it follows that Au[t) G Xq. 
Thus Au(t) is separably-valued. Since it is weakly continuous it is strongly measurable 
(cf. Yosida [96], p.l31) and being bounded, it is Bochner integrable. Prom (1.3.13) 
we have that u{t) is an indefinite integral of -Au{t). For the existence of the strong 
derivative of u{t) almost everywhere on [0,r] and its being equal to —Au(t) almost 
everywhere on [0,!r] we refer to Yosida [96], p. 134. 

(cf. Lemma 4.6 in Kato [56].) 

1.4 An Outline of Dissertation 

In this section we briefly describe the analysis carried out in the ensuing chapters. 

In Chapter 2, we consider the Sobolev type evolution equation in a separable 
Hilbert space and study the approximations of solutions. We first consider an asso- 
ciated integral equation and a sequence of approximate integral equations using the 
projection operators. Then we establish the existence of a unique solution to every 
approximate integral equation using the contraction mapping theorem. Next, after 
proving some estimates for the solution of the approximate integral equation, we prove 
the convergence of the solution of the approximate integral equation to the solution of 
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the associated integral equation. Finally, we consider the Faedo-Galerkin approxima- 
tions of the solution and prove some convergence results. This work is an extension of 
a paper by the author [11] (with Bahuguna). 

In Chapter 3, we consider a strongly damped semi-linear wave equation and re- 
formulated it as a second order semi-linear evolution equation. Then we study the 
approximation of the solution to the second order semi-linear evolution equation in a 
separable Hilbert space by using the similar techniques which are used in Chapter 2. 
First, with the help of a pair of associated integral equations and projection operators, 
we consider the pair of approximate integral equations and prove the existence and 
uniqueness of a solution. We then estabhsh the convergence of the pair of approxi- 
mate integral equations to the pair of associated integral equations and the limit of 
the solution of the pair of approximate integral equations be the solution of the pair of 
associated integral equations. Finally, we consider the Faedo-Galerkin approximations 
of the solution and proved the convergence results. These results include the results of 
a paper by the author [12] (with Bahuguna). 

In the next two Chapters, we study a class of second order semi-linear integro- 
differential equations arising in the study of viscoelastic material with memory. 

In Chapter 4, we consider an abstract second order semi-linear integro-differential 
equation in a Banach space and first prove the existence and uniqueness of a local 
classical solution with the help of the semigroup theory and the contraction mapping 
theorem. Further, under some additional growth conditions on the nonlinear maps, 
we analyze the continuation of this solution, the maximum interval of existence and 
the global existence. This work is an extension of a paper by the author [91] (with 
Bahuguna). 


In Chapter 5, we consider a strongly damped semi-linear integro-differential equa- 
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tion and reformulate it as a second order semi-linear integro-differential equation in a 
separable Hilbert space and study the convergence of the approximate solution using 
similar techniques which are used in Chapters 2 and 3. These results include the results 
of a paper by the author [13] (with Bahuguna). 

In Chapter 6, we consider the first order quasi-linear implicit integro-differential 
equation in a Banach space. First, we prove the existence, uniqueness and continu- 
ous dependence on the initial data of a strong solution by using the method of semi- 
discretization in time in which we discretize the time axis and replace the time deriva- 
tives by the correspond difference quotients. Then we show that these discrete points lies 
in a ball whose radius is independent of the discretization parameters and prove some 
estimates. After defining the approximate solution in terms of these discrete points, 
we prove its convergence. These results include the results of a paper by author [14] 
(with Bahuguna). Next, we establish the existence of a unique local mild solution and 
its regularity under some additional condition with the help of the semigroup theory 
and the contraction mapping theorem under some different conditions. These results 
include the results of a paper by the author [15] (with Bahuguna) . 

In the closing Chapter 7, we consider some conclusions and provide some sugges- 
tions for further research work in this and related areas. 


Relevant references are appended at the end. 



Chapter 2 


Sobolev Type Evolution Equations 


One may be a mathematician of the first rank 
without being able to compute. It is possible to be 
a great computer without having the slightest idea 
of mathematics. 


— Bertrand Russell 


2.1 Introduction 

In this chapter we shall study the approximations of the solution to the nonlinear 
Sobolev type evolution equation considered in a Hilbert space. Such type of equations 
arise in the analysis of a class of functional and neutral functional differential equations 
with unbounded delay. A functional differential equation is a differential equation in 
which the unknown function occurs with various different arguments. If x be unknown 
function of time t, then the delay differential equations (or differential equation with 
retarded argument) is an equation expressing some derivative of x at time t in terms of 
X (and its lower order derivatives if any) at t and at earlier instants. 
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Consider a class of the nonlinear Sobolev type evolutions equation of the form 

^iu(t)+g{t,u{t))) + Au{t)=^f{t,u{t)), t>0, u{0) = 4>, (2.1.1) 

in a separable Hilbert space (if, ||.1|, where the linear operator A satisfies the 
assumption (HI) stated in the next section so that -A generates an analytic semigroup. 
The function / and g are appropriate continuous functions of their arguments in H. 

The existence of a unique regular solution of (2.1.1) has been proved by Hernandez 
[39] under the assumptions that -A is the infinitesimal generator of analytic semigroup 
of bounded linear operators defined on a Banach space X and / and g are appropriate 
continuous functions on [0, T] x 12 into X, where 12 is an open subset of X. Here, we are 
interested in establishing the convergence results for the Faedo-Galerkin approximations 
of solutions to (2.1.1). 

The case of (2.1.1) in which ^ = 0 has been extensively studied in literature, see, 
for instance, the books by Krein [64], Pazy [77], Goldstein [35] and the references cited 
in these books. The Faedo-Galerkin approximations of solution to the particular case 
of (2.1.1), where ^ = 0 and f{t,u) = M{u), has been considered by Miletta [72]. Our 
analysis is motivated by Miletta [72] and Bahuguna, Srivastava and Singh [16]. 

The study of (2.1.1), with linear g, has been initiated by Showalter [86, 87, 88, 89] 
with applications to some degenerate parabohc equations. Brill [25] has considered 
a pseudo-parabolic partial differential equation which appears in a variety of physical 
problems, for example, in thermodynamics [27], in the flow of fluid through fissured 
rocks [19], in the shear in second-order fluids [94] and in soil mechanics [93]. He has 
reformulated it as a semi-linear Sobolev evolution equation in a Banach space and 
established the existence results by using Schauder’s fixed point principle. He has 
used the Sobolev imbedding theorems and the established abstract results for proving 
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some existence theorenas for pseudo-parabolic partial differential equation under certain 
hypotheses. 

The nonlinear Sobolev type evolution equations of the form (2.1.1) arise in the 
study of a partial neutral functional differential equation with unbounded delay which 
can be modelled in the form (cf. [40, 41]) 

d 

— (u{t) + G{t,Ut))=Au{t) + F{t,ut), t>Q, (2.1.2) 


in a Banach space X, where A is the infinitesimal generator of an analytic semigroup 
in X, F and G are appropriate nonlinear functions from [0, T] x W into X and for any 
function u E C((—oo, oo),X), the history function Ut G C((— oo,0],X) of u, is given by 
Ut{9) = u{t+9). The above equation is called an abstract neutral functional differential 
equation (ANFDE) with unbounded delay. 


As a motivational example, for this class of equations, we consider the boundary 
value problem 


dt 

+ 


i(t,x)+ / / b{s — t,r},x)u{s,ri)dr}ds 

J ~co J 0 

/ a{s — t)u{s,x)ds ai{t,x), t>0, 

J — OO 


u{t, 0) = u{t, tt) = 0, t > 0, 

u{9, x) = 4>{9, x), 9 <0, 0 < X < TT, 


dx^ 


u{t, x) + aQ(x)u{t, x) 


1 


0 < X <11, 


where the functions oq, a, ai, bi and 0 satisfy appropriate conditions. 

The organization of this chapter is as follows. In the second section, we state some 
preliminaries and the assumptions required to establish the convergence results. In the 
third section, we consider an integral equation associated with (2.1.1). A solution of 
this integral equation is defined as a mild solution of (2.1.1). With the help of the 
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associated integral equation and the projection operators, we consider a sequence of the 
approximate integral equations and establish the existence and uniqueness of a solution 
to each of the approximate integral equations. In the fourth section, we prove the 
convergence of the solutions of the approximate integral equations and show that the 
limiting function satisfies the associated integral equation. Further, we show in this 
section that the solution can be extended to the maximal interval of existence and it is 
unique. Finally, in the fifth section, we consider the Faedo-Galerkin approximations of 
solutions and prove some convergence results for such approximations. 

2.2 Preliminaries and Assumptions 


In this section we consider some preliminaries and assumptions essential for onr purpose. 
The operator A assume to satisfy the following assumptions: 

(HI) A is a closed, positive definite, self-adjoint, linear operator from the domain 
jD(A) C II of a into H such that D{A) is dense in H, A has the pure point spectrum 

0 < -^0 ^ -^1 ^ ^2 <, • • • 

and a corresponding complete orthonormal system of eigen functions {uj}, i.e., 

Aui = XiUi and (ttj, Uj) = 
where 5ij — 1 ifi = j and zero otherwise. 

These assumptions on A guarantee that —A generates an analytic semigroup, 
denoted by t>0. 

Now, we mention some notions and preliminaries. It is well known that there exist 
constant M > 1 and a real number u such that 


|e"*^l| < Me‘^, t>0. 
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Since -A generates the analytic semigroup t > 0, we may add cl to -A for some 
constant c, if necessary, and in what follows we may assume without loss of generality 
that ||e"‘^|l is uniformly bounded by M, i.e., l|e~‘^|| < M and 0 € p{A). In this case, 
it is possible to define the fractional power for 0 < cr < 1 as closed linear operator 
with domain D{A°‘) C H (cf. Pazy [77], pp. 69-75 and p. 195). Furthermore, D{A^) 
is dense in H and the expression 

IklU = ||^“a;||, 

defines a norm on D{A°‘). Henceforth, we represent by X^, the space D{A°‘) endowed 
with the norm l|.||a- Ini view of the facts mentioned above, we have the following result 
for an analytic semigroup t > 0 (cf. Pazy [77] pp. 195-196). 

Lemma 2.2.1 Suppose that — A is the infinitesimal generator of an analytic semigroup 
e~^^, t > 0 with ||e“*^|| < M for t > 0 and 0 G p{A). Then we have the following 
properties. 

(i) Xq is a Banach space for 0 < o; < 1. 

(a) For 0 < P < a, the embedding ^ X^ is continuous. 

(Hi) A°‘ commutes with e~^^ and there exists a constant > 0 depending on a 
such that 

< ar“, t>Q. 

We assume the following assumptions on nonlinear maps / and g: 

(H2) The nonlinear map /, defined from [0,oo) x Xa into H, is continuous and there 
exists a nondecreasing function Fr from [0, oo) into [0, oo) depending on 77 > 0 such 


that 
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for all {t,u), {t,ui) and (t, U2) in [0, 00) x Bii{Xo,,(j>), where 

Br{Z, zo) = {z e Z \ \\z - zo\\z < R}, 
for any Banach space Z with its norm ||.||2. 

(H3) There exist positive constants 0<a<j3<l, L, R and 0 < 7 < 1 such that g is 
X/3- valued function and the map A^g, defined from [0, 00) x Xa into H, satisfies 

\\A^g{t,Ui) - A^g{s,U2)\\ < T{lt-sp + ||ui -usila}, 

for all (t,wi) and (5,^2) in [0, 00) x BR{Xa,<p) with L11A““^1| < 1. 

2.3 Approximate Integral Equations 

We continue to use the notions and notations of the earlier section. The existence of 
solutions to (2.1.1) is closely associated with the existence of solutions to the integral 
equation 

u{t) = + g(0, (j))) - g{t, u{t)) + [ Ae~^^~^'^^g{s, u{s))ds 

Jo 

+ f e~^^~^^'^f{s,u{s))ds, t>0. (2.3.1) 

Jo 

In this section, we will consider an approximate integral equation associated to (2.3.1) 
and establish the existence and uniqueness of the solutions to the approximate integral 
equations. By a solution u to (2.3.1) on [0,r], 0 < T < 00, we mean a function 
u e Xa{T) for some 0 < a < 1 satisfying (2.3.1) on [0,T], where Xa{T) is the Banach 
space (^([OjTjjWQ) of all continuous functions from [0,T] into Xa endowed with the 


supremum norm 
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By a solution u to (2.3.1) on [OjT), 0 < T < oo, we mean a function u such that 
u e XaiT) for some 0 < a < 1 satisfying (2.3.1) on [0, T] for every 0<T <f. 

Let 0 < To < oo be arbitrarily fixed and let 

We choose 0 < T < To such that 

||(e-«_j)^o(0 + <,(O,P»)|| < (1-rtf. 

ml— a p 

\\A’‘-%Lr' + Gi^a-^(LR+B)- + C„Fii(r„)i < (1-rt^, 

p — ot 1 — a D 

ml-a 

Ci+a-pL— h C'o.F^(To)- < 1 — iJ., 

p — a 1 — a 

where fi = R = ^/R? + \\(l)\\‘i and Ca and Ci^a-p are the constants appear 

in Lemma 2.2.1. 


Let Hn denote the finite dimensional subspace of the Hilbert space H spanned by 

{uo, ui, Un} and let P'^ : H Hn forn = 0, 1, 2, • • • , be the corresponding projection 

operator. For each n, we define 

fn:[0,T]xX^{T)^H by fn{t,u) ^ f{t,P^u{t)) 


and 

5„:[0,T]xX,(T)-.X^(T) by gn{t,u) = g{t,P^u{t)). 

We set 0(t) = 0 for t e [0, T] and forn = 0, 1, 2, • • • , define a map Sn on BR{Xa{T), 0) 
by 


(SnU){t) = e ‘^(^ + Pn(O,0)) - 5n(i,w) + [ ^^'^gn{s,u)ds 

Jo 

Jo 


(2.3.2) 
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Propostion. 2.3.1 Let (H1)-(H3) hold. Then there exists a unique function G 
BR{Xa{T),^) such that SnUn = Un for each n = 0,1,2,....; i.e., Un satisfies the ap- 
proximate integral equation 


Unit) = e + gni0,fi)) - gn{t,Un) + / Ae ^'>'^gn{s,Un)ds 

Jo 

+ / e (* ^^'^fn{s,Un)ds. 

Jo 


(2.3.3) 


Proof. We claim that Sn ■ BR{Xct{T),^) Br{Xoc{T), fi). First, we show that the map 
1 1-> {Sr,u){t) is continuous from [0,T] into with respect to norm H.Hq. For t G [0,r] 
and sufficiently small h > 0, we have 


\\(SnU){t + h) - (5nti)(t)|U 

< ||(e-“ - /)A“e-“||(||^|| + |ls(0,F”.#.)||) 

Jo 

rt-]rh 

+ / ||e-<‘+'-‘)'<yl‘+“-'>|!|1.40(|„(5.u)||<is 

+ f‘ IKe-'*'' - Ii/„(s, «)|i<is 

Jo 

/ t+/i 

||e-‘'+'‘-‘>''yl“||||/„(s,«)|l(is. (2.3.4) 

Using (H3), we obtain 

\\A^gn{t + h,u)-A^gn{Lu)\\ < L{h? P \\PMt h) - P^OWo) 

< L{h? + ||'ix(t -p h) — 'ii(t)||o,) (2.3.5) 


and 



\\e-^^-^^-^'>^A^+^-%\A^gn(s,u)\\ds 


< 


(LRP B)Ci+a-ph^~'^ 
(3 — a 


(2.3.6) 
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since 


\\A^9n{s,u)\\ < \\A^9nis,u)-A^g{s,(l))\\ + \\A^g{s,(j))\\ 

< L\\P^u{s)-cI>\\^ + B < LR + B (2.3.7) 


and 


p.3.8) 

Part (d) of Theorem 1.3.19 implies that for 0 < ■i? < 1 and x E D{A^), 

||(e-‘^-/)xl| < C!,t^\\x\U. (2.3.9) 

Let be a real number with 0 < d < min{l -a,P- a}, then A°‘y E D{A^) for any 
y E D{A°''^'^). For all t,s E [0, T], t>s and 0 < h < 1, we get the following inequalities 

||(e-'^^-J)A“e-‘^|| < C'>’'||A“+V‘^|| < (2.3.10) 


- I)A^eA^-^A\\ < 


ChP 

{t - s)“+’’ 


and 


I - /) < 


ChA 


where C = C'^ max {Ca+tf , Ci+a^^-p}- 


Using the estimates (2.3.7), (2.3.11) and (2.3.12), we get 



j)_^i+a /3g (t s)^||mh^^(5^u)||c^s 


< Ch^iLR + B) 


rpP-{a+-d) 

■{o 


/? - (a + z?) 


(2.3.11) 


(2.3.12) 


(2.3.13) 
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and 


(2,3.14) 

Jo l-{a + 'd) ^ ' 

From the inequalities (2.3.4), (2.3.5), (2.3.6), (2.3.8), (2.3.10), (2.3.13) and (2.3.14), it 
follows that (jS'jj'u)(t) is continuous from [0,2^] into with respect to the norm ||.||a. 
Now, we show SnU G Br{Xo,{T),^). Consider 

||(S'nn)(t) — (/)||q, < II (e — /)24°^(^ + ^„(0, 0)) II 

+ ll^““^IIP^5n(O,0) - A^gn{t,u)\\ 

+ r ||24^+"-V(‘-)^||||AV(3,n)||d5 
Jo 

+ /'||e-C-)-<4“||||/„(s,u)||* 

Jo 

< ii-F)§ + \\A‘-^\\L{r< + itH(t) - ,^|U} 

n^P—a rpi^a 

+ Ci+a^^{LR + B)- + C'qF^(To)- 

P — a i — O' 

— ~ f^)~^ + 0- ~ p)-Q + pR < R- 

Taking supremum over [0, T], we obtain 

\\SnU — ^\\Xa{T) < R- 

Hence, Sn maps BR(Xa{T),^) into Br{Xoc{T),^). Now, we show that Sn is a strict 
contraction on BR{Xa{T),$). For u,ve Br{Xo,{T),^), we have 

\\(SnU){t) - (5nn)(t)|U 

< \\A^~^\\\\A^gn{t,u) - A^gn{t,v)\\a 

+ [ ||A^+““^e“(*~®)^||||A^p„(s,w) - AV(5,y)IMs 
Jo 

+ f l|e“^‘“®^^2l“||||/„(s,n) - A(s,'y)||ds. 

Jo 


(2.3.15) 
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Now, 


\\Al^gn{t,u) - A^gn{t,v)\\ < L\\u{t) - v{t)\\a 

< L\\u-v\\x^^t). (2.3.16) 

Also, we have 

\\fn{s,u) - Us,v)\\ < Ff,{To)\\u{s)-v{s)\\, 

^ “ '^IU<,(T)- (2.3.17) 

Using (2.3.16) and (2.3.17) in (2.3.15), we get 
||(5„«)(t) - (5,^u)(t)|U 

/ qnp-a n^l-a \ 

< (||.4“-'>||L + + C„Fs(rc,)— ) \\u - «|U,(T). 

Taking supremum over [0,T], we get 
\\SnU — S'n^^i|j)!:o(r) 

( nn(3-a '7^1-a \ 

\\A‘-i>{\L + ) 1|« - »|U.(t). 

The above estimate and the definition of T imply that Sn is a strict contraction on 
Bjfi{Xa{T),^). Hence there exists a unique Un € Bii{Xa{T),^) such that SnUn = u„. 
Clearly satisfies (2.3.3). This completes the proof of the proposition. 

Propostion 2.3.2 Let (H1)-(H3) hold. If (j) e D{A°‘) for some 0 < a < /3 < 1, then 
Un{t) e D{A^) for all t 6 (0,r] where 0 <'& < I3 <1. Furthermore, if (p E D{A), then 
u„(t) e D(A^) for all t E [0, T] where 0<d<p<l. 
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Proof. From Proposition 2.3.1 we have the existence of a unique e Br{Xo,{T),^) 
satisfying (2.3.3). Part (a) of Theorem 1.3.19 implies that for t > 0 and 0 < < 1, 

: H — > D{A^) and for 0 < i? < /? < 1, D(A^) C D{A^). (H3) implies that the map 
t A^g{t,Unit)) is Holder continuous on [0,r]. It follows that (cf. Theorem 4.3.2 in 
[77]) 

[ e-^^-^^^A^gn{s,Ur,)dseD{A). 

JO 

Also from Theorem 1.3.9 we have € D{A) if x G D{A). The required result 

follows from these facts and the fact that D{A) C D{A^) for 0 < < 1. 

Propostion 2.3.3 Let (H1)-(H3) hold. Then for any 0 G D{A°‘), 0 < a < p <1 and 
for any to G (0, T], there exists a constant UtQ, independent of n, such that 

||«n(0llt? < Ut„ 0<t9 </?<!, to<t<T. 

Moreover, if (p E D{A), then there exists a constant Uq, independent of n, such that 

i|iin(i)lli? ^ Jo) 0 < < /? < 1, 0<t<T. 


Proof Applying A^ on both the sides of (2.3.3) and using (hi) of Lemma 2.2.1, for 
to <t <T, we have 


Jo 

+ [‘ \\e-‘-‘-‘>^AOnU3,u„)\\da 
Jo 

< + lls.(0,«) + IIA^-OlKiH + B) 

+ Ci+.a-p{LR + + Ci3Fp^{To)r^—^ - 
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Similarly 


< \\A^gn{t, Un) - A^gn{t, Wm)|| + \\A^gn{t, Um) - A^gmit, U„i)|| 

< L[||'Un(t) — 'UmWIU + T^Z^II-^'^'^m(i)i|]- 

Now, for 0 < t'o < ^ 0 ) we may write 


||'an(0 ttni(t)||(^ 

< ||e"‘^A“(^„(O,0) -5m(O,0))|| + \\A^~^\\\\A^ gn{t,Un) - A^gm{t,u^)\\ 



ll^l+a /3g (t - A^gm{s,Um)\\ds 

||^''e“(‘"®)^||||/n(s,Wn) - fm{s,Um)\\ds. 


We estimate the first term as 


||e ^^A‘^{gn{0,^) - 5m(0,^))|| 

< M||A“-^||||A^5(0,P”^) - A^5(O,P"^0)|| 

< MP“-^||L||(P"-P"^)0||a. 

The first and the third integrals are estimated as 

f ||A^+““^e"^‘“®)^||||A^5n(5,«n) - A'^gm{s,Um)\\ds 
JQ 

< 2C^+^.i3{LR + B){to - 


f° ||A“e-(‘-^)^|ll|/„(s,u„) - /„.(s,ti„.)||ds 
Jo 

< 2CaF^{To){to-t'o)-%. 
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For the second and the fourth integrals, we have 


- A»gUs.u„)\\ds 

f 

0 

< (t - (^||a„(s) - <^( 5 ) 11 , + ^ll^'’«:..(s)ll) ds 


< Ci-\-a-pL 


K:kp-oc) 


+ 


pt N 

{i - s)“f^+“~^^||u„(s) - -u,„(s)||ads 


\\A°^e ®^^||||/n(5,Wn) -/m(s,U™)||ds 

f 

0 

< C^Ff,(To) l\t - s)- (^IK(s) - «„{s)|U + 3;;^M'’u™(s)li) ds 


< C^F^iTo) 


^ yl-a pt 


- oc) 


H~ / (t s) |'U7t,(5) 1^772(5) (xds 


Therefore 


|tiTl(^) ^Tn(OI|Q: 

< ML\\A^-^\\\\{P^ - P^)<f)\U 

+ L\\A°-<^\\ (^||Uu(i) - Umit)\\a + ^ 

^ / Ci+,_^(L^ + 5) C'„F^(To)\ t/tj, 

^ ^ V ^ (to - t',r ) ° 


‘ ^ CgFfSTo) Ci+g-^T 

'f, \ [t-sY {t-sY+'^-P 



jpl—a rj-'P—a 

Ca,p = ^aP'Ri'^o)YZYt 


where 
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Since ||A“ ^\\L < 1, we have 




< 


(1 - ||A-^||L) 




U, 


tk 


X-d-a 

Uf 

On 


I p I ^i-\-cx-p{.LR + B) C'o;F^(ro) \ . p “tp 

[ (to (to-t;)”j “''’ajt' 


+ 



^qF^(Fo) ^ Ci^a- pL 


{t — s)°‘ {t — 




Consider 


f C^F^jTo) \ 

y (t — s)“ (t — 5)^+“”^y 


< Ni 


+ 


1 


(t — s)“ (£ — 


where 


Ni — max{C'aF^(ro)) C'i+a-/3F}. 
Using this estimate in (2.4.1), we get 


Un{t) - Um(i)l|a 


< 


(1-|1A-/^1|L) 


M\\{P--P^)cP\\^ + \\A°‘-^\\L 


C'l+g- ^(LF + F) . CaF^(ro) 


\ 1?— a 


, r, I ^ I +/ , /-r ^ ‘'O 

+ 2 A|F 


(2.4.1) 
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Hence from lemma 1.3.27, there exists a constant K such that 


I ^m(^)l|Q; 


< 


(1 - ||A-/>||L) 


M||(P" - P^)4,\U + (M“-0||L + 


+ 2 


f c^+a-^(LE + B) C^F^(To) ] \ 


Taking supremum over [to, T] and letting m — > oo, we obtain 


lim sup \\Un{t) - Um{t)\\a 

{n>m, fe[£o,ri} 

< 2 ( c^^^.^{LR + B) C^Ff,{To) \ 

- (1 - ||A-/^||L) {to-t'orj 

As to is arbitrary, the right hand side may be made as small as desired by taking t'o 
sufficiently small. This completes the proof of the proposition. 


Corollary 2.4.2 If 4> ^ 'then 

lim sup \\Un{t) - Um(t)lla = 0. 

{n>m, 0<t<T} 

Proof. Propositions 2.3.2 and 2.3.3 imply that in the proof of Proposition 2.4.1 we may 
take to = 0. 

For the convergence of the solution Un{t) of approximate integral equation (2.3.3) 
we have the following result. 

Theorem 2.4.3 Let (H1)-(H3) hold and let (p 6 D{A°‘). Then there exists a unique 
function u G X^iT) such that Un u as n oo in XofT) and u satisfies (2.3.1) on 
[0,r]. Furthermore, u can be extended to the maximal interval of existence [0, tT„o®), 
0 < tmax < oo satisfying (2.3.1) on [0,tmax) eind u is a unique solution to (2.3.1) on 
[O , tmax ) • 
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Proof. We first assume that (p € D{A). Corollary 2.4.2 implies that there exists 
u e XaiT) such that Un converges to u in Xa- Since G BR{Xa{T),^) for each 
n, u is also in BR{Xa{T), 0) . Further, we have 

\\gn{0j)-9{0,m < M0,P^<P)~giO,m < L||(P”-J)0|U (2.4.2) 

and the right hand side of (2.4.2) tends to zero as n oo. Also, 

\\fn{t,Un) - f{t,u{t))\\ 

< \\fn{t,Un) - fn{t,u)\\ + j|/„(t, u) - /(t, w(t)) || 

< ^«(ro)(lkn(t) - ^(t)iu + ll(i"" - iHt)U 

Taking supremum over [0,r], we get 

sup \\fn{t,Un) - f{t,u{t))\\ 

0<t<T 

< fs(ro)(IK-«lk(T) + ll(P”-/)«lk(T)). (2.4.3) 

The right hand side of (2.4.3) tends to zero as n — > oo. Prom (H3), we get 

\\A^gn{t,Un) - A^y(t,u(t))|| 

< \\A^gn{t, Un) - A^Qnit, u)|| + \\A>^gnit; u) - A^ g{t, u(t))|| 

< L(iK(t)-«(t)iu + ii(P"-j)u(t)y. 

Taking supremum over [0, T], we get 

sup \\A^ gn{t,Un) - A^ g{t,u{t))\\ 

0<t<T 

< P(||un-ulU.(T) + ll(^”--fHUa(T))-^0 as n->oo. (2.4.4) 
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Using (2.4.2), (2.4.3), (2.4.4) and the bounded convergence theorem in (2.3.3), we get 

rt 

u{t) = e~^^{^ + g(0,(l)))-g{t,u{t))-\- Ae~^^~^^^g{s,u{s))ds 

Jo 

+ [ e-^^-^^^f(s,u(s))ds. ( 2 . 4 . 5 ) 

Jo 

Thus, for (j) e D{A), there exists a unique function u € XaiT) such that u as 

n oo and satisfies integral equation (2.3.1) on [0,T]. 

Now, let us assume that 0 G D{A°‘). Since, for 0 < t < T, A°‘Un{t) converges 
to A°‘u{t) as n — > oo and ^^^(O) = tt(0) = 0 for all n, we have, for 0 < t < T, 
A°‘Un{t) converges to A°‘u(t) in H as n oo. Since Un G BR{Xa{T), it follows that 
u G BR{Xct{T), 0) and for any 0 < to <T, 

lim sup ||un(t) — 'u(t)||Q, = 0. 

{to<t<T} 

Also, 


sup \\Jn{t,Un) - f{t,u{t))\\ 

to<t<T 

< ■^H(^'o)(||un-u|U,(T) + ||(i"”-/)u|U.{r))^0 as n ^ oo 


and 


sup \\A^gn(t,Un)-A^g{t,u{t))\\ 
to<t<T 

< T(||w„ - ■u||x„(T) + ||(P” - d)w||xe,(T)) 0 as n-^oo. 
Now, for 0 < to < t, V7e may rewrite (2.3.3) as 


Un{t) = e-*^{4> + gn{0,^))-gn{t,Un)+ Ae ^^^gn{s,Un)ds 

^ /*^o rt\ 

+ ] e~^^~^^^fn{s,Un)ds. 

^Jo Jto / 
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The first and third integrals are estimated as 


"to 


Ae M)ds 



rto 

< 

1 


Jo 

< Ci-p{LR + B)T^-Ho, 


^to 


e ^^'^fn{s,Un)ds 


< MF^(ro)to. 


Thus, we have 


u. 


n{t)-e *^(^ + 5n(O,0)) +5„(t,u„) - / Ae ^'^^gn{s,Un)ds 


'to 


e "'^^fn{s,Un)ds 


'to 


< {Ci.g{LR + B)T^-^ + MF^{To))to. 


Letting n — >■ oo in the above inequality, we get 


u 


(t) - e + g{0, 4>)) + g{t. u{t)) - Ae u{s))ds 


'to 


e ®^'^/(s,u(s))ds 


'to 


< {Ci-0{LR + B)T^-^ + MF^(To))to. 


Since 0 < to < T is arbitrary, we obtain that u satisfies the integral equation (2.3.1). 

If u satisfies (2.3.1) on [0,r], then we show that, u can be extended further. Since 
0 < To < oo, was arbitrary, we assume that 0 < T < Tq. We consider the equation 


J 

—{w{t) + G(t, w{t))) + Aw{t) = F{t,w{t)), 0 < t < To < oo, 

(JLL 


w(0) = u{T), 
where F ; [0, To — T] x ^ iJ is defined by 


F{t,x) = f{t + T,x), 
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and G : [0, To - T] x ^ Xp is defined by 

= 9{t + T, x), 

for (f, x) E [0, To — T] X Xa- We note that F and G satisfy (H2) and (H3), respectively, 
where To is replaced by To — T. Hence, there exists a unique w € C([0, Ti],Xa), for 
some 0 < Ti < To — T, satisfying the integral equation 

w{t) = e-^^iu{T) + G{0,u{T)))-Git,w{t)) 

+ [ Ae~^*~^^'^G{s,w(s))ds+ f e~^^~^^"^F{s,w{s))ds, 0<t<Ti. 

Jq Jo 


Now, we define 


u{t) = 


u{t), 0 < t < T, 

w{t-T), T<t<Ti+T. 


Then, u satisfies the integral equation 


u{t) 


e + g{0,(f})) - g{t,u{t)) + [ Ae ^^^g{s,u{s))ds 

Jo 

+ f f{s,u{s))ds, 0 < t < Ti + T. 

Jo 


(2.4.6) 


To see this, we need to verify (2.4.6) only on [T, Ti + T]. For t € [T,Ti + T], 


u{t) = w{t — T) 

= e-(‘-'")^(«(T) + G(0, u(T))) - G{t - T, w{t - T)) 

+ [ Ae~^^~^~^'>^G{s,w{s))ds + [ e~^^~'^~^'>^F{s,w{s))ds. 

Jo Jo 
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Putting T + s = 77, we get 


u{t) = e ^^{4> + g{0,4>))-g(T,uiT)) 

+ f Ae ^^'^g{s,u{s))ds + ( e“^^“'®)^/(s,u(s))cZs} 

Jo Jo 

+ GiO,u(T)))-G(t-T,w(t-T)) 

+ Ae-^^-'^^^Gi'n - T,w(r] - T))dg 

+ j e-^^-^^^F{r)-T,w{r]-T))ds 

fT 

= e-*^{4> + p(0, 0)) - g{t, w{t -T))+ / u{s))ds 

JO 

+ J Ae-^^-^^^g(s,w{s-T))ds + e-^^-^^^f{s,u{s))ds 
+ J e~^*~^'>^f{s,w{s-T))ds, 

as G{0,u(T)) = g{TyU{T)), G{t—T,w{t — T)) = g{t,w{t—T)) and F{t — T,w{t — T)) = 
f{t,w(t — T)). Hence, we have 

u{t) = e"*^(</> + g{0, 0)) - g{t, u(t)) + / u{s))ds 

Jo 

+ [ e“^‘“®)^/(s,tt{s))ds, 

Jo 

for t G [0,Ti+r]. Thus, we see u{t) satisfy (2.3.1) on [0,ri+T]. Hence, we may extend 
u{t) to maximal interval [0,imax) satisfying (2.3.1) on [0, tmoa:) with 0 < tmax < oo. 

Now, we show the uniqueness of solutions to (2.3.1). Let ui and U 2 be two solutions 
to (2.3.1) on some interval [0, T], where T be any number such that 0 < T < tmax- Then, 
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for 0 < i < T, we have 


||Kl(t) -M2 W||c 

Jo 

+ [ ||e-M'‘A«l|l/(s,u,(s)) - fisMmds 

Jo 

+ C,+._^L f (t-<i)-(‘+«-«||ui(s)-U2{s)|Uds 

Jo 

+ CaF^{tmax) f {t - s)““||«i(5) - U 2 (s)|U'^S. 

Jo 

Since, ||A“''‘^||L < 1, we have 


\\Ui(t) - W2Wlla 
1 


< 


< 


C'l+a— ffL C(xF ^{t 

max') 


+ 


ii-\\A-qL)J, \{t-sY+--p ■ 
N2{t^-^ + 1) 


||ui(s) - W2(s)llarf5 


(1 - IIA-^IIL) 


r ^ 

Jo i^-s) 




where 

N2 = m.aX-[CQ.i^^(ima2)) ^l+a—pFY 
Using Lemma 1.3.27, we get 

||ui(t) -U2(t)l|Q = 0 

for all 0 < t < T. Prom the fact that 

i|ai(t)-«2(t)|i < ^\\ui{t) - U2{t)\\a, 

it follows that ui = U 2 on [0,T]. Since 0 < T < tmax was arbitrary, we have ui = U 2 on 
[Oi^moi)- This completes the proof of the theorem. 
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2.5 Faedo-Galerkin Approximations 


For any 0 <T < tmax, we have a unique u G Xa{T) satisfying the integral equation 

u{t) = e-^^{4) + g{0,(l)))- g{t,u{t))+ f Ae~^^~^'>^g{s,u{s))ds 

Jo 

+ [ e~^^~^'>'^f{s,u{s))ds. ( 2 . 5 . 1 ) 

Jo 

Also, we have a unique solution Un G Xc,{T) of the approximate integral equation 

Unit) = e~‘^(0 + ^n(O,0)) -5n(t«n) + / 

Jo 

+ f e~^^~^^^fn{s,Un)ds. (2.5.2) 

Jo 

If we project (2.5.2) onto Hn, we get the Faedo-Galerkin approximation u„(t) = P"'Un{t) 
satisfying 


Unit) = e-^^iP-(f> + P^giO,P^<j>))-P^git,Unit)) 

+ / Ae-(‘-*)^P"p(s,Un(s))ds+ f e-^^-^^^P^fis,Unis))ds (2.5.3) 
Jo Jo 

The solution u of (2.5.1) and of (2.5.3), have the representation 

OO 

uit) = Y^aiit)ui, aiit) = {uit),Ui), z = 0,r,...; (2.5.4) 

i=0 


“n(i) = y~^Q:”(t)ut, a^it) = iunit),Ui), i = 0, 1, ...n. 
i=0 


(2.5.5) 
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Using (2.5.5) in (2.5.3), we obtain the following system of first order ordinary differential 
equations 


_d 

dt 




<(«))) + A<<(() = f”{t,a;{t), 

a?(0) = (2.5.6) 


where 


,<(()) 


and (f)i = {(j), Ui) for z = 1, 2, n. 




i=0 



) 




i=0 


\ 

Ui 

J 


The system (2.5.6) determines the a”(t)’s. Now, we shall show the convergence of 
Oiit) — >■ ai{t). It can easily be checked that 


A°‘[u{t) — u{t)] 




^(Q:i(t) - a2{t))ui 


1=0 


5^ A“(ai(t) - a^{t))ui. 


i—0 


Thus, we have 

||A“[«(t)-a(i)jlP > 

i=0 

We have the following convergence Theorem. 


Theorem 2.5.1 Let H1)-(H3) hold. Then we have the following . 


(a) If (f) G D{A°‘), then for any 0 <to <T, 


Af(t>,(t) -<{())=' 

.i=0 


= 0 . 


lim sup 
to<t<T 
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(b) If (j) E D{A), then 


n 


lim sup 

n-*oo o<t<T 


-<(«))“ 


= 0 . 


The assertion of Theorem 2.5.1 follows from the facts mentioned above and the following 
proposition. 


Propostion 2.5.2 Let (Hl)-(HS) hold and let T be any number such that 0 < T < 
traax, then we have the following. 

(a) If 4> E D{A°‘), then for any 0 <tQ <T, 

lim sup \\A°‘[un{t) - •u„i(t)]|| = 0. 

{n>m, to<t<T} 

(b) If f) E D{A), then 

lim sup ||A“[un(t) - 'hm(t)]|| = 0. 

{n>m, 0<t<r} 


Proof For n>m,V7e have 

- F”u„(t)]|| 

1|P“K(() - U„(t)l||„ + ||(P” - 
ll^n(i) ~ A ^m||- 


< 

< 


114> E D{A°‘), then the result in (a) follows directly from Proposition 2.4.1. If 0 G D{-^) 
(b) follows from Corollary 2.4.2. 



Chapter 3 


Second Order semi-linear Evolution 
Equations 


An approximate answer to the right problem is 
worth a good deal more than an exact answer to 
an approximate problem. 


— John Tukey 


3.1 Introduction 


In this chapter, we shall study the approximation of the solution to a second order 
semi-linear evolution equation. We first show that a strongly damped semi-linear wave 
equation may be written as a second order semi-linear evolution equation. 


Let n be a bounded domain in with sufficiently smooth boundary dQ, and let 



be a symmetric second order strongly elliptic differential operator in Ct. We are con- 
cerned with the following initial-boundary value problem for a strongly damped semi- 
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linear wave equation, 

f dxi\ 

{x,t) + (aL + bl) j (a:, t) + {cL + dl)u{x, t) 

3il 

= h{x,t,u{x,t),~{x,t)), € Q X (0,r), (3.1.1) 

dll 

u{x,o) = (f){x), —{x,o)='tp{x), xen, 

with the homogeneous Dirichlet boundary conditions, where a > 0, 6, c and d are 
constants and h is a smooth function of its arguments. 

Under suitable conditions on the data h, 0 and ■0, we reformulate (3.1.1) as the 
following initial value problem for the abstract second order semi-linear equation, 

^(t) + A (t) -f Bu{t) = f(t, u{t), ^(t)), t > 0, 

w(0)=no, ^(0)=xi, (3.1.2) 

in the separable Hilbert space H = L^{Q,) where the linear operator A with the domain 
D{A) = n Hq{Q,) is given by 

Au = aLu, u € D{A), 

and the operator B is such that D{B) = D{A) with Bu = Ci.A -t- C 2 / for some constant 
Cl and C 2 . The function / is defined from [0, 00 ) x H x if into if given by 

f{t,u,v) = h{t,u,v)—bv. 

D. Henry [38] has considered the coupled system 

dv d*}! 

— + Av = f{t,u,v), — = g{t,u,v), (3.1.3) 

and proved the existence of an invariant manifold (collection of solutions) and estab- 
lished the stability theory for this system under certain assumptions on A, f and g. We 
may write (3.1.2) as a system of the first order equations like (3.1.3), where g{t, u, v) = v. 
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The problem (3.1.1), for the particular case L = -A, the n-dimensional Laplacian, 
is dealt with by Duvaut and Lions [31] and Glowinski, Lions and Tremolieres [34] in 
the context of some viscoelastic materials using the method of variational inequalities. 

Sandefur [84] has considered the problem (3.1.2), for a special case where the 
operators A and B are such that the equation can be written in the following factorized 
form. 


+ A 2 ) + A^A^u = f{t,u), t e (o,r], 

/ V du V . 

u(0)=a:o, — (0) = 2 :i. (3.1.4) 

Using the theory of semigroups together with the method of successive approximations, 
Sendefur [84] applied the results to various physical problems, for instance, damped 
and strongly damped semi-linear equations, the telegraph equation and the equation of 
motion for a thin panel. 

Aviles and Sandefur [3] have modified the techniques of Sandefur [84] by changing 
the order of integration in the definition of the mild solution to (3.1.4) and applied 
the results to the Klein-Gordon equation, the von Karman equation and the vibrating 
beam equation. 

Webb [95] and Ang and Dinh [2] have considered the following initial boundary 
value problem, 

Utt — XAut + f{u) = 0, {x, t) e 9, X (0, T), A > 0, 

u{x,t) = 0, .(x, t) e 50 X [0,T), 

u(a;,0) = mo(a:), Ut(x, 0) = Wi(a:), (3.1.5) 


where O is a bounded domain in R'^ with sufficiently smooth boundary 50. For 1 < 
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N <2, (3.1.5) governs the motion of a linear Kelvin solid (a bar if = 1, a plate if 
N = 2) subject to nonlinear elastic constraints. Webb [95] has studied the existence 
and asymptotic behavior of solutions to (3.1.5) for N = 1,2,3. Ang and Dinh [2] 
have generalized the results of Webb [95] for (3.1.5) using the method of successive 
linearizations. 

Engler, Neubrander and Sandefur [32] have used the theory of semigroups to prove 
the existence and uniqueness of mild solutions of (3.1.2). Balachandran and Park [17] 
have also proved the existence of a mild solution by using the theory of strongly con- 
tinuous cosine families and the Schaefer fixed point theorem. Bahuguna [6, 7] has 
established the existence, uniqueness, continuation of a solution to the maximal inter- 
val of existence, and the global existence of a strong solution and a classical solution 
under different conditions. 


In this chapter we have shown the convergence of the Faedo-Galerkin approxima- 
tion to the solution of the abstract second order semi-linear evolution equation (3.1.2). 
Our aim is to extend for second order semi-linear equation the ideas and techniques of 
Bahuguna, Srivastava and Singh [16] which are applied to obtain the Faedo-Galerkin 
approximation of the solution to a first order semi-linear integro-differential equation 
of the form 


'{t) -f Au{t) = f{t, u{t)) + f a{t — s)g{s, u{s))ds, 0 <t <T < oo, 

Jo 

u{0) = (j). (3.1.6) 


In [16] , authors have used the idea of Miletta [72] to establish the convergence of the 
Faedo-Galerkin approximations of the solution to (3.1.6). Miletta [72] has established 
results for the first order semi-linear evolution equation 


u\t) + Au{t) — M{u{t)), u{0) = (f), 
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in a separable Hilbert space and A satisfies similar conditions as considered in this 
Chapter and M is a. nonlinear map defined on D{A°‘), for some o;, 0 < o: < 1, which 
satisfies a Lipschitz condition on ball in D(A°‘). 

In [97], Zarubin find estimates for the rate of convergence of the Faedo-Galerkin 
method for the Cauchy problem u' + Au + K(t)u = h{t), u{0) = 0, while in [98], for 
the problem u'{t) + A{t)u{t) + K{t,u(t)) = h{t), n(0) = 0, when A{t) is a strongly 
continuous and differentiable operator with a domain that does not depend on t, and 
K{t, .) is a nonlinear Frechet differentiable operator, that is subordinate to A(t) with 
order less than one. 

Despite the widespread use of the Faedo-Galerkin method (in many applications 
it is referred to as the method of harmonic balance) , the convergence behavior in many 
cases is not known. Bazely [20, 21] showed the uniform convergence of the approxima- 
tions to solutions of the non-linear wave equation 

+ Au{t) + M{u{t)) = 0, u(0) = (j), u'{0) — Ip, 

on any closed subinterval [0, T] of existence of the solution. 

In [30], Doronin, Lar’kin and Souza have proved the existence of a global general- 
ized solution to initial boundary value problem for the wave equation with nonlinear sec- 
ond order dissipative boundary conditions with the help of the Faedo-Galerkin Method, 
apriori estimates and compactness arguments. Bessaih and Flandoli [23] investigated 
a two-dimensional Euler equation subject to a stochastic perturbation (a noise). An 
existence and uniqueness result is proved under some assumptions of spatial regularity 
on the noise by using Faedo-Galerkin method. Bellery and Fata [22] have considered 
the strongly damped semi-linear equation in R^, 


Utt — Aut — Au-h g{x,u) + (p{x)ut = f{x,t) in R^x(0, T], 
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with the initial conditions, u(x,0) = uq{x), Ut{x,0) = vo{x), and the condition that 
lim|i|—oo = 0 3,11 t G [0,T]. In this, proof of the existence and uniqueness 

results is carried out via standard Faedo-Galerkin approximation scheme. 

The organization of this chapter is as follows. In the second section we state 
the preliminaries and the assumptions required to establish the convergence results. 
We consider an equivalent form of (3.1.2) in a Hilbert space also in this section. In 
the third section we consider a pair of associated nonlinear integral equations to this 
equivalent form. The solutions to these integral equations are defined as mild solutions 
to (3.1.2) or to the equivalent problem. Using the pair of associated nonlinear integral 
equations and the projection operators, we consider a pair of the approximate nonlinear 
integral equations and we show the existence and uniqueness of solutions to this pair of 
approximate integral equations using the Banach contraction principle in this section. 
In the fourth section we establish the convergence of solutions and convergence of the 
pair of approximate integral equations to the pair of associated integral equations, limit 
of the solutions to the approximate integral equations being the solution of the pair of 
associated integral equations. Further, we show, in this section, that the solution can 
be extended to the maximal interval of existence and it is unique. Finally, in the fifth 
section we consider the Faedo-Galerkin approximations of solutions and prove some 
results concerning the convergence of such approximations. 

3.2 Preliminaries and Assumptions 

In this section we consider some preliminaries and assumptions essential for our purpose. 
The operator A and the nonlinear map / assume to satisfy the following assumptions. 


(Hi) A is a closed, positive definite, self-adjoint, linear operator fi:om the domain of 
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A, D(A) C H into H such that D{A) is dense in H. A has a pure point spectrum 

0 <! Aq ^ Aj ^ A 2 • 

and the corresponding complete orthonormal system of eigen functions {uj}, i.e., 

Aui = XiUi and {ui,Uj) = Sij, 
where 6ij = 1 i = j and zero otherwise. 

If (HI) is satisfied, then —A generates an analytic semigroup in H which we denote by 
t > 0. It follows that the fractional power A“ of A for 0 < a < 1 are well defined 
from D(A°‘) C H into H (cf. Pazy [77], pp. 69-75). D{A°‘) for 0 < o; < 1 is Banach 
space endowed with the norm 

llxlU = llA"x||, xeD{A^). 

Henceforth, we represent by Xa, the space D{A°‘) endowed with the norm |1 .||q; for 
0 < a < 1. 

(H2) The map / is defined from [0, 00 ) xXi x into H and there exists a nondecreasing 
function Fr from [0, 00 ) into [0, 00 ) depending on i? > 0 such that 

\\f{t,u,v)\\ < FR{t), 

\\fit,Ui,Vi) - f{t,U2,V2)\\ < Fii(t){||ui-ii2||i + il'yi-^^2|U}, 

for all (t,u,v) and {t,ui,Vi), {t,U 2 ,V 2 ) in [0, 00 ) x Br{Xi x Xq, (xq, 2 : 1 )), where 

X Xo,, (2:o,2:i)) = {(a;,y) € Xi X Xq, : jja; - xo||i + jjy - a:i||a < R}. 
Without loss of generality, we concentrate on the following problem, 

ti(0) = Xo, «'(0) = XI, 


(3.2.1) 
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in a separable Hilbert space H, since the other terms in (3.1.2) may be transferred on 
the right and merged with / and the changed / still satisfies (H2). 


3.3 Approximate integral equations 


The existence of solution to (3.2.1) is closely associated with the following pair of integral 
equations (cf. Bahuguna [6, 7]), 


u{t) = XQ-\-{e + f (e - I){-A) ^ f{s,u{s),v{s))ds, t>0, 

Jo 

(3.3.1) 

v{t) = e~^'^xi+ f e“^‘“'®^'^/(s,u(s),ti(s))ds, t > 0. 

Jo 


In this section we will consider a pair of approximate integral equations corresponding 
to (3.3.1) and establish the existence and uniqueness of the solutions to the pair of 
approximate integral equations. By a solution (u, v) to (3.3.1) on [0,T], 0 < T < oo, we 
mean a function {u,v) G Xi(T) x Xa{T) for some 0 < a < 1 satisfying (3.3.1), where 
Xi{T) X Xa{T) is the Banach space C([0, T], X\ x X^) of all continuous functions from 
[0, T] into Xi X Xa endowed with the supremum norm 

ll(u,'y)IUi(r)xX„(T) = liw||xi(T) + 


where 


||w|Ui(r) = sup llAM(t)|| = sup l|«(t)||i 

0<t<T 0<t<T 


and 


iklU„(T) = sup P“u(t)|| = sup \\v{t)\\a. 

0<t<T 0<t<T 

By a solution {u,v) to (3.3.1) on [O.T), 0 < T < oo, we mean a function {u,v) E 
Xi{T) X Xa{T) for some 0 < a < 1 satisfying (3.3.1) in [0,T] for every 0 < T < T. 
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Since -A generates the analytic semigroup t > 0 we may add cl to -A for a 
suitable constant c, if necessary, and assume without loss of generality that ||6~*'^|| M 
and -A is invertible. Furthermore, it follows that commutes with and there 
exists a constant Cq > 0 depending on a such that 

Let 0 < To < oo be arbitrary fixed and 

HR) = (1 + R)F^{T„), 

where 

R=- R 1 + R 2 , Ri = 1^/22 + 11 x 0 III and R 2 = yjR? + ||a;i||2. 

Let 0 < T < To be such that 

o|’g>^{||(e"‘^-^)2:i|| + ||(e-‘^-/M'"^^ < I 

and 

T < min|To,|(M + l)-'L(it:)-\[|c-Hl-a)i:(i?n^ 

Let H-n denote the finite dimensional subspace of H spanned by {%, Ui, . . . , u„} and for 

n= 1,2, . . ., let P” : H Hn be the corresponding projection operator. For each n, 
we define 

/„:[0,T]xZi(r)xX,(T)^P, 

such that 

fn{t,u,v) = f{t,P^u{t),PMt))- 

We set xo{t) = xo and xi{t) = xi for t G [0, T]. Let Wr = Pr(Xi(T) x Xa{T), {xo, fi)), 

where 


MMT) X Xo.{T),{x~o,xR)) 

= {{yi,y2) E Xi{T) X Xa{T) : ||yi — iolUi(r) + II2/2 - ^ilUa(r) < R}- 
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Define a map Sn on Wr such that Sn{u,v) := {u,v) with 

Jo 


v{t) = e + 





(3.3.2) 


Propostion 3.3,1 Let (HI) and (H2) hold. Then there exists unique {un,Vn) E Wr 
such that Sniun,Vn) = {un,Vn) for each n = 1,2,..., i.e., (u„,Un) satisfies the pair of 
approximate integral equations 


Un{t)=Xo-{e *^-I)A ^xi- [ (e - I)A 7n(s,u„,u„)(is, 

Jo 

(3.3.3) 

Vn{t) = e“‘^Xi + f e~^^~^'>^fnis,Un,Vn)ds. 

Jo 


Proof. First, we claim that Sn : Wr — >• Wr. For this we need to show first that the 
map t t-> {Sn{u,v)){t) is continuous from [0,T] into Xi x Xa with respect to norm 
||.||i + ||.||a. For ti,t 2 E [O,!"] with ti < t^, we have 


\\uit 2 ) - fi(tl)lli + Ilfi(t2) - 'y(^l)||a 

< - e-‘^^)xi|U 

+ / - m\\fnis,u,v)\\ds 

Jti 

+ f -e"(*i“®)^||||/n(s,u,u)||ds 

Jo 

+ r||e-(‘-^)^yl“||l|/„(s,u,u)||ds 
Jti 

+ / ll(e“^‘^"^^ -e~^*^""))A“|lll/„(s,u,u)llds. 
Jo 


(3.3.4) 
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Also, we have 

r - I\\\\U{s,u,v)\\ds < {M+l)F^{To){h-h), (3.3.5) 

Jti 

f" \\e-(^^-^)^Am\Us,u,v)\\ds < Cai^K(ro)— ■ (3.3.6) 

Jh 1 ^ 

Part (d) of Theorem 1.3.19 implies that for 0 < < 1 and x € D{A'^), we have 

||(e-^^-/)a:|| < C,t%x\U. (3.3.7) 

If 0 < < 1 is such that 0 < o: + 1 ? < 1, then A°‘y € D{A^) for any y G D(A“+’’). 

Therefore, for t, s G [0,r], we have 

|l(e-*^ - /)A“e-*"^a:|| < C'^t^\\A^e-^^x\\^ = C!gf\\A^^e-^^x\\ 

< C';a+,5t’’s-(“+’’)||:rl|. (3.3.8) 


From (3.3.8), we get 




< C;C'„+,(t2-ti)’'(ti-5)-^“+’’). 


Hence 



- e ^^^)A“llll/„(s,w,u)||ds 

rpl'—(cx+'&) 

< - *>)'• 


(3.3.9) 
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Also, from (3.3.8), we have 


Therefore 


< C'Mh-kYih-s)-^. 


Hence 


Jo 

< c;c,F^(r„)^(<2 - h)". (3.3.10) 

Prom inequalities (3.3.4), (3.3.5), (3.3.6), (3.3.9) and (3.3.10), it follows that {Sn{u, v)){t) 
is continuous from [0, T] into Xi x with respect to norm ||.||i + ||.|1 q,. Now, we show 
that Sn{u,v) G Wr, i.e., {u,v) G Wr. Consider 

llu(t)-Xolll + ||u(t) - Xilla 

< 11(6-“ - + ||(e-“ - ;)A“ii|| 

+ [ lle-C-*'' - /||||/„(6.«,t-)|l<is + ||e-<‘-“.4“||||/„(s,«,»)||d6 

< ? + (M + i)Fs(ro)r + c„Fj(ro) / (t-s)'“<fa 

o Jo 

< j + {M + l)L{R)T + CaL{R)^^<R. 


Taking supremum over [0, T], we obtain that Sn maps Wr into Wr- Now, it remains 
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to show that Sn is contraction on Wh- For (ui, vi), {u 2 , V 2 ) € Wr, we have 
lluiCt) - M2(i)|li + - 'y2(i)lla 

< f - I\\\\fn{s,Ui,Vi) - fn{s,U2,V2)\\ds 

Jo 

+ f - fn{s,U2,V2)\\ds. 

Jo 

From (H2), we have 

\\fn{s,Ui,Vi) - fn{s,U2,V2)\\ 

= \\f{s,P^U,is),P%^{s)) - f{s,P^U2{s),P'^V2{s))\\ 

< FR{n){\\P^m{s) - P^n 2 {s)h + \\P^v,{s) - P^V 2 {s)\U) 

^ ■F'^(^o)(||wi(s) - ^^2(3)111 + ||ui(s) - W2(s)IU) 

^ ■F!r(^o)(1|Mi - W2|lxi(T) + 11^1 - ■y2|Ua(T))' 

Hence 

||ui(t) -'U2(01ll + llt'i(i) -'^sWIU 

< |(M + 1)Fr{To)T + Ca^FRiTo) 

x(||ui - U 2 \\xi{T) + Ibl - 'V 2 \\Xa{T)) 

x(||ui — U 2 \\xi{T) + 11^1 ~ V 2 '^Xcc{T)) 

2 

^ gdl^l — '*^2llxi{T) + 11^1 “ ‘y2|iXc<(T))- 
Taking supremum over [0,T], we get 

2 

11*1 - *2||xi(T) + 11*1 - *2l|xc(r) < “ “2||xi(r) + II^^I - ■y2l|xa(T))- 
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Thus, Sn is a strict contraction on Wr. Hence, there exists a unique («„, Vn) 6 Wr such 
that Sn{un,i^n) — Clearly, {un,Vn) satisfies (3.3.3). This completes the proof 

of the proposition. 

Propostion 3.3.2 Let (HI) and (H2) hold. If (xo,Xi) e D{A) x D(A“) for some 
0 < a < 1, then {un{t),Vn{t)) e D(A) x D{A^) for all t € (0,r] where 0 < < 1, 

Furthermore, if (a:o,a;i) € D{A) x D{A), then (u„(t),u„(t)) e D{A) x D{A^) for all 
t G [0, T] where 0 < i? < 1. 

Proof. From the Proposition 3.3.1 we have the existence of a unique (u„, Vn) G Wr 
satisfying (3.3.3). Clearly, Un{t) G D{A). Part (a) of Theorem 1.3.19 implies that 
g-tA . ^ D{A^) for t > 0 and 0 < i? < 1. Also, from Theorem 1.3.6, we have 

e~^'^x G D{A) if X G D{A) and for x £ H, jQe~^"^xds G L)(A). The result of the 
proposition follows from these fact and the fact that D(A) C D{A^) for 0 < ?? < 1. 

Propostion 3.3.3 Let (HI) and (H2) hold. Then for any Xi G D{A°‘), 0 < a < 1 and 
any to G (0,T], there exists a constant independent ofn, such that 

lbn(i)l|i? < Vt^, 0 <'&<!, to<t<T. 

Moreover, if Xi G D(A), then there exists a constant Vo independent ofn, such that 

lbn(t)ll,5 < VJ), 0<i9<l, 0<t<T. 

Proof. Applying A^ on both the sides of second integral equation of (3.3.3) and 
using part (c) of Theorem (1.3.19), we have for to<t <T, 

IKWIIs < 

Jo 

< c„r'’||xi|i + c,fs(r„)^ 

< O„t^0\\x,\\+CtL(R)PH < Vt,. 
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Since xi £ D{A) implies that Xi e D{A'^) for 0 < t? < 1, we get 

||«»(t)IU < + C.,L{R)^ 

< M||/l'’xi|| + C<L(ii)^ < Vo. 

1 — V 

This completes the proof of the proposition. 

3.4 Convergence of approximate solutions 


In this section we establish the convergence of the solution (un,u„) € Xi{T) x Xa(T) 
of approximate integral equations 




(e-*^ - I)A-^x, - - l)A-^fn{s,Un, 

Jo 


Vn)ds, 


(3.4.1) 


i(t)=e ‘^Xi+ / e "^^fn{s,Un,Vn)ds, 
Jo 


to a unique solution {u,v) of (3.3.1). For proving the convergence, we need following 
stronger assumption than (H2) on the nonlinear map /. 


(H2') The map / is defined from [0, oo) x Xi x Xa into for 0 < a < ^ < 1 and 
there exists a nondecreasing function Fr from [0, oo) into [0, oo) depending on i? > 0 
such that 

\\f{t,u,v)\\i3 < Fuit), 

\\f{t,Ui,Vi)-f{t,U2,V2)\\p < FR(t){||ui -U2lii+ llui -U2|la}, 

for all {t,u,v) and (t,ui,Vi), {t,U 2 ,V 2 ) in [0, oo) x Br{Xi x Xq., (xo,a:i))- 

We can easily check that the changed / in (3.2.1) still satisfies (H2) with FR{t) = 
Ci-^R(t) + C 2 for some constant Ci and C 2 but this changed / does not satisfy (H2 ). 
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Propostion 3.4.1 Let (HI) and (H^) hold. //(xo,Xi) € D{A) x D{A°‘), 0 < a < 1, 
then for any to E (0, T] 

lim sup {\\un{t) - u^(t)l|i + \\vn{t) - u„(t)|lQ} = 0. 

{n>m, to<t<T} 


Proof For n > m, we have 


\\fn{i)'^n)'^n) fTn{t-,Ura,Vm)\\ 

— ll/n(^) ^n) /n(^) '^mi '^m) || "h il/n(^) ’^rrij '^m) fm(l) ‘^mi ’^m) || 

< \\f{t, P^Unit), PXW) - /(t, P^Um{t),P^Vmm\ 

+ l|/(t, P^Vmit)) - fit, P^Umit), P^^Vmitrn 

< F^iTo)[\\P^Unit) - P^Umimi + IIPXW - P^vUt)\\a 
+ IKP'^ - PnUm{t)\\l + ||(P" - PnVmit)\U] 

< PR(^o)[llWn(0 —Umit)\\i + ll't^n(i) “ 'ymWlla 
+ 1|(P" - P'")n^(t)l|i + |1(P" - Pnvm{t)U 


Also, 


ll(P"-p-K(t)lU 


= p«(p"-p>„(t)ll 

= \\A^-'^{P^-P^)A%mit)\\ 

m 

- ’ 


where a < 1 } < 1. Thus, we have 


||/n(^5 '^n) 


/m(^j '^m) II 

< L(P)[l|'Un(^) “ + ll^n(i) — '^mit)\\a 


+ |l(P"-P”^)u^(i)l|l + 


\i9-q: ■* 
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For convenience, we denote 




For 0 <t'o< to, we may write 


^m,n{t) ^ 



e (‘ - I\\\\fn(s,Un,Vn) - fm{s,Um,Vm)\\ds 

\e~^^~^'>^A°‘\\\\fn{s,Un,Vn) “ fm{s,Ura,Vm)\\ds. 


We estimate the first and the third integrals as 


\e - I\\\\fn{s,Un,Vn) - fmis,Um,Vm)\\ds < 2L{R){M + l)t'o, 





v) - /m(s, Wm, •ym)|M5 < 2L{R)Ca{to - t'o) H'q. 


We estimate the second integral as 


-{t—s)A 


•^11 yTn(^) ^m) 




< (M + l)L(R) J‘ (^Un(.s) + ||(P“ - P”)«„{s)lli + ■ 

< {M + 1)L{R)( sup ||{P"-P”M«)llir+T^T+ 

\t'Q<s<T Jt’c. 

< (M + l)L(i?) i^rn^nT ^m,n(s)'j ) 


Bm,n= sup ||(F"-P”‘)«m(s)||l + 


t'o<s<T 


Vf 

H 

X^-a' 


where 
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For the fourth integral, we have the estimate 


lie "^^A'^\\\\fn{s,Un,Vn) - fmis,Um,Vra)\\ds 

< CaL(R) (t - s)-“ + i|(P" - P"‘)u^(5)||i + 

< C„i(fl) j d3. 


X'd-a 


ds 


Therefore, we have 


U,n{i) < 2L(P)[(M + 1) + C^{to - + C{R,T)Bm,n 

+ L(R) J\(M + l) + -^]UAs)ds 
< 2L(P)[(M + 1) + C^{to - 1' )-“]To + C{R, T)Bm,n 
+ Ni{T'^ + l) 


where 


and 


C{R,T) = L(R) 


{M+1)T + Ca 


rj'l-a ■ 

1 — a 


Ni = max{(M + l),Ca}L{R). 

Hence, from Lemma 1.3.27, there exists a constant K such that 


U,nii) < {2L(P)[(M+l) + C',(to-g-1t^ + (7(P,T)P,„,„}pr. 


Bm,n — > 0 as m — » oo provided IKP'^ — — > 0 as m — ^ oo for 0 < t < T, we 

get 


lim sup 

{n>m, to<t < T} 


^ 2L{R)[{M + 1) + Ca(to “ ^o) 
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As fo is axbitraxy, we observe that the right hand side may be made as small as desired 
by taking sufficiently small. 

Now, we show that for 0 < t < T, ||(P- - ^ 0 as m oo. We can 

easily check that ioi x E H and r] <0, 

WA^iP" - P’^)x\\ < AS.||(P”-P’»)i||<AJ.||j||. ( 3 , 4 . 2 ) 

From equation (3.4.1), we get 

||4l{P”-P’“)a„(t)ll < ll(i’"-P”)Aco|| + (M + l)||(P"-P'")i,|l 

+ (M + l)f \\{P" - F")Us,u^,Vr,)\\ds. ( 3 . 4 . 3 ) 

J 0 

We note that in (3.2.1), f(t,u,v) = f^(t,u,v) — Bu where represents original / 
appears in (3.1.2). Also, we can check ||Pu|| < C||Au|| for u 6 D[A) and some constant 
C. Therefore, we have 

ll(^”-0/m(s,U^,U^)|| < \\{P'^-P^)fZ{s,Um,Vrr:)\\ 

+ ||A(P--P>^(5)||. (3.4.4) 

Since A^f^{s, Um, Vm) e H, hence from (3.4.2), we have 

H{P" - P")fS.(s,Vr.,v„)\\ = \\A-I>(P" - P"')A^ 

Atti 

< AhiTo). ( 3 . 4 . 5 ) 

Am 

Using (5.4.10) and (3.4.5) in (3.4.3), we get 

ll(^"-J^^)«m(i)||i < \\iP^-P'^)Axo\\ + {M+l)[\\{P^-P^)xi\\ 

+ ArACro) + f‘ ll(P” - P”)ti.„(s)||i*]. 

Am Jo 
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Gronwall’s inequality implies that 

||(f>“-P"')u„(()||i < (||{P"-P'").4i„|j+(Af+l)[||{P”-P'»)a;,|| 

+ ^rPji{T„)l)e<"+»’' 

Am 

which tend to zero as m -> oo for 0 < t < T. This completes the proof of the 
proposition. 

Corollary 3.4.2 Let (HI) and (H2! ) hold. If {xo,Xi) 6 D{A) x D{A), then 
Urn sup {||u„(f) - Um{t)\\i + \\Vn{t) - Um{t)lla} = 0. 

Proof. Proposition 3.3.2 and 3.3.3 imply that in the proof of Proposition 3.4.1, we may 
take to = 0- 

For the convergence of the solution {un{t),Vn{t)) of the pair of approximate equa- 
tion (3.4.1), we have the following result. 

Theorem 3.4.3 Let (HI) and (HS! ) hold and let (a:o, 2 ;i) G D{A) x D(A°‘). Then 
there exists function (u,v) G Xi(T) x Xa{T) such that {un,Vn) —>■ {u,v) as n oo in 
Xi{T)xXa{T) and {u,v) satisfies (3.3.1) on [0,T]. Furthermore, {u,v) can be extended 
to the maximal interval of existence [0,tmax)) 0 < tmax < oo satisfying (3.3.1) and {u, v) 
is a unique solution to (3.3.1 ) on [0, tmo®)- 

Proof First, we assume (xq, Xi) G D{A) x D{A). Corollary 3.4.2 imphes that there 
exists {u,v) G Xi{T) x Xa(T) such that converges to {u,v) in XfiT) x XdT). 
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Since ^ each ti, (u^v') is also in Further, we have 

\\fn{t,Un,Vn) “ f {t , u{t) , V {t))\\ 

+ \\f{t,PMi),PMt)) - 

^ FniToMP'^Unit) - P"n(i)||i + ||PXW - P"t;(t)||a 

+ ll(P"-/Mt)||l + ||(P--/)z;(t)y 

< ^«(2o)[||'lin(t) - u{t)\\i + ||z;„(t) - y(t)||^ 

+ ||(P”-/)u(t)||i + ||(p--j)^{t)||j. 

Taking supremum over [0,T], we get 

sup \\fn{t,Un,Vn) “ f {t, u{t),v(t))\\ 

0<t<T 

— ■^.R(^o)[||Un — u||xi(r) + ||u„ — v\\x^{T) 

+ IKP" - n«IU.(7-) + ||(P“ - /)»|U.(7-)], (3.4.6) 

and the right hand side of (3.4.6) tends to zero as n tend to infinity. Using (3.4.6) and 
the bounded convergence theorem in (3.4.1) , we get 

u(t) = xo- (e-‘^ - /)(44)-^a:i - [ - I){A)-^f{s, u(s),v{s))ds, 

Jo 

v{t) = e-‘V+/ e-^^-^^^f(s,u{s),v{s))ds. 

Jo 

Thus, for (a:o)a;i) G D{A) x D{A), there exists a unique pair of functions {u,v) G 
^i(T) X Xa{T) such that (u„,Un) u) as n oo in Xi{T) x ^^(T) and (u,t;) 

satisfies (3.3.1) on [0,T]. 
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Now, let us assume that (a;o,2;i) € D{A) x D{A°‘). Since for 0 < t < T, 
{Aun{t),A^Vnit)) {Au{t),A^v{t)) and also (un(0),t;n(0)) = (u(0),u(0)) = (xo,xi) 
we have, for 0 < t < T, {Aunit), A'^Vnit)) — > (^u(t),.A‘^u(t)) in H. Furthermore, since 
each (un, ^n) is in Wr, we have (u, v) G Wr and for any 0 < to < T, 

lim sup {llun(t)-u(t)||i + |lun(t)-u(t)lU = 0. 

oo fo<i<T 

To prove {u,v) satisfy integral equation (3.3.1), it is sufficient to show that v satisfies 
second integral equation of (3.3.1). We have 


sup \\fn{t,Un,Vn)- f{t,u{t),v{t))\\ < F^(To)[|lu„ - ulU,(T) + K “ t^ll;f„(T) 

to<t<T 

+ IKP“-J'Mk(T) + ll(P”-J>lk(T)l 


and the right hand side of the above inequality tends to zero as n tends to infinity. 
Now, for 0 < to < t, we may rewrite the second integral equation of (3.4.1) as 


The first integral can be estimated as 

rto 


e ")^/„(s,u„,u„)ds 


< ML{R)to. 


Therefore, 


Vn{t)-e ^ ^^'^ fn{s,Un,Vn)ds = 6 ®^^/n(s, «n, Wn)di 

< ML{R)tQ. 


< ML{R)to. 


Letting n -h- oo in the above inequality, we get 

v(t) — e~^^xx - f eA^~^'>-^f(s,u{s),v{s))ds 
Jto 

Since 0 < to ^ T is arbitrary, we get that v satisfies the second integral equation of 
(3.3.1). 
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If (w, v) satisfy (3.3.1) on [0,T], then we show that (u,v) can be extended further. 
Since 0 < To < oo was arbitrary, we assume that 0 < T < To. We consider the equation 

U"(t) + AU'{t) = F{t, U{t), U'it)), t > 0, 
f/(0)=u(T), U'{0)=u'{n 

which can also be written as system of equations 


U'(t) = V(t), U(0) = n{n 
V'{t) + AV{t) = F(t,U(t),V{t)), V{0) = v{T). 


where F : [0,To — T] x Xi x Xa H are defined by 

F{t,x,x) = f{t + T,x, x), 

for (t, X, x) € [0, To — T] X Wi X Xa- F satisfies (H2) for To replaced by To — T. Hence, 
there exists a (JJ, V) G (^([O, Ti], X\ x Xa), for some 0 < Ti < To — T, satisfying the 
integral equation 


U{t) = u{T)F{e-^^-I){-A)-\{T) 

+ / - I){-A)-^F{s, U{s), V{s))ds, 

Jo 

V{t) = e-^\{T)+ [ e-^^-^^^F{s,U{s),Vis))ds, ' 

Jo 

with 0 < t < Ti. Now, we define 




{u{t),v{t)), 0<t<T, 

{U{t-T),V{t-T)), T<t<Ti + T. 
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Then, for t e [0, Ti + T], (u, v) satisfies the integral equations 

u{t) = XQ-\-{e 

+ f\e-C-’}^-I){-A)-^f(s,u(s),v{s))ds, 

Jo 

(3.4.7) 

v{t) = + f e~^^~^'>^f{s,u{s),v{s))ds. 

Jo 

To see this, we need to verify (3.4.7) only on [T, Ti + T]. We have 


(t) = U{t-T) 

= u{T) + - I){-A)-^v{T) 

+ [ - r){-A)-^F{s, U{s), V{s))ds 

Jo 

= Xo- - I){A)~'^xi - f - I){A)''^f{s, u{s),v{s))ds 

Jo 

— — /)(yl)“^ + J e~^'^~^^^f{s,u{s),v{s))ds 


rt-T 


(e' 


-(t-T-s)A 


I){A)-^F{s,U{s),V{s))ds. 


Putting T + s = rj, we have 


u{t) = Xo- (e - I){A)~'^xi- f {e~^^ - I){A) ^f{s,u{s),v{s))ds 

Jo 


/; 


„-(t-s)A _ 


I){A)-'^F{s - T, U{s - T),Vis - T))ds. 


Since F{s - T, U{s - T), l/(s - T)) = f{s, U{s - T), V{s - T)), we have 

ft(t) = a;o-(e“‘^-J)(A)“^ 2 :i- f - I)iA)~'^f{s,u{s),v{s))ds. 

Jo 
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Since 


\\f{s,Ui{s),Vi{s)) - /(s,U2(s),W2(s))i| 

^ F nit max )[ll^l('5) - W2(s)||i + 11^1 (s) -n2(s)|y, 


we get 


\\ui{t) - ^ 2(0111 


+ |lt>i(t) - y2(t)llQ 

< {M + l)Fii{tmax) [ll-iii - ^211x1(5) + ini - •y 2 llx„(s)]cis 

+ Fft,{t„^ax)Ca J [t — s) ^'[llui — ti2||xi(s) + 11^1 — n 2 ||xa(s)]i ^5 

< N2{T° ^ / {t - s)°‘ ” '“ 2 |Ui(s) Hr Ibi - n2llx„(5)](^S, 


where 

N2 = max{{M + l),Ca}Fii{tmax)- 

Hence, from Lemma 1.3.27, we get 

\\ui{t) - U 2 (t)||i + ||ni(i) - ■U 2 (t)iU = 0) 

for all 0 < t < T. Prom the facts that 

IK(i)-W2WI| < -«2(t)l|i 

Ao 

and 

11^1^-^2(011 < •4llt^i(0-''^2(0iU. 

it follows that = ( 112 , 112 ) on [0,r]. Since 0 <T < tmax was arbitrary, we have 

(ui,i;i) = (^^ 2 , 112 ) on [0, tmax)- This completes the proof of the theorem. 
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3.5 Faedo-Galerkin approximations 


For any 0 < T < tmax, we have a unique pair {u,v) G XiiT) x Xa{T) satisfying the 
integral equations 


u{t) — xq (e — [ (e — /)(7l) ^f{$,u{s),v{s))ds, 

Jo 


v{t) = e [ e "^^f{s,u{s),v(s))ds. 

Jo 


(3.5.1) 


Also, we have a unique solution (K„,ti„) € Xi{T) x Xc{T) of the approximate integral 
equations 


u. 


(t) = xo-{e-^^-I)A-^x^- [\e-<-^-^^^-l)A-^f{s,P^Ur,{s),P^Vr,{s))ds, 

Jo 

(3.5.2) 

Vnit) = e~^^Xi+ [ e“^*"*)^/(s,P"un(5),PX(s))c?s- 

Jo 


If we project the equation (3.5.2) onto Hn, we get the Faedo-Galerkin approximations 

(un(t),Vn(t)) = (P’^Un(t), P’^Vn(t)) Satisfying 

Unit) = P^Xo - (e-‘^ - /)s4-'P"Xi - - I)A-^P^f(s, Un{s),Vn{s))ds, 

Jo 

(3.5.3) 

Vn{t) = e“‘^P"Xi-f /" e“^‘“^^^P"/(s.*n(5),fy(5))(^'S- 

^0 

The solution (u,u) of (3.5.1) and {un,Vn) of (3.5.3), have the representation 

OO 

uit) = '^aiit)ui, Q:t(t) = (n(t),Ui), i = 0, 1,...; 


1=0 


’(t) = '^Piit)ui, Piit) = iv{t),Ui), i = 0, 1, 


2=0 


(3.5.4) 
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and 


u. 


,(t) = o:'y{t) = {un{t),Ui), z = 0,l,...,n; 


2=0 


(3.5.5) 




2=0 


Using (3.5.5) in (3.5.3), we obtain the following system of first order integro-differential 
and differential equations, 


daji't) 

dt 

dm) 

dt 


+ 


JO 

(3.5.6) 

+XiPm = 


with the initial conditions 


af{0) = (f>i, Pi{0) = ^i, 


where 




Xli 


2=0 


2=0 


and (j)i = (xo, Ui), ijji = (a;i, Ui) for z = 1, 2, n. 

The systems (3.5.6) determine the Q:"(i)’s and /5f(t)’s. Now, we shall show the 
convergence of (af, pf-) to (a, P). It can be easily checked that 


A[u{t) — u(t)] = A 


^(ai(t) - a^(t))ui 


, 2=0 


^Ai(Q:i(t) -a"(t))ui 


2=0 
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and 


OO 

.i=0 

= T,Km)-m)w 

i=0 


A°‘[v{t) - v{t)] = 


Thus, we have 

MK*) - flWir > X^A?(ai(()-a"(f))^ 


i=0 

We have the following convergence Theorem. 


Theorem 3.5.1 Let (HI) and (H2) hold. Then, we have the following. 


(a) If {xq, Si) € D{A) X D{A°‘), then for any 0 <to <T, 


lim sup < f = 0- 

n^<=Oto<t<T ^ 1 


(b) If {xo,Xi) G D{A) X D{A), then 


lim sup f = 0. 

n-*ooo<t<r 1 ^ ^ I 


The assertion of Theorem 3.5.1 follows from the facts mentioned above and from 
the following proposition. 


Propostion 3.5.2 Let (HI ) and (H2) hold and let T be any number such that 0 < T < 
imax, then we have the following. 

(a) If {xo,xi) G D{A) X D{A°‘), then for any 0 <to <T, 

lim sup {l|A(u„(t) - Umimi + = 0. 

{n>m, to<t<T} 
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(b) If {xo, a:i) G D{A) x D{A), then 

— - ««.(0)ll} = 0 . 

Proof. For n > m, we have 


- Um{t))\\ + |iA“(t)n(t) -«m(i))i| 

= \\A(P^Un{t) - P^Umirn + - P^Vmitm 

< \\AP^{Un{t) - li^(t))|| + \\A{P^ - P^)Umit)\\ 

+ \\A^P^(Vn{t) - Vm{t))\\ + \\A^{P^ - P^)Vm{t)\\ 

< \\Unit) - U^(t)||i + ||t;„(t) - 

+ ||(P"-P’”)u„(t)l|. + ^PVil 

If (xo,Xi) 6 D{A) X D{A°‘), then the result in (a) follows directly from Proposition 
3.4.1. If {xo,xi) E D(A) X D(A), then (b) follows from Corollary 3.4.2. 



Chapter 4 


Second Order semi-linear 
integro-differential Equations - I 


Obviousness is always the enemy of correctness. 
Hence we must invent a new and 
difficult symbolism in which nothing 
is obvious 


— Novalis 


4.1 Introduction 

In chapters 2 and 3, we have established the results for the approximation of the solution 
of a Sobolev type evolution equation and a second order semi-linear evolution equation, 
respectively. For both kind of evolution equations, the existence and uniqueness results 
had been established by Hernandez [39] and Bahuguna [6], respectively. Now, we are 
interested to establish similar results for a class of second order semi-linear integro- 
differential equation. For such kind of evolution equation the existence and uniqueness 
results have not been proved. In this chapter, therefore, we shah establish the existence 
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and uniqueness of a solution for an abstract second order semi-linear integro-differential 
equation in a Banach space. Such type of equations arise in the study of viscoelastic 
material with memory. 

Consider the following initial boundary value problem for the strongly damped 
integro-differential equation 

■+ (aL H- 5/) {x,t) -f- {cL + dl)u{x,t) = h{x,t,u{x,t),^{x,t)) 

du 

-f / k{t - s)h{x,s,u{x,s),—{x,s))ds, (x,t) € Q x (0,T), 

du 

u{x,to) = Xq{x), —{x,to) = Xi{x), XEQ, (4.1.1) 

where be a bounded domain in with sufficiently smooth boundary dfl and Lu = 
^ ^ symmetric second order strongly elliptic differential operator 

in ri. Here a > 0, b, c, d are constants and h and g are smooth nonlinear functions and 
fc is a locally p-integrable function for 1 < p < oo. 

We can treat (4.1.1), as a special case of the following initial value problem for 
abstract second order semi-linear integro-differential equation in a Banach space X 

= / (t, M(t), ^{t)) + Kt- s)g{s, u{s),.~{s))ds, t > to, 

u{to) = xq, u'(to) = xi. (4.1.2) 

For (4.1.1) we take A = aL, B = (cL + dl) and f[t,u,v) = h{t,u,v) — bv. We assume 
that —A generates an analytic semigroup T(t) in X. The nonlinear map / and g satisfy 
Assumption (F) and Assumption (G), respectively, and the kernel k satisfies (K) stated 
in the next section. 
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In this chapter, we concentrate on the study of the equation 

u"{t) + Au'{t) = f{t,u{t),u'{t)) + f k{t — s)g{s,u{s),u'{s))ds, t>to, 

Jto 

u(to) — iTo, u (to) ^ii (4.1.3) 

as we can merge Bu with / so that / still satisfied (F). 

Duvaut and Lions [31], Glowinski, Lions and Tremolieres [34] have studied par- 
ticular case of (4.1.1) in which L = -A and A: = 0, in the context of the theorem of 
viscoelastic materials. Sandefur [84] studied the second order semi-linear differential 
equation 


u''{t) -1- Au'{t) -f Bu{t) = f{t,u{t)) 

u{0) = (t)y u'{Q) = -0, (4.1.4) 

in a Banach space X under the assumptions that the linear operators A and B can 
be decomposed as —A = Ai + A 2 and B = A 2 A 1 , where A^ generates a co-semigroup 
Tk{t), k = 1,2 ; and / satisfies a locally Lipschitz condition. He established the local 
existence and uniqueness of a mild solution to (4.1.4), i.e., there exists a continuous 
function u on [0, c] for some c > 0 such that u satisfies the integral equation 

u{t) = Ti(t)^ -f- 

+ [ [ Ti(t-T)T2it - s)f{s,u{s))dsdT, 

Jo Jo 

where 0 € D(Ai). Aviles and Sandefur [3] studied the well-posedness of (4.1.4) under 
similar conditions. 


/ 


Ti{t - T)r2(r)(V’ - Ai(l))dT 


Engler, Neubrander and Sandefur [32] proved the local existence and uniqueness 
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of a mild solution to special case of (4.1.3) in which /c = 0, that is, 
u"{t) + Au'{t) = f{t,u{t),u'{t)), 

u{to) = a:o, u'(t) = Xi, ( 4 .I. 5 ) 

under the assumptions that —A generates an anal 3 dic semigroup T(t) in X and / 
satisfies a condition similar to Assumption (F), where a mild solution on [to,ti), for 
some ti > to, to (4.1.5) is the first component of a solution {u{t),v{t)) of the integral 
equations 

u{t) = xo + {T{t-to) - I){-A)~^xi 

+ f {T(t- s) - I){-A)~^f{s,u{s),v{s))ds, to<t<ti, 

Jto 

v{t) = T{t — to)xi+ / T{t — s)f{s,u{s),v{s))ds, tQ<t<ti. 

Jto 


Bahuguna [ 6 ] improved the results of [32], by showing that (4.1.5) has a unique 
classical solution locally, that is, there exists a unique u e C^{[to,ti) : X) n C^{{to,ti) : 
X) and satisfies (4.1.5) on [toi^i) for some ti > to- Further, he has discussed the 
continuation of this solution, maximal interval of existence and the global existence. 

In this chapter, we show that (4.1.3) has a unique classical solution locally, i.e., 
there exists a unique u € C^([to, h) : X) D C^{{to,ti) : X) and satisfies (4.1.3) on [to, ti) 
for some ti > to. Further, we discuss the continuation of this solution, maximal interval 
of existence and the global existence. We achieve these objectives by extending the ideas 
and techniques used in the proof of Theorems 6.3.1 and 6.3.3 in Pazy [77] concerning 
the semi-linear equations of the first order to (4.1.3). For the global existence, we 
require the modified version of Lemma 4 . 4 . 2 , stated and proved in fourth Section which 
is originally stated in Pazy [ 77 ] as Lemma 5.6.7. 
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4.2 Preliminaries and Assumptions 

Let X be a Banach space and let -A generate the analytic semigroup T{t) in X. we 
note that if -A is the infinitesimal generator of an analytic semigroup then -{A + al) 
is invertible and generates a bounded analytic semigroup for a > 0 large enough. This 
allow us to reduce the general case, in which -A is the infinitesimal generator of an 
analytic semigroup, to the case where the semigroup is bounded and the generator is 
invertible. Hence, for convenience, without loss of generality, we assume that T{t) is 
bounded, that is ||T(t)|| < M for t > 0 and 0 € p(-A), i.e., -A is invertible. Here 
p(-A) is the resolvent set of -A. It follows that, for 0 < a < 1, A“ can be defined as a 
closed linear invertible operator with its domain D(A^) being dense in X. We denote 
by Xc, the Banach space £»(A“) equipped with the norm 

l|2^IU = ||A“a:|| 

which is equivalent to the graph norm of A-. For 0 < a < p, we have X^ C X, and 
the embedding is continuous. 

We consider the problem 

u"it) + Au\t) = f{t,u{t),u'it))+f k{t-s)9it,u{t),u'{t))ds, t>to, 

J to 

w(to) = Xq, u'(io) = Xi. (4.2.1) 

On the kernel k we assume the following condition: 

(K) Kernel k e oo) for some 1 < p < oo is locally Holder continuous on (0, oo), 

i.e., 

|A:(t) - A;(s)| < foj. 5,ie(0, oo) and 0 < p < 1. 

The nonlinear functions / and g are assumed to satisfy the following assumptions: 
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Let U be an open set in R+ x Xi x Xa- 

Assumption (F): A function / said to satisfy the Assumption (F), if for every {t, x, x) € 
U, there exists a neighborhood V C U and constant L/>0, 0<i?<l, such that 

\\f{ti,xi,xi) - f{t2,X2,X2)\\ < Lf[\ti - t2\^ + ||xi - a;2l|i + llxi - X2tlQ,], (4.2.2) 


for all {ti, XiXi) € V. 


Assumption (G): A function g said to satisfy the Assumption (G), if for every 
{t,x,x) € U, there exists a neighborhood V C U and a nonnegative function Lg € 
LfociO, oo), where l<g<oo, ^ + ^ = 1, such that 

\\g(t,xi,xi) - g{t,X 2 ,X 2 )\\ < Lg{t)[\\xi - X 2 \\i + \\xi ~ X 2 \\a], (4.2.3) 


for all {ti,XiXi) G V. 

By a local classical solution to (4.2.1), we mean a function u G (^^([io,^!) : X) fl 
: A), satisfying (4.2.1) on [toiii), for some > to. By a local mild solution 
to (4.2.1) on [toi^i)) for some ti > to, we mean the first component of a solution {u,v) 
to integral equations 


u{t) 


v{t) 


xo + {T{t - to) - I){-A) ^xi+ / {T{t - s) - I){-A) ^[f{s,u{s),v{s)) 

JtQ 

+ / k{s — T)g{r^u{T)^v{T))dT]ds^ tQ<t<ti, 

JtQ 

T{t-to)xi+ / T{t- s)[f{s,u{s),v{s)) 

JtQ 

+ I k{s — T)g{T^u{r),v{r))dT]ds^ to<t<ti. (4.2.4) 

Jtcs 
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4.3 Local Existence of Solutions 


As we have already pointed, without loss of generality, that the semigroup generated by 
-A can be assumed to be bounded and A is invertible. Under these conditions imposed 
on A, we prove the following local existence and uniqueness theorem. 

Theorem 4.3.1 Suppose that -A generates the analytic semigroup T{t) such that 
||r(t)|| < M and 0 e p(-A). If the maps f and g satisfy Assumption (F) and As- 
sumption (G), respectively, and kernel k satisfy (K), then (4.2.1) has a unique local 
classical solution. 


Proof. Fix (to,Xo,Xi) m U and choose F, > and 5 > 0 such that (4.2.2), with 

some fixed constant L/ > 0, 0 < < 1 and (4.2.3) with nonnegative function L„(t) 

hold on the set 

V' = {(t,l,i)€£f ||x-ao||, + p-ii||„ < {}, 

Let 

and 

C(5) — [Lf + -i-Bf-f- Bg\\k\\u,^t^ i>^){t[ - to)i. 

Choose ti > Iq such that 

lir(t-to)a:i -a:i|| + ||r(t-to)A“2;i -A“2 ;i|| < - 

3 

and 
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where is a positive constant depending on a and satisfying 


\\A°‘Tit)\\ < for t>0. 


(4.3.1) 


Let Y - C{[to, ti\:XxX). Then, y G T is of the form y = (y^, y. ^ c{[to, ti] : X), 
z = 1, 2. Y, endowed with the supremum norm, 


\{yuy2)\\Y= sup [|lyi(t)|| + ||y2(t)||], 


is a Banach space. We define a map F on y by Fy = F(yi,y 2 ) := (y^, y^) with 


yi{i) 

y2{t) 


= Axo - {T{t - to) - I)xi - f\T{t -s)- I)Fy{s)ds, 

JtQ 

= Tit~to)A‘^xi+ [ T{t-s)A°^Fy{s)ds, 

•/to 


(4.3.2) 


where 

^yii) = f{i,A ^yi{t),A~°‘y2{t)) + f k{t-T)g(r,A~\i{T),A~°'y2{T))dT, 

J to 

for t G [io 5 ^] • 


For every y eY, Fy{to) — {Axo,A°‘xi). Now, we show that the map 1 Fy(t) is 
continuous function from [to, ti] to X x X. For s<t, 

Fy{t)-Fy{s) = (yi(t),y 2 (/))-(yi(s),y 2 (s)) 

= (yi(/) - 1/1(5), y2(t) - y2{s)). 
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Taking norm on both side, we have 

||Fy(t)-T’?/(s)|| = ||yi(i) -yi(s)|| + ||y2 (t) -y 2(s))|| 

< mt-to)-T(s-to)|||Ia;i|l 

+ f ll^(^-^)-r(s-r)||||Fj,(T)||dr 

to 

+ l P(*-^)--f|IK(r)||dT 

+ ll^(^ - to) - T{s - to)||||A“xi|| 

+ / -r)- T{s - r))^“||||f;(r)|ldT 

J to 

+ l Wnt - r)A‘‘\\\\F,(T)\\dT. 


We have 


< \\fit,A-^y^{t),A-°‘y2{t)) - /(t,Xo,Xi)|| + ||/(i,Xo, Xi)|| 

< Lf[\\y^{t) - Ax4 + \\y^{t) - A^xi\\] + Bf 

< Lf sup [||yi(t)-yla;o|| + ||y2(t)-^“Xi||] + 5;. 


Mt,A-\,{t),A-‘^y 2 {t))\\ 

< Lg{t) sup ^ [||yi(t) - AxoW + \\y 2 {t) - A“xi||] + Bg. 


(4.3.3) 
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Therefore, for Fy{t) we have the estimate 

\\Fy{t)\\ = \\f(t,A-^yi{t),A~^y 2 {t))\\ 

+ / \k{t-r)\\\g{T,A-^y^{T),A-^y2{T))\\dT 

to 

< (-^'Z + P||LP(Jo,t')||-f'gliL9(to,t'i)) sup [||yi(t) - Axo\\ + \\y2(t) -A^Xill] 
+ (Bg + Bg\\k\\Lp^to,t[){t'i - := Ni, 


since ?/i, j /2 S ^([to.ii] : X). Thus, we have estimates for the second and the fourth 
integrals 

J ll^(^ ~ ~ -^||||-^j/('r)IMT < {M l)Ni[t — s), (4.3.4) 

j\\T(t-r)A’'\\\\Fy(T)\\dT < (4.35) 

Part (d) of Theorem 1.3.19 implies that for 0 < ^ < 1 and x 6 D{A^), we have 

ilTO-/)x|| < C'/\\x\\,. 

If 0 < /? < 1 IS such that 0 < a + /5 < 1, then A°‘y e D{A^) for any y 6 D{A°^+^). 
Therefore, for t,s € [0,r], we have 

||(T(t)-/)A“T(s)rr|| < C't^p“r(s)a:|l^ = C'^t^p“+^T(s)a;|| 

< (4.3.7) 

From (4.3.7), we get 

||(r(e-r)-r(s-r))X“|| = ||(T((-s)-/)yl“r((-7-)|| 



4.3 Local Existence of Solutions 


113 


Hence 


r II (r(* - 

Jto 


- T(s - r))A'>\\\\F,(r)\\dr 


< 


NxC'pCo^^^(t-sf 


1 - (a + ^) 


Also, from (4.3.6), we have 


(4.3.8) 


||{r(f) - /)T(s)l|| < C),t«||T(s)i||j = C;i'’l|A'’r(s)i|| 
< C'i,CgtOs-’’\\x\\. 


Therefore 

lir(t - T) - r(s - r)|| = ||(T((-s)-/)r(s-r)|| 

< c'fi0(t-sf(s-TrK 

Hence 

rilT(f-T) - T(s-T)ll||F,(r)|MT 

JtQ 

< (4.3.9) 

Inequalities (4.3.3), (4.3.4), (4.3.5), (4.3.8) and (4.3.9) implies that F :Y —^Y. Let S 
be a nonempty closed and bounded set given by 

5 = {j/ € y I y = iyi,y2),yi{to) = Axo,y2(io) = A°‘xu ||yi(t)-Aa;o||+l|t/2(t)-44“2:i|| < 5}. 
Let y = {yi,y 2 ) be any element of S. Prom (4.3.2), we have 

l|yi(0 ~ 44a:ol| + ||y2('^) “ 44“a;i|| 

< ll(ni-^o)-/)a:i||+ f\\T{t-s)-I\\\\Fyis)\\ds 

Jto 

+ \\{T{t - to) - /)A“ii|| + f ||A“r(t - s)||||F„(s)||ds. 

Jto 


(4.3.10) 
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To find the estimate for Fy{s), adding and subtracting fis,xo,x^) and g{s,xo,Xi) and 
using (F), (G) and (K), we obtain 


||^y(s)ii < 11/(5,^ “2/2(5)) -/(s,a;o,rci) II + B/ 

< [Lf + PllLP(to,ii)ll^5llL.(fo,f;)]^ + + SJ|A:|Up(t„,j,)(t' - to)-^ 


< C{6). 


(4.3.11) 


estimate (4.3.11) and the fact that ||T(t)|| < M and (4.3.1), we get from (4.3.10) 


||yi(t) — x42;o|| + ||y2(t) — A“a:i|| 


< 

~ 3 
< 5. 


{M + l)C'(5)(t - to) + P^*^(^)(^-^o)^~ 

1 - a 


Hence, F : S S. Now, we show F is a contraction on S. Let (yi,y 2 ) and (zi, z^) be 
any two point of S. From (4.3.2), we have 

Il2/i(i) — i'i(t)|| + \\y2{t) — Z2{t)\\ 

< / l|r((-s)-/||||Ji(s)-nW||* 

•j to 

+ f ||r(t - s)/l“||||F,(s) - f;(s)||*. (4.3.12) 

tQ 

Using (F), (G) and (K), we get 
ll^y(5) - F,(5)|| 

< \\fis,A-^y^{s),A~^y2(s)) - f{s,A-h,(s),A-‘^Z2{s))\\ 

+ / |a(s - T|||^(r,24-Vi(r), A-“y 2 (r)) - g(r,A-^zi{T), A-^Z 2 {T))\\dT 

to 

< [Lf + P||LP(to,t/)||L5||i,(t„,t/j]||(2/i,y2) - (2^i,^2)||y 
^ ~^\\{yuy2)-{zi,z2)\\Y- 


(4.3.13) 
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Using (4.3.13) in (4.3.12), we get 


llSi(t) - + 

< 

< 


\\y2ii)-Mi)\\ 

r{M + l)C'(<5)(t-to) aC(5)(t-to)'-“l 

5 6{l - a) 

2 

■^\\iyi,y2) - {zi,z2)\\Y- 


ll(yi)?/ 2 ) — (2:1,22) 11 k 


Taking supremum over we have 

\\{yi,y 2 )-(zuZ 2 )\\Y < |||(?/l,y2) - (2l,22)l|y. 

Thus, F is a contraction on S. Therefore, it has a unique fixed point in S. Let 
y = ( 2 / 1 , 2 / 2 ) e 5 be that fixed point of F. Then, 

yi{t) = Axo-{T{t-to)-I)xi- f {T{t - s) - I)Fg{s)ds, 

Jto 

Mt) = T{t-to)A^Xi+ [ T{t- s)A^Fy{s)ds, (4.3.14) 

Jto 

where 

P'yii) = /(^i ^~Vi(0) ^~“y 2 (i)) + [ k(t- T)g{T, A~^yi{r),A~°'y 2 {T))dr. 

Jto 

We note that {u,v) = {A~^yi, A''°‘y 2 ) is the unique solution of the integral equations 
(4.2.4) on [toj^i]- We can easily check that Assumption (F) and the continuity of yi 
and y 2 on [to, ti] imply that the map t i— »• Fy{t) is continuous and hence bounded on 
[toj^i]- Let |lFy(t)|| < N fox to < t < ti. We wiU now show that t Fy{t) is locally 
Hdlder continuous on (to,ii]- For this, we first show that yi and y 2 are locally Holder 
continuous on (to, ti]. Prom Theorem 2.6.13 in Pazy [77], for every 0 < P < 1 — a and 
every 0 < /i< 1, we have 

l|(T(/i)-/)A‘^r(t-s)|| < C/j/i^||A“+^r(t-s)|| 

< C/i^(t - 


(4.3.15) 
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Now 

||S2(t + ft)-ff2(<)ll < ||(r(fe) - 

+ [‘ \mh) - I)A-T{t - s)f 5 (s)|lds 

Jto 

/ t+/l 

||.i 4 T(t + h — s)Fy(s)||(is = Ii -j- 12 -{■ /a- 

We use ( 4 . 3 . 15 ) to get 


h < < Mih^, 

7. r\tAk-s)-^ = ^ 

Jt 1 “ 


h < NCa 


< M^hK 


Here, Mi depends on t and increases to infinity as t | toi while M2 and M3 can be 
chosen independent of t. From the above estimates, it follows that there exists a positive 
constant C such that for every tg > ^o, 

11^2(0 - ^2(5)11 < C\t - sf for to < to < t, 5 < ti. 

Similar result holds for yi (if we take a = 0 in above consideration). For s, t G (to,ti] 
with t > s, we have 


IIFjO) - F5(s)|| < ||/(t, A-'ftO). A-“»2(f)) - f{s, >i-“!r2(s))ll 

+ f \k{t - t) - a{a - T)\\\g{T,A-'-yi{T),A-'‘ff2(T))\\dT 

Jto 

+ f \k{t-r)\\\g{T,A~^yi{T),A~'^y2{T))\\dT. 

J $ 

Since k is holder continuous with exponent /i, we have 



\k{t-T) -k{s-r)\\\g{T,A ^yi{T),A “2/2('r))||dr 


< iV(ti-to)|t-sp 


( 4 . 3 . 16 ) 
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and 


[ \k(t-T)\\\g{r,A ^yi{T),A ^y2{T))\\dT < ArA:o(ti - to)“|t - (4.3.17) 

J S 

where ko = maxto<t<ti \k{t)\. The local Holder continuity of Fy(t) on (to,ti] follows 
from the Assumption (F), and the local Holder continuity of yi and y 2 on (to,ti] and 
from estimates (4.3.16) and (4.3.17). 

Consider the inhomogeneous initial value problem 


dv{t) 

dt 


+ Av{t) = Fy{t), 


vito) = Xi. 


(4.3.18) 


By the Corollary 1.3.20, this problem has a unique solution v 6 C^((to,fi] : X) given 
by 


v{t) = T{t-tQ)xi+ [ T(t ~ s)Fy{s)ds, (4.3.19) 

Jto 

for t > to- Each term on right hand side belongs to D{A) and hence belongs to D(A°‘) 
since D{A) C D{A°‘), 0 < a < 1. Operating on both sides of (4.3.19) with A°‘, we find 

A^^vit) = T(t - to)4l“a:i + f T{t - s)A‘^Fy{s)ds, (4.3.20) 

Jto 

By (4.3.14), the right hand side of (4.3.20) equals to y 2 {t) and therefore A°‘v{t) = y 2 {t) 
or v{t) = A~°‘y 2 {t). Let u{t) = A“^yi(t), then we have u(t) = xq + J^v{s)ds, which 
yields u{t) E C^([to,ti) : X) nC^((to,ti) : X). Thus, u satisfies (4.2.1) on [toi^i)- This 
completes the proof of Theorem 4.3.1. 
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4.4 Global Existence of Solutions 


In this section, we will prove under additional growth condition on the nonlinear map 
/ and g the following global existence result. 

Theorem 4.4.1 Let 0 6 D(-A) and -4 be the infinitesimal generator of an analytic 
semigroup T(t) such that ||r(t)|| < M fort>0. Let J,g : ( 0 , co) x x X. X satisfy 
Assumption (F) and Assumption (G), respectively, and k satisfies (K). If there exists a 
nondecreasing function a, : [to, oo) ev and a nonnegatwe function a, € i^(0, oo) 
where q is same as before, such that 

\\f{t,x,x)\\ < a/W[H- ||rc||i + pllj for t>tQ,{x,x)eXixXos, 

< a 5 (i)[H- |la:||i + for t>tQ,{x,x)eXixXce, 

then for each {xq, Xi) e Xi x X^, (4-2.1) has a unique classical solution u which exists 
for all t > to- 


For the proof of this theorem we shall require the following lemma. 

Lemma 4.4.2 Let (/){t,s) > 0 be continuous on 0 < s < t < T. If there are positive 
constants A, Sj and /3, such that 

s) < ^ + 0 (a, s)da + ^2 ^ (t - V(c 7 , s)da, ( 4 . 4 . l) 

for 0 < s < t < T, then there exists a positive constant C such that fit, s) < C for 
0 < s < t < r. 


Proof We have 


Is I ^fir,s)dTda = (^j (t - s)dr 


(4.4.2) 
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and a well-known identity 



- ry-^da = (t - 


i mm 

vifi+iY 


( 4 . 4 . 3 ) 


Iterating (4.4.1) n - 1 times, using (4.4.2) and (4.4.3) and majorating (t - s) and (t - t) 

by T, we get 



Choosing n sufficiently large so that n/5 > 1 and replacing {t — by we get 

Cl C 2 f (f){cr, s)do', 

J S 

where Ci and C 2 are the positive constants independent of s. The required result then 
follows from Gronwall’s inequality. 

Proof of Theorem f.f.l. Let [toi^) be the maximal interval of existence for the 
solution u to (4.2.1) guaranteed by Theorem 4.3.1. It is suffices to prove that 

im^)iii + ikWIU < C 

on [to, T) for some fixed constant C > 0 independent of t. 

Now, since u{t) is a solution of (4.2.1) on [to, T), it is also a mild solution to (4.2.1) 
therefore from (4.3.14), we have 

Au{t) = Axo - (r(t - to) - I)xi - f{T{t -s)- I)F{s)ds, 

Jto 

A^u\t) = r(t - tQ)A^xi + frit - s)A°‘F{s)ds, 

Jto 


(4.4.4) 
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where 


+ I k(t- T)g{T,u(T),u'(T))dT. 

J to 

From (4.4.4), we have 

[1 + ||u(7?)||i + ||«'(77)1|J = [1 + ||Am( 7?)|| + pV(7?)||] 

< '^ + \\Axo\\ + {M + l)\\xi\\ + [M + 1) f ||F(s)||(is 
+M ||a:i|U + a(77 - s)-^F{s)\\ds. ( 4 . 4 . 5 ) 

Our assumptions on /, g and k imply that 

l!^(s)!( < \\f{t,u(t),u'(t))\\+ f |/:(s-r)l||5(r,u(r),u'(r))||dr 

^ to 

< (af{T) + ||A:l|iP(to.r)||a3lU«(to,r))^ sup [1 + ||u(r)||i + i|u'(T)||^]. 


(4.4.6) 


Using (4.4.6) in (4.4.5), we get 


[1 + ||u(r/)||i H- \\u (r 7 )||Q,] < C1 + C2 f sup [1 + ||u(r)||i + \\u'(t 

J to 


+ C3 j {rj-s) “ sup [1 + 

Jto to<T<S 


)lla]ds 




Hence, we have 


+ + < C1 + C2 r snp[l + \\uiT)\U + \\u'{T)Uds 

Jto to<T<S 

I iv-s)-°‘ sup [l + ||u(r)(|i + ||u'(r)yds. 


+ Cs 


Using Lemma 4.4.2, we have 


sjip [i + y(^)lli + IK(?7)y < a 

^o<r}<s 


Hence, we get the required result. This completes the proof of the theorem. 



Chapter 5 


Second Order Semi-linear 
Integro-difFerential Equations - II 


/ have tried to avoid long numerical 
computations, thereby following Riemann’s 
postulate that proofs should be given through 
idea and not voluminous computations. 

— David Hilbert 


5.1 Introduction 


In this chapter, we shall continue the study of a second order semi-hnear integro- 
differential equation. In chapter 4, using the semigroup theory and the contraction 
mapping principle, we have proved the existence of a unique local classical solution of 
the abstract second order semi-linear integro-differential equation 
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in a Banach space and discussed the continuation of this solution, the ma vi m al interval 
of existence and the global existence under additional growth conditions on the nonlinear 
maps / and g. In this chapter we establish the results for the approximation to the 
solution of this equation in a Hilbert space. 

Let n be a bounded domain in with sufficiently smooth boundary dCl. Con- 
sider the following initial boundary value problem for the strongly damped semi-linear 
integro-differential equation 

^{x,t) + (aL + bI) {x,t) + {cL + dI)u{x,t) = h{x,t,u{x,t),^{x,t)) 

+ k(t - s)g(x,s,u{x,s),—{x,s))ds, (x,t) G x (0,T), 

Jto 

3x1 

u{x,0) = Xo{x), —{x,Q)=Xi{x), xeQ, (5.1.1) 

with the homogeneous Dirichlet boundary conditions, where Lu = ^ (%(^)^) 

be a synametric second order strongly elliptic differential operator in f2, a > 0, 6, c, d 
are constants and h and g are smooth nonlinear functions and kisa. locally p-integrable 
function for 1 < p < oo. 

In the abstract form, we may write (5.1.1) as the initial value problem 

= fit,u{t),^{t)) + k{t-s)g{s,u{s),^{s))ds, t>to, 

w(0) = xo, ^(0) = rri, (5-1-2) 

in the separable Hilbert space H = L^{Q), where the linear operator A with the domain 
D{A) = H\n) n ffi(Q) is given by 


Au = aLu, u £ D{A) 
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and the operator B is such that D{B) = D{A) with Bu = ciA^-c^I for some constant 
Cl and C 2 . The function / is defined from [0, oo) x if x F into if given by 

f{t,u,v) = h{t,u,v) ~ bv. 

In this chapter again we are concerned with the convergence of the Faedo-Galerkin 
approximation of the solution to (5.1.2). Without loss of generality, we concentrate on 
the abstract second order semi-linear integro-differential equations of the form 

u"(t) + Au'{t) = f{t,u{t),u'{t)) + f k(t - s)g{s,u{s),u'{s))ds, 

Jo 

u{0) = xq, u'{0) = Xi, ( 5 . 1 . 3 ) 

in a separable Hilbert space (if, ||.|1, (., .)), since the other terms may be transferred on 
the right hand side and merged with / as the condition satisfied by the changed / still 
remains the same. We extend the technique used in Chapter 3 for a class of second 
order semi-linear evolution equations to integro-differential equation (5.1.3). 

The organization of this chapter is as follows. In the second section we state some 
preliminaries and the assumptions required for establishing results. In the third section 
we consider a pair of associated nonlinear integral equations to (5.1.3). The solutions to 
these integral equations are defined as mild solutions to (5.1.3). With the help of the pair 
of associated nonlinear integral equations and the projection operators, we consider a 
pair of approximate nonlinear integral equations and show the existence and uniqueness 
of a solution to this pair of approximate integral equations using the Banach contraction 
principle in this section. In the fourth section we estabUsh the convergence of solutions 
and convergence of the pair of approximate integral equations to the pair of associated 
integral equations, limit of the solutions to the approximate integral equations being 
the solution of the pair of associated integral equations. Further, we show in this section 
that the solution can be extended to the maximal interval of existence and it is unique. 
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Finally, in the fifth section we consider the Faedo-Galerkin approximations of solutions 
and prove some results concerning the convergence of such approximations. 


5.2 Preliminaries and Assumptions 

In this section, we consider some preliminaries and assumptions essential for our pur- 
pose. The operator A assume to satisfy the following assumptions; 

(HI) A is a closed, positive definite, self-adjoint, linear operator from the domain 
I?(A) C H of A into H such that D[A) is dense in H, A has the pure point spectrum 

0 Aq ^ Aj ^ A 2 ^ 

and a corresponding complete orthonormal system of eigen functions {ui}, i.e., 

Aui = XiUi and (ui,Uj) = <5^, 
where 5ij — 1 if 2 = j and zero otherwise. 

If (HI) is satisfied, then —A generates an analytic semigroup in H which we denote 
by e~‘^, t>0. 

Now, we mention some notions and preliminaries. It is well known that there exists 
constant M > 1 and a real number oj such that 

|le-‘"^ll < Me^*, t>0. 

Since — A generates the analytic semigroup t > 0, we may add cl to —A for some 
constant c, if necessary, and in what follows we may assume without loss of generahty 
that l|e~*^|| is uniformly bounded by M, i.e., < M and 0 e p{A). In this case, 

it is possible to define the fractional power A°‘ for 0 < a < 1 as closed linear operator 
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with domain D{A°‘) C H (cf. Pazy [77], pp. 69-75 and p. 195). Furthermore, 
is dense in H and the expression 

\Ma = P“a:||, 

defines a norm on D{A°‘). Henceforth, we represent by X^, the space D{A°‘) endowed 
with the norm ||.||q;. In view of the facts mentioned above, we have the .following result 
for an analytic semigroup t > 0 (cf. Pazy [77] pp. 195-196). 

Lemma 5.2.1 Suppose that -A is the infinitesimal generator of an analytic semigroup 
e-‘^, t > 0 with ||e-*^|| < M for t > 0 and 0 E p{A). Then we have the following 
properties. 

(i) Xq, for 0 < Q; < 1 is a Banach space. 

(ii) For 0 < P < a, the embedding X^ Xg is continuous. 

( Hi) A°^ commutes with and there exists a constant Ca > 0 depending on a 
such that 

||A“e-*"^|| < CJ-^, t>0. 

The kernel k and the nonlinear maps / and g are assumed to satisfy the following 
assumptions. 

(H2) k E LFioJfi, oo) for some 1 < p < oo be a bounded function. 

(H3) The map / is defined from [0, oo) xXi xXq, into H and there exists a nondecreasing 
function Fr from [0, oo) into [0, oo) depending on i? > 0 such that 


||/(t,w,'u)|| < Fi?(i), 

\\fit,Ui,Vi) - f{t,U2,V2)\\ < FH(i){||^il -«2||l + ||U1 



5.3 Approximate Integral Equations 


126 


for all it,u,v) and {t,U 2 ,V 2 ) in [0,oo) x Br{Xi x (xo,a:i)), where 


BR{XiXXc,,{xo,Xi)) = {{x,y)eXixXa : l|x - Xolli + |ly - < i?}. 


(H4) The map g is defined from [0, cxd) x x X^ into H and there exists a nonnegative 
bounded function Gr € oo) depending on i? > 0, where l<a<oo 1 + 1 = 1 

such that 

||^(t,u,u)l| < GR[t), 

\\9{t,ui,vi) - g{t,U2,V2)\\ < C?n(i){|ltii-U2|li-l-||ui-U2iU}, 
for a.e. t e [0,oo) and all (u,u), {ui,Vi) and (^2,1)3) in Br{Xi x Xa,{xo,Xi)). 

5.3 Approximate Integral Equations 


we continue to use the notions and notations of the earlier section. The existence of a 
solution to (5.1.3) is closely associated with the following pair of integral equations 


u{t) 


v{t) 


xo + (e-‘^ - /)(-A)-Ui + - /)(-A)-^[/(s,n(s),n(s)) 

Jo 

+ / k{s - T)g{T,u{T),v{T))dT]ds, 

Jo 


(5.3.1) 




/* e (i «)^[/(s,u(s),u(s)) + f k{s-T)g{r,u{T),v{T))dT]ds, 

Jo Jo 


for t > 0. In this section, we shall consider a pair of approximate integral equations 
associated to (5.3.1) and establish the existence and uniqueness of a solution to the pair 
of approximate integral equations. By a solution {u,v) to (5.3.1) on [0,T], 0 < T < 00 , 
we mean a function {u, v) € Xi(T) x XaiT) for some 0 < a < 1 satisfying (5.3.1) where 
■^i(T') X Xa{T) is the Banach space C([0,r], D{A) x D{A°‘)) of all continuous functions 
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from [0, T] into D{A) x D{A°‘) endowed with the supremum norm 

UK v)\\xi(T)xXc{T) = |klUi(T) + ll?^|Uc(r), 


where 


llwllxi(r) = sup llAu(t)|| = sup \\u{t)\\i 

0<t<T 0<t<T 


and 


ll«^IU<,(T) = sup ||A“i;(t)|| = sup ||u(t)|U. 
o<t<r o<t<r 

By a solution {u,v) to (5.3.1) on [0,T)) 0 < T < oo, we mean a function {u,v) e 
Xi{T) X Xa(T) for some 0 < o; < 1, satisfying (5.3.1) in [0, T] for every 0<T <f. 


Let 0 < To < oo be arbitrary fixed and 

L(R) = (1 + i2)[F^(To) + llA:||Lp(o,ro)||<2;^||x,9(o,ro)], 

where 

R = R, + i? 2 , Ri = ^R^ + \\xo\\l and R 2 = y/W+^a- 
Let 0 < T < To be such that 

sup {IKe-''’- 7)1.11 + ||(e-“- /MV, 11} < j 

0<t<T 3 

and 

T < minjr„,|(M + l)-M(JJ)-M|cy(l-a)7{K)-'lA|. 

Let Hn denote the finite dimensional subspace of H spanned by {uq, «i , . . . , Un} and for 
n = 1, 2, . . .; let F” : i? — Hn be the corresponding projection operator. For each n, 
we define 

fn, gn:[0,T]xXi{T)xX^{T)-^H, 


such that 


fn{t,u,v,) = fit, P'^u(t),P'^v{t)), 

Quit, u, V, ) = git, P‘^u{t),P'^v{t)). 
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We set xo{t) = xo and f'i(t) = a:i for t € [0, T]. Let Wr = Br{Xi{T) x X^{T), (fo, xi)), 
where 


Br{Xi{T) X Xc.{T),(xo,x\)) 


= {(yi. y 2 ) e Xi(T) X X^(T) ; \\yi - fo|Ui(r) + ||y 2 - ^'ilUc(r) < H}. 


Define a map <S„ on Wr such that Sn{u, v) := {u,v) with 


u{t) = xo-{e ^^-I)A ^xi- / (e - I)A-\f^{s,u,v) 

JQ 

+ / k{s -T)gn{T,u,v)dT\ds, 

Jo 


( 5 . 3 . 2 ) 


0(t) = 


e"‘^xi + 




^^[/n(s, M, u) + / fc(s - T)i?„(r, ti, v)dT]ds. 


Propostion 5.3.1 Let (H1)-(H4) hold. Then there exists unique {un,Vn) G Wr such 
that Sn{un, Vn) = {un, u„) for each n = 1, 2, ..., i.e., (u„, Vn) satisfies the pair of approx- 
imate integral equations 


Un{t) = Xo-(e ^^-I)A ^Xi - [ (e ^[/n(s,MTi, ^^n) 

Jo 

+ / /i:(s -T)5f„(T,ii„,Pn)d'r]ds, 

yp 


( 5 . 3 . 3 ) 


Un(t) = e *^Xi+ / e ®^'^[/n(5,'Un,'yn)+ / A:(s - r) 5 n(T, ■?;n)dr]ds. 

Jo Jo 


Prop/. First, we claim that Sn '■ Wr — >• Wr. For this we need to show first that the 
map t {Sn{u,v))(t) is continuous from [0, T] into Xi x X^ with respect to norm 
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|.||i + ||.||q. For ii, ^2 G [ 0 ,r] with ti < we have 


Il'u(t2) - «(^l)||l + \\v{t2) - ^(tOllQ; 


-{t2-s)A _ g-(«l-. 


+ / ||e-(‘2-^)A _ /|| 

'ti 

+ r „e 

Jo 

rt2 

+ / 

Iti 

+ r II (e 

Jo 


11^ ^ -^11 fnis,U,v) + j k{s - T)gn{T,U,v)dT 

1 ^ 

fn{s,U,v)+ k(s-T)gn{T,U,v)dT 

J 0 

f ||e fnis,U,v) + [ k{s -T)gn{T,U,v)dT ds 

"'ii Jo 

^ )^°'ll fnis,U,v)+ f k{s — T)gn{T,U,v)dT 
Jo 


{t2-s)A _ 


Now, we estimate 


/n(s, u, t;) + / k{s - T)gn{r, u, v)dT 
Jo 

< ll/n(5,u,z;)|| + f \k(s - T)\\\gn{T,U,v)\\dT 

Jo 

^ -^^(^0)+ [ \k{s-T)\Gj^{T)dT 
Jo 

^ -^r(^o) + ||A:|lLP(o,ro)||<?fllU9(0,to) < L{R)- 


Also, we have 


r*2 


/ \\e-(^2-s)A_^ 

Jti 


fn{s,u,v)-\- / k{s-r)gn{T,u,v)dr 
Jo 

< {M+l)L{R)(t2-k), 


ds 


ds. 

(5.3.4) 


(5.3.5) 


Jti 


fn{s,U,v)+ f k{s-T)gn{T,U,v)d- 

Jo 

< L{R)Ca f\h-s)-^ 

Ju 


ds 

ds 


= L{R)Co, 


{h — k) 

1 — a 


l-a 


(5.3.6) 
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Part (d) of Theorem 1.3.19 implies that for 0 < < 1 and x e D{A^)^ we have 




(5.3.7) 


If 0 < -i? < 1 is such that 0 < a + < 1, then e D{A^) for any y e 
Therefore, for t,s E [0, T], we have 




(5.3.8) 


From (5.3.8), we get 

Hence 


f 

Jo 


:,-(t2-s)A 


,-(ti-s)A^^a|| 


fn{s, u,v)+ / k{s - T)gn{T, u, v)dr 
Jo 


ds 


< C!,Ca+^L{R){t2 - h)^ f\h - s)-(“+’')ds 

Jo 






(5.3.9) 


Also, from (5.3.7), we have 


li(e-*^-/)e-"^a:|| < C'^f\\e-^^xU = C'^f\\A^e-^^x\ 


< CLC^fs-^ X . 
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Therefore 

||g-(t2-3)A ^ ||(e-(*2-‘iM_ j)g-{ti-sM|j 

< C'^C^{t2-hY{tr-s)-^. 


Hence 


|g-(t 2 -s)A _ g-(tl-f 


fn{s,U,v)+ f k{s-T)gn{T,U,v)dT 

Jo 

< C!gC^L{R){t2-hY Tih-sy^ds 

Jo 

< c;c,i(fi)^(t2-(i)« (5.3.10) 


From inequalities (5.3.4), (5.3.5), (5.3.6), (5.3.9) and (5.3.10), it follows that {Sn{u,v)){t) 
is continuous from [0, T] into D(A) x D{A°‘) with respect to norm ||.||i + ||.||q,. Now, we 
show that Sn{u,v) e Wr i.e., (ft, t)) € Wr. Consider 


P(^)-2:ol|l + |lf)(t) - Xilla 

< ll(e-‘^-/)xill + ll(e-‘^-J)A“:rill 

I 


+ / - /|| 
lo 

rt 


Jo 
R 


fn{s,U,v)+ [ k{s - T)gn{T,U,v)dT 

Jo 

fn{s,U,v)+ f k{s-T)gn{T,U,v)d 
Jo 


ds 

ds 


< j + {M + l)L{R)T + CaL{R) J {t-s)-^ds 

< ^ + {M + l)L{R)T + CaL{R)^^ < R. 


Taking supremum over [0,T], we get 


- £o|Ui(t) + 11^ - ^llUa(T) ^ 
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which implies that S„(«, v) £ W^. So, S„ maps Wn into W^. Now, it remains to show 
that is contraction on VFh- For (ui, ^;i), ^ 

||wi(t) — ^2(0111 + ll^l(i) - ^ 2 (i)|la 

^ lo ~ ~ 9n{r,U2,V2)\\dT}ds 

+ ^‘||e-<‘-*'">l“||{||/„(s,nr,t,r)-/„(s,«s,ns)|| 

+ / -9»(T,ti2,t)2)||dr}ds. 

(5.3.11) 


FVom (H3), we have 


||/n(-s, Wi, r^i) — fn{s, 112, 1)2)11 

= \\fis,P^Ui(s),P^Vi(s)) ~ fis,P^U2{s),P%2{s))\\ 

^ P'R(TQ)i\\P^Ui{s) - P”«2(s)||l + ||F"Ul(s) - P”?;2(5)||a) 

- FR(Fo)(||^il(s) - U2(5 )||i + llt'l(s) - ■y2(5)||a) 

^ Fr(Fo)(||ui - 'U2|Ui(r) + ||t^i - 'f^2|Uc(r))- 

Similarly, from (H4), we get 

||£^n(<S, Wi,ni) - ^n(s,W2,U2)|| 

- ^a(^)(II“i “ 'f^2|Ui(T) + ||ui - n2|U„(r))- 
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Therefore 


||/„(s, til, Pi) - fn{s, U2,V2) II + / \k{s - r)| ||g„(r, Ui,Pi) - g^{T, U2, V2)\\dT 

J 0 

< [Fr(7'o) + il^lUp(0,To)||G^||L«(0,ro)](||'Ul - W2|Ui(r) + ||Pi - 'y2|U„(T)) 


< 


L(H) 

R 


(l|txi - U2|Ui(r) + ||Pi - U2lUc.(r))- 


(5.3.12) 


Hence, using the estimate (5.3.12) in (5.3.11), we get 


l|ui(£) -U2(t)||i + ll^i(£) -P2(i)lia 

< |(M + l)L(i?)r + CaL{R){To) J\t - s)-“ds| 

X (||ui — ti2lUi(r) + \\vi — P2|U„(r)) 

< i|(M + l)L(i?)T + C,L(Ji;)^| 

X (||ui — U2||xi(r) + ll^^i - 'y2|Ua(r)) 

2 

< 3(11^1 “ '“2 |Ui(T) + ll^^l - 'y2lU„(T))' 

Taking supremum over [0, T], we get 

2 

ll«l - «2|Ua(T) + 11^1 - P2|Ua(T) < -^{Wui “ U2\\xi(T) + 1^1 “ 'y2lUc(T))- 

Thus, Sn is a strict contraction on Wr. Hence, there exists a unique {un,Vn) € Wr such 
that Sn(un,Vn) = (wn,t^n)- Clearly, {un,Vn) satisfies (5.3.3). This completes the proof 
of the Proposition. 

Propostion 5.3.2 Let (H1)-(E4) hold. If {xo,xi) e D{A)xDiA), theniun{t),Vn(t)) € 
D{A) X D{A'^) for all t € [0,T] where 0 < < 1. 


Proof Prom Proposition 5.3.2 we have the existence of a unique (u nj ^n) € Bk{X,(T) X 
^a{T), {xq,xi) satisfying (5.3.3). By Theorem 1.3.9, we have for x E H, f^e xds E 
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D(A) and if a: 6 D{A) then e D(A). The result of the proposition follows from 
these facts and the fact D(A) C D{A^) for 0 < < l. 

Propostion 5.3.3 Let (H1)-(H4) hold. Then {or x, £ D[A) there exists constant Vo 
independent of n, such that 

l|■^n(^)||13 < Vq, 0<1?<1, 0<t<T. 

Proof Applying on both the sides of second integral equation of (5.3.3) and taking 
norm, for 0 < t < T, we have 

ll«»Wlle < P‘’e-“xr||+ 

Jo 

+ \K^ -T)\\\gniT,Un,Vn\\dT]ds. 

J 0 

Since Xi G D{A) implies that xy G D(A^) for 0 < < 1, we get 

ikn(i)||t? < Ar||a;i||,j + C',sL(i?)- < Vq. 

1 — 

This completes the proof of the proposition. 

5.4 Convergence of Approximate Solutions 

In this section, we establish the convergence of the solution G Xi{T) x Xoc{T) 

of the pair of approximate integral equations 


Unit) = Xo- {e-^^ - I)A-^Xy - f\t 
s 

+ / kis - T)gn{T,Un,Vn)dT]ds, 
Jo 


_ n 4-1 




(5.4.1) 


Unit) = 


Xi+ [ e ^^^[fnis,Un,V n)+ [ kis - T)gniT,Un,Vn)dT]ds, 

Jo Jo 
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then there exists a constant K such that 

f{t) + f f{r)dT < AK for 0<t<T. 
Jo 

Furthermore, 

f{t) < AK for 0<t<T. 


Proof Integrating the inequality (5.4.2) with respect to t for 0 < t < T and change the 
order of integration in the double integral on the right hand side, we have 


[ f{s)ds 
Jo 

< 

AT + B [ 
Jo 


< 

AT + B [ 
Jo 


< 

AT+B f 
Jo 

Estimating (t — 

sy- 

■“ by f( 

f f{s)ds 
Jo 

< 

BT^ 

AT +- — 


t nr) 


t / rt 


^dr] 


/( 5 )+ f f{T)di 

Jo 


T ds 


(5,4.3) 


Adding the inequalities (5.4.2) and (5.4.3), we have 


fit) + f{r)dT < A(r + 1) + 15 ^ ^/(s) + ^ /(r)dr^ ds 

< A{T + 1) + .B(r“ + 1) ^ Jo 


where 


B = max 


r BT^ 

f'— 


l-a 

a 


The required result follows from Lemma 1.3.27. 
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Propostion 5.4.2 Let (HI), (H2), (H3) and (H4') hold. //(xo,Xi) e D{A) x D{A), 
then 


lim sup {\\un{t) 

{n>m, 0<4<r} 


Mm (t) 111 + ||u„(t) - Um(t)||a} = 0. 


Proof. For n>m, we have 


||/n(t) Myii Mjj) 

< \\fn{t,Un,Vn) - fn{t,Um,Vm)\\ + \\fn{t,Um,Vm) - fm{t,Um,Vm)\\ 

< Wfit, P^Unit),P%n{t)) - fit, P^Umit), P^Vmitm 

+ \\fit,P^Umit),P^Vm{t)) - f{t,P^Umit),P^Vmit))\\ 

< F^{To)[\\P-Un{t) - P-uUt)\\i + WP^Vnit) - P"u^(t)|U 
+ ||(P" - P”^)u^(t)||i + ||(P” - P”^)t;„(t)|U. 


Also, 


ii{p" 


< 


||A“(F* - F")»m(<)ll = - P”).4’’»„(t)|| 

^|1(P”-p")aV(()|| < 


For convenience, we denote 

Cm, nit) = llUnW “ MmWIIl + lbn(t) -''ym(i)|U- 


Thus, we have 


1 1 fn (t 5 '^n ) '^n ) 


fm(t^ '^m) I1 

< Fk(To) + ll(P" - P’”)««.(*)lli + • P-4-4) 
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Similarly, we have 


||5n(^) ^ni ^n) 

< Gs(t) ("{„,„(() + |l(/>“ - F“K(()||, + . (5.4.5) 

From (5.4.5) and Proposition 5.3.3, for n > m, we have 



I ^n) ’^n) 9m(j'i‘^m.^'^m)\\dT 


< P!Up(0,To)l|G^IU»(0,ro) ( sup ||(P” - P”^)Um{T)\\i + 

Vo<T<s A"-“y 

+ / l^(s-r)lG^(r)^^,„(7-)dT. (5.4.6) 

Jo 


Hence, from (5.4.4) and (5.4.6), we get 


|l/n(s, Wji, Un) /m(s, ^02, Vrn) II 

+ [ |/c(s-r)|llp„( T^Un)Vfi) ^rri) '^m) 

Jo 

< L(R) ( sup ||(P“ - P'")«™(t)||, + ^) 

Vo<T<s Ki°‘J 

+ (p«(Io)W(«) + ^ - ^IGsWfm.sW*) 

< L(R) ( sup ||(F*-P">)«..(r)lli + ^) 

\0<r<5 / 

+ N Un{r)dr'^ , - (5-4-7) 


where N = max{F^(To), BkBa^} and Bk and Bq^ be bounds for k and G^, respectively 
as these are bounded functions of t. 
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Now, from the pair of integral equations (5.3.3), we have 

J 0 


+ 


[ 1^(5 -'r)||kn(r,ii„,'y„) - 

Jo 


^rri) )l|dr](is 


+ X - fra{s,Urn,Vrn)\\ 

^ lo ~ '^)\\\9n{r,Un,Vn) - gm{r,Um,Vra)\\dT]ds. 


Using the estimate (5.4.7), we get 


U.W < ((M+l) + ^^^(^gpJ|(P«-P".)«„M||, + ^)i3 

+ N f ({M+l) + ,^)(u,„ (s) + f ^m,n{'T)d'l^ ds 


(t - sY 


where 


< C{R, T)B^n (5”>."(s) + im,n(T)dT\ ds, 


B„„= sup ||(P”-P'")u4r)|h + ^, 


0<t<T 


CiR, T) = L{R) ({M + 1)T + 9^^ 


\' ' L-a 

Ni = N{T^ + 1) max{(M + 1), Ca}. 

Hence, from Lemma 5.4.1 there exists a constant K such that 


^m,n(t)+ [ im,n{r)dT < KC{R,T)Bmn, 

Jo 


which implies that 

Unit) < KC{R,T)B^n. 

Taking supremum over [0, T] and the limit as m oo on both sides, we get the required 
result since Bjnn — > 0 as m oo provided ||(P" — P™)'Um(^)lii ^ 0 as m ^ oo for 
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Now, we prove that for 0 < t < T, \\{P- - 0 as m -> oo. We can 

easily check that for every x e H and r] <0, 

l|.4’'{P" - P”)I|| < A:!,||(P“ - P”)x|| < A^||x||. (5.4.8) 

From the equation (5.4.X), we get 

Jo 

+ / \Ks -r)\\\{P'^ - P'^)gn{T,Un,Vn)\\dT]ds. 

Jo 

(5.4.9) 

We note that in (5.1.3), f{t,u,v) = f'^{t,u,v) - Bu where represents original / 
appears in (5.1.2). Also, we can check 11 jBu|| < C||Au|| for u G D{A) and some constant 
C. Therefore, we have 

WiP^ - Pnfm{s,Um,Vm)\\ < \\{P^ - P"^)f^{s,Um,Vm)\\ 

+ \\A{P^ - P^)Um{s)\\. (5.4.10) 

Since A^f°{s ,Ujn,Vm) G H, hence from (5.4.8), we have 

< A^r(T«)- 

Xm 

Similarly, as A^gm{s, Um, Vm) G H, we have 
ll(P”-P’»)s„(T,«„,t,„)|l < 

Am 
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Using these inequalities in (5.4.9), we get 

||(P«-P'")„„(t)||l < ||(P’‘-P”)Al„|| + (M+l){||(P«-P™)l;|| 

+ PIU’’(0,To)ll^filU«(0,To)) 

+ f‘ WiP" - P'")u„(s)\Uds}. 

Jo 

Gronwall’s inequality implies that 

l|(P”-P”)«„(t)lli < [||{P”-P”Mio|| + (M+l){l|(P"-P“)n|| 

^ ^r(Pfi(?o)ll^ilLp(o,r(,)l|<5,;l!Lfl[o.ro)}}]e^^^^*^ 

which tend to zero as m oo for 0 < t < T. This completes the proof of the 
proposition. 

For the convergence of the solution u„(t)) of the pair of approximate integral 
equations (5.4.1), we have the following result. 

Theorem 5.4.3 Let (HI), (H2), (H3) and (H4') hold and let {xq,xi) G D{A) xD{A). 
Then there exists function {u,v) G Xi{T) x Xq{T) such that {un,Vn) — > {u,v) as 
n oo in Xi{T) x Xa{T) and {u,v) satisfies (5.3.1) on [0,T]. Furthermore, {u,v) 
can be extended to the maximal interval of existence [0, tTnai); 0 < < oo satisfying 

(5.3.1) and {u,v) is a unique solution to (5.3.1) on [0, tmox)- 

Proof Proposition (5.4.2) implies that there exists {u,v) G Xi{T) x Xa{T) such that 
{un,Vn) converges to (u, u) in Xi{T) x Xo,{T). Since {Un,Vn) £ Wr for each n {u,v) is 
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also in Wr. Further, we have 

= ll/(<, P’Mt), - f{i. «(t),u(«))|| 

< ll/(t, ■P'tinW. F“o„(t)) - f{t, P”u{«), F'v{t))\\ 

+ |l/{«.P"u(t).P“»(i)) - /(t,«((),u(«))|l 

< PR{To)l\\P"Mt) - F‘u(t)|h + l|PX(i) - P«„(l)||„ 
+ ||(P“-/H()|l, + ||(P"-;)u(t)||J 

^ ^R(ro)0i«n{f) - ll(*)||l + ||w„(t) - k(()|U 
+ ||(P“-i)«{i)||, + ll(P"-;M«)iy. 

Taking supremum over [0, T], we get 

sup \\fn{t,Un,Vn) “ f{t,u{t),v{t))\\ 

0<t<T 

^ •f’R(^o)[|lWn - ullxi(T) + ll^^n “ ^lUaCT) 

+ |l(P"-/HU,(T) + ||(P”-/>lk(r)]. 

The right hand side of (5.4.11) tends to zero as n tend to infinity. Similarly 

\\gn{t,Un,Vn) “ g{t, u{t) ,v{t))\\ 

= \\g{t, P^'Unit), P"n„(t)) - g{t, u{tl v{t)) II 

< \\g{t, P^Un{t), PX(t)) - 9{i, P”^;(t))ll 

+ \\g{t,PMt),PMt))-9itMt)At))\\ 

< G^{t)[\\P^Un(t) - P"«(t)||i + ||PX(t) - PMt)\U 
+ ii(p"-j)u(t)iii+ii(p"-x(t)y 

< G^(^)[yn(t) - u{t)\\i + Il^;n(i) - v{t)\\a 
+ ||(P"' — /)u(t)||i + |1(.P” "" .^)^(^)l|a]- 


(5.4.11) 


(5.4.12) 
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5.4 Conv ergence of Approximate Solution s 
Prom (5.4.12), we have 

sup / IKi- s)\\\gn{s,Un,Vn)~g{s,u{t)Mt))\\ds 
< L(i?)[||Un - «|Ui(r) + l|Un - ylU^(r) 

+ - I)u\\xi{T) + |1(P” - /)uiU„(r)] -^0 as n-^oo. (5.4.13) 


Hence, using (5.4.11), (5.4.13) and the bounded convergence theorem in (5.4.1), we get 


u{t) 


v{t) 


xq — (e — I){A) ^xi — [ {e — I){A) ^[/(s,'u(s),t;(s)) 

Jo 

+ k{s -r)g{T,u{T),v{T))dT]ds, 

Jo 


e'^^xi + 



e ")^[/(s,n(s),u(s)) + 



T)5(r,'u(r),v(r))dr]ds. 


Thus, for ( 2 ;o,a:i) € D{A) x D{A), there exists a unique pair of functions (u,u) € 
Xi{T) X Xo,(T) such that (u„,Un) — > (u, u) as n oo in Xi{T) x Xq(T) and {u,v) 
satisfies (5.3.1) on [0,T]. 

If {u, v) satisfy (5.3.1) on [0, T], then we show that (u, v) can be extended further. 
Since 0 < To < oo was arbitraxy, we assume that 0 < T < Tq. We consider the equation 


U"{t) + AU'{t) 


F{t,U{t),U\t)) + 



U{ 0 )=u{T), U’(0) = u\n 


s)G{s,U{s),U'{s))ds, t>Q, 


which can also be written as system of equations 


U\t) = V{t), U{0)=u{T), 

V\t) + AV{t) = F{t, U{t), V{t)) + f k{t- s)G{s, U{s),V{s))ds, 7(0) = v{T). 
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where F, G : [0,To ~ T] x X. y Y u- 

are defined by 


F(t,x,x) _ f{t + T,x,x) + h{t), 
G{t,x,x) = g{t + T,x,x), 


where 


for ft *) s ( 0 , 21 - rj X V. X X.. F and G satisfy {H 3 ) and 

replaced by r„ - r Hence tl, ■ ^ ' 3 ) and (H 4 ), respectively, for r„ 

y 2 „ r. Hence, there exists a (ft V) e C(( 0 , r,), D(A) x i 5 M.il f„ 

0 < T, < 71 _ T' ooM f. • , ''1 ’ IJ) X )) for some 

- 0 ^ satisfying the integral equation 


U{t) = 


V{t) = 


- <r)^+(e 

+ - !){~AriF(,,u(s),V(s)) 

+ 1 k{,-T)G(T,U(s).V(,))dT]ds, 

+ 1 ft(s), V(s)) 

/ *^(«-r)G(r,ft(s),l/(s))<,^j*^ 


with 0 < t < 71. Now, we define 

5(i)) = / 0 < t < r, 

1 lu{t~T),V{t-T)), T<t<Tj + T. 

, for t s [ 0 , 71 + rj, (a, satisfies the integral equations 

+ (a- - 7)(_^)-i,. , 

- 'r)y(v, u{T),v(T))dT]ds, 


Jo 


( 5 . 4 . 14 ) 

)^(’',«(T),5(r))dr]ds. 


5(0 = 
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To see this, we need to verify ( 5 . 4 . 14 ) only on fr T ri p 

> + Forte [T,Ti+T], 
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u(t) = U{t-T) 

+ ^ ~ r)a{T,V{T),V(T))dr]ds 

Jo l/(5,«(s),t>(s)) 

+ k(s - T)g(T, u(t), v{T))dT]ds + (e-(^-T)A _ 

/*T 


+ e 
Jo 

rt-T 


+ 

+ 


k(s - T)g(T,u(T),v(r))dT]ds) 

I - /)(-.4)-'[f (s, c/(s)_ i/(s)) 

I '=(" - ^)G(rV(r), y(r))rfrl*. 


Putting s + T = 7 ] and r + T = ^, 


we get 


u(t) 


= .0 + (e- - ^ __ 

+ ^ Ks~T)g(T,u(T),v(T))dT]ds 

+ £ __ 

p-T 

+ J^ Kn-T~T)G(T,U(r),V{T))dT]ds 

= Zo + (e-« - 

ns ^ 

+ y ^(^-^)^(^,^(r),t;(7-))dr]ds 
+ y^ A:(i; - - T, 17 (e - T), F(e - T))di]dg. 



Since 


F(s - T, U(s - T). V(s - D) = /(., U(s - n V(s - T)) + - T) 


and 


we have 


G(s - T. U(s - T), V(s~ D) = g(s, U{s - T), Vi 


s-r)), 


u{t) = xo + (e 


—tA 


+ 1 - n(-A)-'f(s, uis), v{,))ds 

^ i I - rUr, uir),vir))drds 

^ ir Jo * * ~0(-A) 'k{s-T)g(T,u(r),v(T))dTis 

^LL - t )«( t , fi(r), S(r))drds 

~ s(,)) 

+ / k{s- T)g{T, u{t), v{T))dT]ds. 

J 0 


Now, for n, we have 


H-T 


v{t) = V(t-T) 

= + f e-»-^-’->^[F(s,U(s),Vis)) 

J 0 

+ f k{^-T)G{T,U(T),V{T))dT]ds 

J 0 

= + / e-(^-^)^[/(5,-n(5),n(s)) 

J 0 

+ / k{s-T)g{T,u{T),v{T))dT\ds) 

J 0 


+ 


Jo 


.-(t-T- 


")^[F(s, U{s),Vis)) + f k{s- r)G{T, U{t), V{T))dT]ds. 
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Again, putting T + s 7] and T" + r ^ and using the definitions of F and G, we get 


v{t) 


e + I e ® 


/ e ^ *^"^[/(S)n(s),u(s)) -f- f k{s — T)g{T,u{T),v{T))dT]ds 
-'0 , Jq 

+ f e-^^-^^^[F{rj-T,U{r,-T),V{n-T)) 

Jt 

rv-T 

+ / k{r]-T-r)G{r,U{T),V{T))dT]ds 

Jo 

e f e u(s),u(s)) + / k(s - T)g{T,u{T),v{r))dT]ds 

■JO Jo 

+ /, V(7, - T),V{r, - T)) + h{n - T) 

+ Kri - OGK - T, U(i - T), V(i - T))di]dri 


= / e-»-‘^^f{s,u{s),v(s))ds 

/o 


= e 


+ f / e ^^'^k{s — T)g{T,u{T),v{T))dTds 
Jo Jo 

+ [ [ e~^*'~^^'^k{s — T)g{r,u{T),v{T))dTds 
Jt Jo 

■F f f e~^^~^^'^k{s — T)g{T,u(T),v{T))dTds 

Jt Jt 

j:i+ f e~^'^~"'>^[f{s,u{s),v{s))+ f k{s - T)g{T,u{T),v{T))dT]ds 
Jo Jo 


Thus, we see (u(t),v{t)) satisfy (5.3.1) on [0, Ti +r]. Hence, we may extend {u{t), v{t)) 
to the maximal interval [0,tmax) satisfying (5.3.1)on [0,tmax) with 0 < < oo. 

Now, we show the uniqueness. Let and (^ 2 ,^ 2 ) be two solutions to (5.3.1) 

on some interval [0, T], where T be any number such that 0 < T < tmax- 

For convenience, we denote 

i{t) = \\ui{t) - ■U2(i)||l + Ibl(^) - t'2(t)||a- 
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Then, for any 0 < f < T, we have 

(5.4.15) 

Since 

ll»(s.«i{5),0i(s))-5(s,u,{s),„j(s))|j < Gs(s)f(s), 

we have 


Ilf(^’‘^i(^),'^i(s))-f(s,u2(s),v2(s))ll 

+ ^(^-^)lb(T«i(T),yi(r))-^(r,n2(r),n2(r))j|dr 

^ ^R{tmax)({s) + |A;(5 - r)|G^(r),f (r)dr 

- + 1 (Mdr'j . (5.4.15) 


where 


^' = ^ax{F^{tma.),BkBG^}. 

Using (5.4.16) in (5.4.15), we get 


m < N' f ({M + l) + 


a 


< ^2 


{t-s) 

lo Jo 
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where 

N 2 = N'{T°‘ + 1) max{(M + 1), Co}- 
Hence, from Lemma 5.4.1, we get 

CW-h / ^(r)dr = 0, 

Jo 

which implies that 

^(t) = 0 for 0<t<T. 

Prom the facts that 

l|ui(t) -tiaWII < ^llui(t)-U2(t)l|i, 

Ao 

l|Ul(f) -U2(t)|| < ■A||uj(t)_U2(t)||o, 

^0 

it follows that {ui,Vi) — (u 2 ,n 2 ) on [0,T]. Since 0 <T < tmax was arbitrary, we have 
(ui,Ui) = (^ 2 ,^ 2 ) on [0, tmax)- This completes the proof of the theorem. 

5.5 Faedo-Galerkin Approximations 

For any 0 < T < tmax, we have a unique pair (u,v) € Xi{T) x X^iT) satisfying the 
integral equations 

u(t) = xq- (e-*^ - I){A)-^Xi - f - /)(A)-^[/(s, u{s), v{s)) 

Jo 

+ [ A:(s-r)^(r,P”un(r),P"(T))dr]ds, 

Jo 

v(t) = e~^^xi+ f e“(‘“"^^[/(s,u(s),w(s)) 

Jo 

+ r k{s - r)g(r, P”Un(r), P'^{r))dr]ds. 

Jo 


(5.5.1) 
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Also, we have a unique solution (un.Vn) E Xi{T) x XaiT) of the pair of approximate 
integral equations 


1X71 (t) 


Vn{t) 


s ‘/O 

+ / k{s - T)g{T,P^Un{T),P%^{T))dT]ds, 

Jo 

(5.5.2) 

e-*^xi+ f e-(‘-^)^[/(s,P"u„(s),PX(5)) 

Jo 

+ f k(s- T)g{T, P”u„(t), P"'u„(r))dr](is. 


If we project the equations (5.5.2) onto Hn, we get the Faedo-Galerkin approximations 
{un(t),Vn{t)) = (P”u„(t),P”'nn(t)) Satisfying 


Unit) = 


Vn{t) = e' 


P^xo - (e-‘^ - /)A-'P"a:i - - /)A-'P"[/(s, u„(s), t)„(s)) 

Jo 

+ / A:(s-r)i?(r,tin(r),f)„(r))dT](is, 

Jo 

(5.5.3) 

P"a;i+ f e'^*"^)^P"[/(s,u„(s),?)„(s)) 

Jo 

+ / fc(s -r) 5(r,«n(T),Un('r))dr]ds. 

Jo 


The solution (tt,u) of (5.5.1) and (un,Vn) of (5.5.3), have the representation 


u{t) = y~^ai(t)ui, Q:i(t) = (u(t),Ui), i = 0, 1, 


i=0 


(5.5.4) 


;(t) = 5^A(t)ui, A(t) = (u(t),«i), i = 0,l, 
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and 


= Y. cx 7 {t) = {Un{t), Ui), i = 0, 1, . . . , n; 

1=0 

(5.5.5) 

n 

^n(t) = Pt{t) = {Vn{t),Ui), i = 0,l,...,n. 

2=0 


Using (5.5.5) in (5.5.3), we obtain the following systems of the first order integro- 
differential equations, 


da'^d) /■* 

-^ + A,<(f) = A* + *+ / 

+ /' k(s - r)Gr(r,a;(r) Q^Cr), A"{r), ffi{T))dr)ds, 

Jo 

(5.5.6) 


dt 




+ [‘ k{t- s)Gt{s,a’S(s),...,c^„{s),PS{a),...,/3:(s))ds, 

Jo 


with the initial conditions 

a?(0) = 0i, 

where 

<(t),)9y(t), ...,«(()) = 

and 4>i = {xo,Ui), V'i = for j = 1,2,.. 

The systems (5.5.6) determine the a" 


ff(0) = *, 


/(t,^o;7(t)ui,^/3r(t)ni),Ui j , 

2=0 ' 2=0 / 



,n. 

’s and Pi{tys. Now, we shall show the 
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convergence of (af , ) to (a, /?). It can be easily checked that 


A[u{t) — = A 


1=0 


1=0 


and 


= A° 




2=0 


= 'EKm)-mt)w 


i=0 


Thus, we have 


and 


M[u(t) -*(()] IP > ^A?(a,(t)-c”(t)p 


i=0 




i=0 


We have the following convergence Theorem. 


Theorem 5.5.1 Let (HI), (H2), (H^) and (HJ() hold. If {xo,Xi) e D{A) x D{A), 


then 


lim sup 

n-^oo o<t<T 


^ Xj(a,{t) - am? + E [ = 0. 


i=0 


i=0 


The assertion of Theorem 5.5.1 follows from the facts mentioned above and from the 
following proposition. 


Propostion 5.5.2 Let (HI), (H2), (H3) and (Hf!) hold and letT be any number such 
that 0 < T < tmax, then If (xq, xi) 6 D(A) x D{A), then 

lim sup {||A[«n(t) — Wm(0]ll + ~ = 0. 

0<t<T} 
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Proof. For n > m, we have 


\\A{Un{t) - U^(i))|| + \\A<^{v,(t) - i)„(i))|| 

= \\A{P^Un{t) - + ||A“(P"r;„(t) - P^Vm{t))\\ 

< \\AP'^{Unit) - ?/^(t))|| + p(pn _ 

- 'y,n(^))|| + ||A“(P^ - P"^)w^(t)|| 

^ ll'*^n(i) - Wm(i)||l + ||yn(0 “ Um(i)||Q: 

+II(p”-J="KWII. + -^II4V|| 

If (xo,Xi) G D{A) X D(A), then the result follows directly from Proposition 5.4.2. 



Chapter 6 


Quasi-linear Integro-differential 
Equations 


Geometry may sometimes appear to take 
the lead over analysis, but in fact precedes 
it only as a servant goes before his master 
to clear the path and light him on the way 

— James Joseph Sylvester 


6.1 Introduction 

In. this chapter, we consider the following quasi-linear implicit integro-differential equa- 
tion in a Banach space X, 

d/ll 

~{t) + A{t,u{t))u{t) = f(t,u{t),G{u){t)), 0<t<T, u{0) = uo, (6.1.1) 

where 0 < T < oo and C is a no nlin ear Volterra integral operator. In the first section, 
we prove the existence, uniqueness and continuous dependence on the initial data of a 
strong solution by using the method of semi-discretization in time which is also known 
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as the Rothe method or the method of lines. In the second section, we first establish the 
existence of a unique local mild solution and then prove the existence of a local classical 
solution using the semigroup theory and the contraction mapping theorem under some 
different appropriate assumptions. 

For applications of the theory of analytic semigroups to related quasi-linear evo- 
lution equations, we refer to Amann [1], Lunardy [69] and for applications of the fixed 
point theorems the reader may refer to Kartsatos [48, 49], Kartsatos and Parrott [50] 
and references cited therein. Kartsatos [49] and Kartsatos and Parrot [50] have estab- 
lished the existence and uniqueness of a solution and used the Galerkin method for the 
approximation of the solution to the quasi-linear functional problem 

u'{t) -(- A{t,ut)u{t) =0, t G [0,T), u(0) = <f>, 

in a Banach space X, where A(t, is m-accretive in v £ D (fixed) for every pair 
(t, € [0, T) X Co, J9 is a subset of X, Co is a closed subset of the space of all 

continuous functions 0 on (— oo, 0] into D with ||</>||co < ^ for some fixed r > 0 and 
Uf £ Co, given by Ut{s) = u{t + s), s £ (-oo, 0]. Murphy [73] has constructed a family 
of approximate solutions to the homogeneous quasi-linear evolution equation 

u'{t) + A(t,u{t))u{t) = 0, t£{0,Tl u{0) = (P. (6.1.2) 

He showed that the approximate solution converges to a limit solution and this limit 
solution becomes a unique solution to this homogeneous case of (6.1.1). Kartsatos [48] 
established a theorem concerning the existence of a unique strong solution to (6.1.2) 
under the assumption that A(t,u)v is Lipschitzian in t, u and m-accretive in v. In the 
second section, the condition on A(t,u)v is motivated by Kartsatos [48]. 

For the study of a particular case of (6.1.1), in which f{t,u,v) = f{t,u), that is. 
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the abstract quasi-linear evolution equation 

u'{t) + A{t,uit))u{t) = f{t,u{t)), te(0,T], u(0) = <^, (6.1.3) 

we refer to T. Kato [61, 58, 60, 62], S. Kato [54], Amann [1] and references cited therein. 
The crucial assumption in these papers is that the linear operator depending 

on t and the unknown u, has the dual property of being the negative generator of a 
Co-semigroup but not necessarily of an analytic semigroup and at the same time, a 
bounded operator on Y into X. Moreover, it is assumed in [61, 58, 60, 62] that there is 
an isomorphism S oiY onto X with the property that 

SA{t,w)S~^ = A(t,ui) + 

where B{t,w) is in the space B{X) of all bounded linear operators from X into X. In 
[54], it is assumed that there is a family {S(t,w)} of isomorphisms of Y onto X such 
that 

S{t,w)A{t,w)S(t.,w)~^ = A{t,w) -1- B{t,w), 
where B(t,w) € B{X) and S{t,w) satisfies the Lipschitz-like condition 

||S'(t,Wi) - S'(s,u;2)|ly,x < K\t-s\ + \\wi-W2\\x), 

where is a constant. T. Kato [61] has proved the existence, uniqueness and continuous 
dependence on the initial data of a solution of an abstract quasi-Unear and showed 
that these results are applicable to different kinds of quasi-linear equations such as 
symmetric hyperbolic systems of the first order, the wave equation, the Korteweg-de 
Vries equation, the Navier-Stokes and the Euler equation, the magnetohydrodynamics 
equation, the coupled Maxwell and Dirac equations, etc. The results in [61] are based 
on the theory of linear ‘hyperbolic’ equations which was developed in [57, 59]. S. Kato 
[54] has proved the existence of a strong solution to (6.1.3) under certain conditions on 
A{t,u) and f{t,u) for {t,u) E J xW similar to that of Crandal and Sougandis [28]. 
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Amann [1] has treated various cases of (6.1.3) in interpolation spaces using the theory 
of analytic semigroups. 

Using the Rothe method, Kacur [47] has proved the existence and uniqueness 
results for the nonlinear implicit integro-differential equation 

^{t) + Au{t) = G{t,K{u)it)), a.e. te(0,T), 
u = <f) on [-q,T], q>Q, 

in a Hilbert space H with the assumption that A \ Y — ^ V* is a coercive maximal 
monotone operator where U is a reflexive Banach space with the dual V* such that 
U n is dense in V and H and the function / and the Volterra operator K{u){t) 
satisfy certain Lipschitz-like conditions. 

The existence, uniqueness and continuous dependence on the initial data of a solu- 
tion to some abstract nonlinear implicit and exphcit integro-differential equations have 
been studied by Bahuguna and Raghavendra [8, 9] by using the Rothe method under 
some appropriate conditions. Bahuguna [6] has also proved similar results by using the 
Rothe method for the abstract quasi-hneax exphcit integro-differential equation 

^(t) + A{u{t))u{t) = f a{t — s)k{s,u{s))ds + f{t), 
dt Jo 

for 0 < t < T with the initial value uq in a real reflexive Banach space X, whose dual 
is uniformly convex, A{u) is a linear operator in X for each u in an open subset W of 
y where Y is also a real reflexive Banach space which is continuously and compactly 
embedded in X, a is a real-valued and / is a T-valued function defined on [0, T] and 
fc is a y- valued map defined on [0, T]xW while using the semigroup theory together 
with the contraction mapping theorem, Bahuguna [7] has established the existence of 
a unique classical solution to the quasi-hnear integro-differential equation in Banach 
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space X, 


du 

—(t) + A{t,u{t))u{t) = a{t-s)k{s,u{s))ds + f{t), 
u{0) = uo, 


where A{t, u) is a linear operator in X for each u in an open subset W of Y, maps a, f, 
k and spaces X, Y have same properties as in [6]. Also using similar techniques of paper 
[6], in [12], we have established the existence, uniqueness and continuous dependence 
of a strong solution to the quasi-hnear “implicit” integro-differential equation 

dti 

—(t) + A{u{t))u{t) = /(t, u{t), G(u){t)), 0<t<T, u(0) = Uo, 

in a reflexive Banach space whose dual is uniformly convex. 

6.2 Method of Semi-discretization in Time 

In this section, we shall prove the existence of a unique strong solution to the quasi- 
hnear integro-differential equation and its continuous dependence on the initial data. 

Let X and Y be two real reflexive Banach spaces such that the embedding Y ^ X 
is dense and continuous. Consider the following abstract quasi-linear implicit integro- 
differential equation 

^(t) + Ait, uit))uit) = fit, uit), Giu)it)), 0<t<T, u(0) = Uo, (6.2.1) 

in a Banach space X. Here 0 < T < cx), A(t,u) is a linear operator in X for each 
it, u) in [0, T] X W, where W is an open subset of H, G is a nonlinear Volterra integral 
operator defined from G([0,T];X) into G([0,r];X) and the nonlinear map / is defined 
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from [0,T] xY xY into F. By a strong solution to (6.2.1) on [O.T'], 0 < T' < T, we 
mean an absolutely continuous function u from [0,r^] into X such that u{i) G W for 
almost every (a.e. in short) t e [0,r] and satisfies (6.2.1) almost everywhere (a.e. in 
short) on [0, T']. 

In this section, we shall use the Rothe method to establish the existence and 
uniqueness results. The Rothe method, introduced by E. Rothe in 1930, is a very 
efficient tool for proving the existence and uniqueness of a solution to a linear, nonlin- 
ear parabolic or a hyperbolic problem of higher order. It has been further developed 
by Kacur [47], Necas [74], Rektorys [81], Kartsatos and Ziglar [53], Bahuguna and 
Raghavendra [6, 8, 9] and others. This method consists in replacing the time deriva- 
tives in an evolution equation by the corresponding difference quotients giving rise to 
a system of time independent operator equations. With the help of the theory of semi- 
groups and the theory of monotone operators, each of these systems is guaranteed to 
have a unique solution. An approximate solution to the evolution is defined in terms 
of the solutions of these time independent systems. After proving apriori estimates for 
the approximate solution, the convergence of the approximate solution to the unique 
solution of the evolution equation is estabhshed. In these works, either global Lips- 
chitz conditions or local Lipschitz conditions with some growth conditions on nonlinear 
forcing terms have been assumed. 

In the present work, we assume only local Lipschitz conditions on the nonlinear 
maps / and G. We first prove that the discrete points lie in a ball in X of fixed 
radius R where R is independent of the discretization parameters. Then, using the 
local Lipschitz continuity, we establish apriori estimates on the difference quotients. 
With the help of these apriori estimates, we prove the convergence of a sequence of 
approximate solutions, defined in terms of the discrete points, to a unique solution of 
the problem. 



6.2 Method of Semi-discretization in Time 


160 


Below, we state the preliminaries, the assumptions and the main result Theorem 
6.2.1. After that we state and proved the basic lemmas which are used in the proof of 
the main result. We conclude this section with a proof of Theorem 6.2.1. 

Let X and T be as earlier assumed. Let ||x||^ denote the norm of an element x 
belonging to a Banach space Z. For r > 0, let Bz{x, r) denote the open ball {z e Z : 
l|z — x||2'<7’}of the radius r and let Bz(x,r) be its closure. For an interval J of the 
real numbers, we denote by C{J\Z), WC{J\Z), Lip{J-Z) and ABS{J;Z) the spaces 
of all continuous, weakly continuous, Lipschitz continuous, and absolutely continuous 
functions from J into Z, respectively. 

For a real jS, N(Z, (3) represents the set of all densely defined linear operators L in 
Z such that if A > 0 and Aj0 < 1, then (J -f XL) is one-to-one with a bounded inverse 
defined everywhere on Z and 

||(/ + AL)-'||b(z) < (l + A/3)-‘. (6.2.2) 

where I is the identity operator on Z. The Hille-Yosida (1.3.11) theorem states that 
L G N(Z,P) if and only if -L is the infinitesimal generator of a Co-semigroup e“‘^, 
t > 0, on Z satisfying ||e'"*^||B( 2 ) ^ t > 0. A linear operator L on D{L) C Z into 
Z is said to be accretive in Z if for every u G D{L), 

< Lu,u* > > 0 for some u* G F{u), 

where < z, 2 * > is the value of z* G Z* at z G Z and F : Z -4 2^* is the duality map 
given by 

F(z) = {z* eZ* : < z,z* > = Pill = P*ll|.}- 

Here 2^* denotes the power set of Z*. If L € N{Z, 13), then the Lumer-PhiUips Theorem 
1.3.14 implies that (L + pI) is m-accretive in Z, i.e., {L + pI) is accretive and the range 
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R{L + AJ) — Z for some A > /?. If Z* is uniformly convex, then F is single valued and 
uniformly continuous on bounded subsets of Z. 

We assume, in addition, that the embedding Y ^ X is compact and the dual X* 
is uniformly convex. Further, we make the following hypotheses. 

(HI) There exist an open subset W of 7 and /? > 0 such that UoEW and 

A:[0,T] xW-^N{X,P). 

(H2) For each t,w E [O.T] x VF, we have A{t, w) G B{Y, X) for each w € W, t A{t, w) 
is the Lipschitz continuous with the Lipschitz constant La in B(Y,X)-norm.There exist 
the positive constants ij,a and 7 ^ such that for all Wi, W2 ^W, t E [0, T] and u G 7, 
we have 7 C D{A(t,w)), 

||(A(t,wi) - A(t,u;2))t;|lx < -•u;2lUI|'yl|r, 

||A(t,u;)u||x < LaMW- 


(H 3 ) There exist a linear isometric isomorphism 5:7 — > X, a map P : [ 0 , t] x W 
B{X) and the positive constants fxp and jp such that for all w, Wi, W 2 E W and 
t G [0,t], 

SA{t,w) = A{t,w)S + P{t,w)S, \\P{t,w)\\x < Ip, 

\\P{t, wi) - P{t, W2)\\x < /iplki - w^ 2 ily- 

(H 4 ) The nonlinear map G : C{[Q,T\,X) C(0,T-,X) satisfies 

(a) For all u,v E Bc{[o,T];X){Fo,r), 

||G(u) -G(u)|lcao.r];^) ^ /^G(01i'“-^llc([O.Tl;X), 
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where Hair) is a nonnegative, nondecreasing function and uo € ^([O.r];^) be defined 
by uo{t) = uq for all t G [0,T]. 


(b) For all t, s E [0,r] and u E Lip(jQ,T]]X) H 5c([o,7’];X)('ao)^')) 

|iG(u)(t) — G(u)(s)||x < 7 g( 0K ~ ^i(i + lM^/^^||L°°([o,ri;-’f))) 
where 'ycir) is a nonnegative and nondecreasing function. 

(c) Furthermore, the subspace G([0, T]; Y) of space G([0, T]; X) is an invariant subspace 
of the map G, i.e., G : G([0, T];y) — j- G([0, T"]; Y) which satisfy 

ll<^(^)WI|y < for ueSyfyo.r), 

where is a noimegative and nondecreasing function. 

In particular, we may take the operator G as a Volterra operator defined by 

G{u){t) = f a{t — s)k{s,u{s))ds 
Jo 

in which the function a is a real valued continuous function defined on [0,r] and k is 
defined on [0,T] x Y into Y and ||fc(t,iy)||y < G* for every {t,w) E [O.T] x Y, then 
map G satisfies hypotheses (c). 

(H5) The nonlinear map / : [0, T] x Y x Y Y is a bounded function: 

\\f{t,u,v)\\Y < A/(r), 

for all {t, u, v) E [0, T] x Y x Y with \\u\\y + |fyl|y < r, where A/(r) is a nonnegative 
function. Also, this map satisfies the Lipschitz condition 

\\fiil,Ui,Vi) - f{t2,U2,V2)\\x < Ai/(0ll'^(il) “ <^(^2)1 + IK - W 2 IU + Ihl -■y 2 lU], 
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for all ti, t 2 € [0, T] and all iXj, Vi € Bx{uQ,r), i = 1,2, where 0 is a real-valued 
continuous function of the bounded variation on [0,r] and /i/(r) is a nonnegative and 
nondecreasing function. 


Let R>0 he such that Wr = By (uq, R) C W. We set 


Bo 


(6.2.3) 

Ri 

= max{B,AG(3B) + lluolly}, 

(6.2.4) 

M 

= A/(B -h ||uoll -1 - Ag(3B)), 

(6.2.5) 


where 6 = j3 ||B||B(.xr) and V {4>) is the total variation of (j) on [0, T]. 

Let Zq eY and To, 0 < To < T be such that 

11‘S'wo ~ •s^olU ^ (6.2.6) 

^ofTAlkolIv + + ^] < Bq. (6.2.7) 

We note that (6.2.6) and (6.2.7) imply that 

(1 -t- e^^^)[llS'uo — aoll -I- ro{7All.Zo||r + 7plkol|x + .^}] < 2R/3. (6.2.8) 

We have the following main result for the existence, imiqueness and continuous 
dependence on the initial data of a strong solution to (6.2.1). 

Theorem 6.2.1 Let (H1)-(H5) hold. Then there exists a unique strong solution u 
to (6.2.1) on [0,To] such thatuE Lip{[0,To\;X). Furthermore, ifvo E By{uo,Ro), then 
there exists a strong solution v to (6.2.1) on [OjTo] with the initial point uq replaced by 
Vo and 

llu(t) — u(t)||x < C^||Mo“WiU) ie(0,To], 
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where C is a positive constant depending only on Tq. 


To apply the Rothe method, we proceed in the following way. Let wq = Suq. Let 
h = Tq/u for all the positive integers n > N, where N is a, positive integer such that 
e{To/N) < 1/2. For n> N, we set = uq, = uq and t] = jh for j = 1,2,..., n. 
We consider the scheme 


Su] -f A(t^,uJ_,)u] = f{tlu]_„G{u'}_,){t])), j = 1,2, 


(6.2.9) 


where, for j = 1, 2, . . . , n, n > iV, 


5u] 


um 




Uq, t = 0, 

Ui-1 + - t2.,){u2 - <_i), t € [ttiA^l 

i = l,2,...,j. 




t e [t^To]. 


( 6 . 2 . 10 ) 

( 6 . 2 . 11 ) 


The following result establishes the fact that Uj E Wr, j = 1,2,. . . ,n, n> N. 


Lemma 6.2.2 For each n> N, there exist a unique € Wr, j = 1,2, ...,n satisfy- 
ing (6.2.9). 


Proof. It follows from Lemma 1.3.28 that there exists a unique eY such that 

u'l -f hA{ti,UQ)ui =uq-\- /i/(ti,uo,G(uo)(ti)). (6.2.12) 

Applying S on both the sides in (6.2.12) using (H3) and putting w^ = Sui, we have 

(m” — ^o) + Wo)('^r ~ -^o) + hP{ti,uo){wi — zq) 

= (wq-zq) - hA{t^,uo)zQ + hP{fl,uo)zQ + hSf{t^,uo,G{uo){f^)). 
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The estimates in Lemma 1.3.28 imply that 

— 2:o||x < {l — h6) ^[11^0 — zollx + h(7A|i-2ol|y + TpII^^oIU + M)]. 

Since hO < 1/2, we have 

Ww^-ZqWx < e^^^[|ko-2;olU + M7A||^o||y + 7p||'2o||x + M)]. 

Therefore 

11^1— mollx < (1 + — -^oll + ^(7A|ko||y + 7p|kolU + Af)] < Ri 

in view of the estimate (6.2.8). Hence u” 6 Wr. Now, suppose that u" 6 Wr for 
i = 1, 2, ... ,y — 1. Again, Lemma 1.3.28 implies that for 2 < j < n, there exist a 
unique €.Y such that 

u] + hA{t], u]_,)u] = -h G{u]_,){t])). (6.2.13) 

Proceeding as before and putting w'/ — Suj, we get the estimate 

\\yjj - zoWx < e^^^[\\wj-i - zoWx + Kia\\zo\\y + IpWzqWx + M)]. 

Reiterating the above inequality, we get 

\\wj - zoWx < e^^^^[\\wo-zo\\x+jK'yA\\zQ\\Y + 'yp\\zo\\x + M)]. 

Using the fact that jh < To, we arrive at 


Ikj-f^olU < (1 -t-e^^^)[||'iao- prolix + 2o(7Apol|y+7pikolU + Af)] < R. 
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The above estimate and the equation (6.2.12) and (6.2.13) gives the required result. 
This completes the proof of the lemma. 

Lemma 6.2.3 There exists a positive constant C, independent of j, h and n such that 

< c, j = l,2,...,n; n>N. 

Proof Putting j = 1 in (6.2.9), we get 

5ui + hA(Jfl,uo){5ui) = —A{ti,uo)uQ — f(ti,UQ,G{uo){t^)). 

Using Lemma 1.3.28, we have 

IIKIU < e^^^MuoWv + M] := Co. (6.2.14) 

Prom (6.2.9) for 2 < j < n, we have 
5u^ + hA{t],u]_f){5u]) 

Using Lemma 1.3.28 again, we have 

Adding and subtracting A(t"_i,n"_ 2 )u"_i, we have 

\\Su]\\x < e2'^®[Pu^ilU + l|A(t”,t^-_i)-A(t^-i.w^i)lU!l«^^ 

+ \\{A{t]_„u]_r) - A(t^_i,u^_2))u^_ilk 
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Now, using (H2), we get 






We have 

iifajiu < t^”'[a+c,h)\\su]_,\\x+c2h 

+ ll/(i”, G(5”_.){!“) - /(t?-,. u"_,. G(5^,)(iJ..))|U] , 

(6.2.15) 


for some positive constants Ciand C2 independent oij, h and n. We note that 


l|G(a?-i)(t?) - G(fi”_,)(tj_,)k 

< f^3(3iJ)/i||5K“_ilU + 7c(3flWl ll'5<IW. (6.2.16) 


< + h\\Su].,\\x + ;*G(3H)fc||i<_,|U 

+ 7G(3fi)A(l + , , IIKIU)]. (6.2.17) 

Using (6.2.16) and (6.2.17) in (6.2.15), we obtain 


||<^w”||x < e^^'^[(l + hC'i)l|5M"_i||x + hC2 + — 0(tj_i)i 

+ /i||duJ_i|U/iG(3H)h||5u^_ilU +7G(3i?)h(l + II<^<1U)}] 


< 


^2he 


(1 + C3/1) max ||5u"|U + ^4^ + C'4|^^(t”) — ^(ty_i)l 
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Hence 


IIKk + - 4>{tU)\ + Csk , 

( 6 . 2 . 18 ) 

where (7$ is another positive constant independent of h and n. Reiterating the 
inequality (6.2.18), we get 

max 1|5<|U < + CshYlWSu'lWx + CsVicf)) + C5T0], 

where V{(p) is the total variation of (f>. Using (6.2.14), we have 

P<|U < e‘^^^+^^^^[Co + C5V{<i>) + C,To\ := C. 




This completes the proof of the Lemma. 


Now, we define the Rothe sequence of the functions {U"'} from Jq into Y by 




.,n. 


( 6 . 2 . 19 ) 


Furthermore, we define a sequence of the step functions from {—h,To] into 
Y given by 


J Uo, t E (— h, 0] 

^ Uj, t E j=l,2,...,n. 


( 6 . 2 . 20 ) 


Remark 6.2.4 We observe that X"(t) G Wr for all t G i-h,To] and n > iV. Also, 
X^{t) - U^{t) 0 in X uniformly on Jo as n 00 and {U"} are in Lip{Jo,X) with 

the uniform Lipschitz constant C. 
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Also, define 

Pit) = 

and 


j = l,2,...,n. (6.2.21) 


Lemma 6.2.5 

(a) A’‘{s,X’‘{s ~ h))X’'(s)ds = no-C/"(i)+ / /"(s)*. 

J 0 


( 6 . 2 . 22 ) 


(6.2.23) 


d 


dt ~ = /”W for iG(0,T]. (6.2.24) 

Proof, (a) Using the definitions (6.2.19), (6.2.20), (6.2.21), (6.2.22) and the equation 
(6.2.9), we have 


f A^(s,X^(s-h))X'^{s)ds 
Jo 

rU rt 

= ^”(5 - h))x^{s)ds + / A”(s, X^{s - h))X^{s)ds 

i=l •'*<-1 Jtj-l 

fU rt 

= zZ I A{ti,Ui-i)uids / A{tj,Uj-i)ujds 

Jti—i 



3-1 




i=l 


= uo-U"(t)+ f F^is)d$. 
Jo 
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(b) For i e (tj-i.t,), 

A"{t, X^t - A))X"(f) = - h))X’'(t) = 


and 


^^TTn(f\ _ ('^j %-l) 

dt ^ ^ h 

Therefore 

The proof of the lemma is complete. 

Lemma 6.2.6 There exist a subsequence {!/’”} ofllT^} and a function u in Lip{Jo, X) 
such that u in C{Jo, X) ( with the supremum norm) as m oo. 

Proof. Since, {X"(t)} is uniformly bounded in Y, the compact embedding of Y, imphes 
there exist a subsequence {X’"} of {X"} and a function u : Jq—* X such that X"^(t) — ♦ 
u{t) in X as m—^ oo. The reflexivity of Y implies that u{t) is the weak limit of X"^ in 
Y, hence u{t) lies in Y (in fact in Wr since X'”(t) is in Wr). Now, X"*(t) - U"^{t) — »■ 0 
in X, so IT^it) u{t) as m — s- oo. The uniform continuity of {[/”} on Jq implies that 
{X^} is an equicontinuous family in C{Jo,X) and the strong convergence of U'^{t) 
to u{t) in X implies that {17'"} is relatively compact in X. We use the Ascoli-Arzela 
theorem to conclude that -^u in C{Jo,X) as m oo. Since are in Lip{Jo, X) 
with a uniform Lipschitz constant, u e Lip{Jo,X). This completes the proof of the 
lemma. 
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Proof of the Theorem 6.2.1. First, we show that 

A{u{t))u{t) in X as m — > oo, where denotes the weak convergence in X. 

A^{t, X^{t - h))X^{t) - A{f u(t))u(t) 

= (X”(f, X”(t - A)) - A(t, u{t)))X'^(t) + A{t, u{t)){X’^{t) - u{t)). 

Now, for t e j= ,m; 

||(,4'”(t,X”(e - h)) - A(t,«(i)))X”'(f)|U 

= \mt]‘,uf_,)-A{t,umuT\\x 

< \mtf,uf_,)-A{tf,um’‘Tiu + \\(A(tf,u{t))-A{tMmy-Tk 

< lkl|r)l|«7_i -xWIU + iAlif -iilluriU 

= + |i«ollr)l|X”{« - ft) - «(t)|U + ij|f” - «lllX"(i)|U ^ 0, 

as m — >• oo, since X'^(t) — > u{t) in X uniformly on Jo and t. Since A{t,u{t)) € 
N{X, p), pi + A(t, u) is m-accretive in X. We use Lemma 1.3.31 and the fact that 

llA(t,«(t))(X"^(t)-u(t))|U < 27Ai?, 

to assert that A{t, u{t))X'^{t) A{t, u{t))u{t) in X and, hence, A'^it, X^{t-h))X'^{t) 
A{t,u(t))u{t) in X, as m ^ oo. To show that A{t,u{t})u{t) is weakly continuous on 
Jo, let {tfc} C Jo be a sequence such that tk fas k oo. Then, u{tk) -> u(t) in X 
ask-^oo and we follow the same arguments as above to prove the A(tk, u{tk))u{tk) 
A{t, u(t))u{t) in X as k oo. For this, consider 

\\A{tk,u{tk))u{tk) - A{t,u{'t)M'^)h ^ l|(>i(ifc,«(ifc)) -^(i,«(i)))w(ife)lk 


+ P(i,'u(t)(u(tfc)-u(t))lix 
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So 


UA{tk,u{tk)) - A{t,u{t)))u{tk)\\x < \\{A{tk,u{tk))-A{t,u{h)))u{tk)\\x 

+ ||(^(t,«(tfc)) - A{t,u{t)))u{tk)\\x 
< LA\tk - t\\\u{tk)\\x 

+ lJ'A\\u{tk) - u(t)||x(i? + ||wo|ir) 0 

as > CO since tk-^t and u{tk) — >■ u{t). We again use Lemma 1.3,31 and the fact that 

to assert that A(t,u(t))u{tk) A{t,u{t))u{tk) in X and, hence, A{tk,u{tk))u(tk) 

A(t, u{t))u{t) in W, as A: 00 . 

The Bochner integrability of A{t,u{t))u{t) can be established in the similar way 
as in lemma 1.3.36. Now from (6.2.23), for each x* G X* we have 

<U'^it),x*> = <Uo,x*>+ f <-A^{s,X'^is-h))X^{s) + r{s),x*>ds. 

Jo 

Letting m — > oo using the bounded convergence theorem, we get 

<u{t),x*> = <uo,x*>+ [ < -A{s,u{s))u{s) + f{s,u{s),G{u){s)),x* >ds. 

Jo 

The continuity of the integrand imphes that < u{t),x* > is continuously differentiable 
on Jo. The Bochner integrability of A(t, u{t))u(t) implies that the strong derivative of 
u{t) exists a.e. on Jo and 

^it)+A{t,u(t))uit) = f{t,uit),G{u)it)), a.e. onJo. 

(JjL 
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Since tt(0) = uq, w is a strong solution to (6.2.1). 

Now, we establish the uniqueness and the continuous dependence on the initial data of 
a strong solution to (6.2.1). 

Uniqueness. Let v be the another strong solution to (6.2.1) on Jq. Let U = u - v. 
Then, for a.e. t E Jq 

+ <{pI + A{t,u{t)))Uit),F{U{t))> 

= P\\Umjc+ < - A{t,vmv{t),F{U{t)) > 

+ < f{t,u{t),G{u)it)) - f{tMt),G{v){t)),F{U{t)) > . 

Using the m-accretivity of pi + A{t,w) and the assumptions on A{t,w) for (^,^y) € 
[0,T] X W and f(t,u,v), we obtain 

where Jj = [0, t] and 

Cr = P + (Xa{R + ll^olly) + + l^ciR))- 

Integrating the above inequahty on (0, t) and taking the supremum, we get 

Applying the Gronwall’s inequahty we get U = 0 on Jo . 

Continuous Dependence. Let vq € By{uq, i?o)-Then 


ll/Suo — ^o||x ^ IlSuo “ 'S'uollx + -^ollx ^ 2i?o- 
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Hence 


(1 + e )[|l5'uo 2;o|| + ro{7^||2o|ly + 7 p||zoI1x + .W”}] < 3(1 + 

= R. 

We may proceed as before to prove the existence of u" G Wr satisfying the scheme 
(6.2.1) with vij and uq replaced by u” and vq, respectively. The convergence of vj to 
v{t) can be proved in a similar manner. Let U =■ u — v. Then, following the same steps 
used to prove the uniqueness, we get for a.e. t G Jo, 

^ ^RW\\c{Jt,x)- 

Integrating the above inequality on (0,t) and taking the supremum, we get 

< \\\u{o)ex+CRj‘m\i^j,j:-,ds. 

Applying the GronwaU’s inequality, we get 

\\u\\c(j.^) < ^^(OlIU, 

where (7 is a positive constant. This proves the required result. This completes the 
proof of the theorem. 

6.3 Regular Solution 


In this section, we shall prove the existence of a unique local mild solution and a local 
classical solution to the quasi-linear integro-differential equation. 
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Let X and Y be two real Banach spaces, not necessarily reflexive as considered in 
the earlier section, such that the embedding Y X is dense and continuous. Again, 
consider the following quasi-linear integro-differential equation in X: 

d/UL 

— (t) + A(t,u(t))u(t) = /(t,u(t),G(u)(t)), 0<t<r, «(0) = uo, (6.3.1) 

where 0 < T < oo, A{t,u) is a linear operator in X, for each u in an open subset 
lY of X, G is a nonlinear Volterra operator defined from C{J,X) into C{J,X), where 
J = [0,T] and the nonlinear map / defined from J x W x W into X. We follow the 
approach of T. Kato [56, 61, 63] to establish the existence of a unique classical solution 
to (6.3.1) under the assumptions (Hl)-(H8) to be stated below. 

In this section, first we mention some notations and preliminaries. Then, we 
estabhsh the existence of a unique local mild solution using the contraction mapping 
theorem and the existence of a local classical solution to (6.3.1) under some additional 
conditions. 

Let X and Y be as mentioned before and Bz(r, 2 : 0 ) be the closer of the open ball 
zq) = {z £ Z \ llz — zo\\z < ?'} with the radius r and the center at zq in a Banach 
space Z. The space of all bounded linear operators from X to Y is denoted by B{X, Y) 
and B{X, X) is written as B{X). Let J denote the interval [0,r]. The space C"”(J, Z) 
represents the space of all m-times continuously differentiable function defined on J 
into X, m = 1, 2, ...; endowed with the supremum norm 

\\u\\cr-iJ,z)= 53 

l<i<m 

where denote the ith derivative of u with = u. Let lY be a subset of X such 
that for every {t,w) E J x W A{t, w) is a infinitesimal generators of Co-semigroups 
St,w{s), s > 0 on X. The family of operators {A(t,u;)}, (t,w) € J x W, is called stable 
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if there exists constants M > 1 and a;, known as stability constants^ such that 

p{A{t, w)) D (w, oo) for {t, w)eJxW, 


where p{A{t,w)) is the resolvent set of A{t,w) and 


j=l 


< for A> 


(jj 


B{X) 


and every finite sequence 


0 < < t2 < .... < ffc < T, Wj e W", l<j<k. 

For a linear operator 5 in X, by the part S of 5 in a subspace 7 of X, we mean a linear 
operator S with domain D{S) = {x 6 D{S) n 7 : 5x € 7} and values Sx = Sx for 
X G D{S). 

Let St,w{s), s > 0, be the Co-semigroup generated by A{t,w), {t,w) e J xW. A 
subset 7 of AT called A(t, w)-admissible if 7 is an invariant subspace of the operator 
St,w{s), s > 0, and the restriction of St,w{s) to 7 is a Co-semigroup in 7. 

For the details of the above mentioned notions, the reader may refer to the chapters 
5 and 6 in Pazy [77]. On the family of operators {A{t,'w), (t,w) G J x W}, we make 
the same assumptions {H1)-{HA) considered in §6.6.4 in Pazy [77] for a homogeneous 
quasi-linear evolution equation. For completeness, we restate those assumptions. 


(Hi) The family {A{t,w)}, {t,w) E J xW is stable. 

(H2) 7 is A{t, i(;)-admissible for {t,w) £ JxW and the family {A{t, w), {t, w) E Jx W} 
of the parts A(t, w) of A{t, w) in 7, is stable in 7 . 

(H3) For (t, w)eJxW, D{A{t, w)) D 7, A{t, w) is a bounded linear operator from 7 
to X and the map 1 1 -> A{t, w) is continuous in B{Y, X) norm [|.[| for every w eW. 
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(H4) There is a positive constant La such that 

\\-A.{t,Wi) — A{t,W2)\\Y-*X < La\\wi—W2\\., 
for every toi, W 2 ^W and 0 <t <T. 

Hue C{J,X) is W-valued and The family {A{t,w),{t,w) e J xW) oi the 
operators satisfies (H1)-(H4), then there exists a unique evolution system Uy,{t,s), 0 < 
s <t <T, in X satisfying 

(i) \\U^{t,s)\\ < (6.3.2) 

for 0 < s < t < T, where M and u are the stability constants; 
d'^ 

(ii) —Uu{t,s)w\t=.s-=A{s,u{s))w, (6.3.3) 

for w e y, 0 < s < T; 

(iii) —Uu(t,s)w = -Uu{t,s)A{s,u{s))w, (6.3.4) 

os 

for cu G y, 0 < s < T. 

Further, there exists a positive constant Cb such that for every u, u G C( J, X) with 
values in W and every y eY, we have 

\\Uu{t,s)y-U,{t,s)y\\x < CMW Mr) - v{T)\\xdT. (6.3.5) 

For details of the above mentioned results, reader may refer to the Theorem 6.4.3 and 
Lemma 6.4.4 in Pazy [77]. 


We further assume : 
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(H5) for every u e C{J, X) satisfying u{t) eWiovteJ, we have 

Uu{t,s)Y cY, t,seJ, s<t 
and Uu{t, s) is strongly continuous in Y iov s, t e J, s < t. 

(H6) The closed convex subsets of Y are also closed in X. 

(H7) The nonlinear map G : C{J,X) C{J, X) satisfy the following. 

(a) For all u, v e Bc{j,x){uo,r), 

||G(u) - G(u)|lc(j,x) < I^G{'r)\\u - u||c(j,x), 

where /laif) is a nonnegative and nondecreasing function and uq € C{J,X) is defined 
by uq = uq for all t £ J. 

(b) The subspace C{J, Y) of the space C( J, X) is an invariant subspace of the map G, 

i.e. G : C'( J, Y) C(J, Y) and satisfies 

liG(w)(t)|ly < Xair) for ueByiuo.r), 

where Xoif) is a nonnegative and nondecreasing function. 

(H8) The nonlinear map f: JxWxW-^X satisfies - 

(a) For {t,u,v) e J x {W OY) x {W OY), f{t,u,v) e Y and 

\\f{t,u,v)\\Y < A/(r), 

for all {t,u,v) £ J xW xW with ||ully + ||u||y < r, where A;(r) is a nonnegative and 
nondecreasing function. 
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(b) In both Z X, Y, the nonlinear map / satisfies the Lipschitz-like condition 

\\fik,UuVi) - f{t^,U2,V2)\\z < /i/(r)[l<^(ti) - (I>{t 2 )\ + \\U, - U2\\z + |h - V 2 \\zi 

for all ti, t 2 G [0)2^] and all Ui, Vi G By(no,r), i = 1,2, where 0 is a real-valued 
continuous function of bounded variation on [0, T] and Hf{r) is a nonnegative and 
nondecreasing function. 

By a mild solution to (6.3.1) on J = [0, T], we mean a function u G C{J,X) with 
values in W satisfying the integral equation 

u{t) = Uy,{t,0)uo + [ Uu{t,s)f(s,u{s),G{u){s))ds, t G J. (6.3.6) 

Jo 

By the Classical solution u to (6.3.1) on J = [0, T], we mean a function u G C{J^X) 
such that u{t) G F fl W for t G (0,T], u G C^((0, T],X) and satisfies (6.3.1) in X. If 
there exists a,0 < T' <T and a function u G C{J',X), J' = [0, T'], such that u is a 
mild (classical) solution to (6.3.1) on J', then u is called a local mild (classical) solution 
to (6.3.1). 

Now, we state and prove the following main result. 

Theorem 6.3.1 Suppose that uq G T and the family {A{t,w)} of linear operators for 
t e J = [0,T] and w eW = {y eY : \\y - uoHr < r}, for fixed r > 0, satisfy the 
assumptions (H1)-(H6) and A{t,w)uo gT with 

\\A{t,w)uo\\Y < C'a, (6-3.7) 

for all (t, w) £ J X W. 

Further, Suppose that nonlinear map G and f satisfy (H7) and (H8), respectively. 
Then there exists a unique local classical solution to (6.3.1). 



6.3 Regular Solution 


180 


Proof. First, we establish the existence of a unique local mild solution to (6.3.1). We 
note that from the assumption (H5), it follows that 


\\Uu{t,s)\\B(Y) < Cl, 


(6.3.8) 


ioT s <t, s,t e J and every u e C( J, X) with values in W. We choose 


To = min 


/t — —\ 

\ ’ 2CaCi' 2Ci\f{Riy 2P ! ’ 


(6.3.9) 


where 

P = Colluolly + Me“^/i/(i?i)(l + ixair)) + Co\f{Ri)T 

and 

Ri = r + lluol|y + AG(r). 

Let S be the subset of C{Jq,X) defined by 

S = {u : u€ C(Jo,X), u{0) = no, u{t) eW iovt E Jo}, 

where Jo = [0,To]. From (H6), it follows that .S' is a closed convex subset of C{Jo,X). 
Next, we define a mapping F : S S given by 

Fuit) = Uu(t,0)uo+ f Uuit,s)f{s,u{s),G{u){s))ds (6.3.10) 

Jo 

and check that F is well defined. Clearly, Fu{0) = no, Fu S C{Jo, X) and (H5) implies 
that Fu{t) G F for t e Jo and now, it remains to show that l|Fu(t) - Molly ^ 
t G Jo- Consider 


Fu{t) — Mo = Uuit, 0)mo — Mo -f 



s,u{s), G{u){s))ds. 


(6.3.11) 
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Integrating (6.3.4) in X from 0 to t, we find 


Uu{t,0)uo - uq — J Uu{t,T)A{T,u{r))uQdT. (6.3.12) 

Hence, using (6.3.7) and (6.3.8) in (6.3.12), we have 

||f/,(t,0)ao- Molly < CiCaT^ < (6.3.13) 


Also, using (6.3.8) and (H8), we have 


s)f{s, u{s),G{u){s))ds 


< CiXf{R^)To < -, 

y ^ 


(6.3.14) 


since |lM(s)||y + ||G('n)(s)|jy < Ri. Using (6.3.13) and (6.3.14) in (6.3.11), we see that 


F is well defined. Now, suppose that u,v E S, then 


Fu{t) - Fv{t) = {Uu{t, 0) - Uv{t, 0))uo 

+ / [Uu{t,s)f{s,u{s),G{u){s)) 

JQ 

- Uy{t,s)f{s,v{s),G{v){s))]ds 

= T 1 + T 2 , (6.3.15) 


where Ti and T 2 represent the first and second terms on the right hand side of (6.3.15). 
We use (6.3.5) to obtain 


l|Ti|U ^ C'olluoilyTollM - m|1c(Jo,^)- 


Now, we use (H7), (H8) and (6.3.5) and get 


2\\X 


[' Uu(t,s)[f{s,u{s),G{u){s)) - f{sMs)>G{v){s))]ds 

Jo 

+ S) - U,it, s)]fis, v{s), G{v){s))ds 

Jo X 

< [Me‘^V/(i?i)(l + fiair)) + CoXf{Ri)T]To\\u - uilccjo.^). 
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since, we have 

\\f{s,u{s),G{u){s)) - f{s,v,{s),G{v){s)\\x 

< /i/(i2i)[||n(s) - u(s)|U + liG(u)(s) - G(u)(s)lU] 
^ ^/(•^i)[||w-'y||c;(jo,x) + ItG(u) - G(u)||c(Jo.a:)] 

< l^f{Ri){l+ ^iG{r))\\u-v\\c(Jo,x)■ 

Rence, from (6.3.15), we have 


\\Fu- Fv\\c{jo,x) < PTolIu - u||c(jo,x) 

1 ,, 

Thus, F is a contraction map from S to S. Since S is closed in X, by the contraction 
mapping theorem, F has a urdque fixed point u € 5 which is the local mild solution to 
(6.3.1). 

Now, we show that u E C{Jo,Y). For s,t E Jo with s < t, we have 

u{t) - u(s) = {Uu{t, 0) - Uu{s, 0))uo 

+ f {Uu{t,v) - Uuis,r]))f{ri,u{T]),G{u){ri))dr) 

Jo 

+ Uu{t,v)f{vMv),G{u){T]))di]. 

Since Uu{t,s) is strongly continuous in F for s,t € <7, s < t so, for every e > 0 there 
exist di > 0, ^2 > 0 such that 

^ 1 ) ^2 £ <7o with |ti — ^ 2 ! ^ ^ l|Gii(t, 0) — Fu(s, O)llB(y) G 


e 

^Xf{Ri)To 


^1,^2 G Jq with |ti — t2| ^^2 ^ l!l7u(t, 7?) Uu{s, 7 ])\\b(y) G 
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Let 


Then, for s,t E Jq 


Thus, u e C{Jo,Y). 




\t~s\ < 5 =4. 


Consider the following linear evolution equation 


dv{t) 

dt 


+ B{t)v{t) = h(t), 0<t<To, 

u(0) = Uq, 


(6.3.16) 


where B{t) = A{t,u{t)) and h{t) = f{t,u{t),G{u){t)) for t E Jq and u is the unique 
fixed point of F in 5'. We note that B{t) iov t E J be the infinitesimal generator of a 
Co semigroup St(s), s>0 on X and satisfies (H1)-(H3) of Theorem 1.3.23. 

We have to prove that h E C{Jo,Y). For s, t E Jo (we assume with out loss of 
generality that s < t), we have 


||h(t)-/l(s)||y = |l/(t,u(t),G(u)(t))-/(s.u(s),C(u)(s))|ly 
+ ||GM(t)-GWW||y]. . 

As ^ is a continuous function of bounded variation on J, u € C{Jo,Y) and G(u) G 
C{Jo,Y) for u E C{Jo,Y). So, for every e > 0, there exist (5i > 0, ($2 > 0 and <53 > 0 
such that 


^ 1,^2 G Jo with |ti — t 2 l < 


\4>{t) - ((>{s)\ < 


e 
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ti,t 2 eJo with \ti-t 2\<62 =» ll^^(i)-'u(s)lly < ^ 

3/x/(i?i) 

ti,t 2 eJo with \ti-t 2 \<S 2 \\G{u){t)-G{u){s)\\Y < — r. 

3/X/(i?i) 

Let 5 = min{(5i, 52 , 53 }. Then, for s,t e Jo, we have 

\t-s\ < 5 ^ ||/i(t)-/i(s)||^ < e 

Thus, h G C{Jq,Y). Theorem 1.3.24 implies that there exists a unique function v G 
C{Jo,Y) such that v G C^(Jo/{0},X) satisfying (6.3.16) in X and v is given by 

v{t) = Uy,{t,0)uo + [ Uu{t,s)f{s,u{s),G{u){s))ds, t G Jo, 

Jo 

where Uu(t,s), 0 < s < t < % is the evolution system generated by the family 
{A(tyU{t))}, t G Jo, of linear operator s in X. The uniqueness of v imphes that v = u 
on Jo and hence, u is & unique local classical solution to (6.3.1). This completes the 
proof of the theorem. 



Chapter 7 


Scope for Further Research 


The outcome of any serious research can 
only he to make two questions grow where 
only one grew before. 


— Thorstein Veblen 

In Chapter 2 we have proved the convergence results for the approximate solution 
to a first order Sobolev type evolution equation under more general conditions than 
the ones considered by [39]. In order to show the continuity of the map 1 1 -> {Snu){t), 
we assume Holder continuity in t on A^g{t,u) only and not on the function f{t,u) 
as assumed in [39]. After proving that Sn is a strict contraction, we conclude that 
the solution Un(t) of approximate integral equation is Holder continuous which in turn 
implies the Holder continuity of the solution u{t). 

To get the approximate solution of this evolution equation we have to solve a 
system of first order ordinary differential equations satisfied by the coefficients af{t). 
Finding the solution of this system would not be so straightforward since derivatives of 
ctf are not given explicitly. It would depend on the forms of g and / . We may solve this 
system by using numerical techniques for solving ordinary differential equations. For 
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further research work we plan to consider some initial boundary value problems which 
may be reformulated as Sobolev type evolution equations considered in the present 
work. We would like to study the behavior of solution obtained numerically with the 
help of corresponding systems of first order ordinary differential equations. 


Neutral Functional Differential Equations 


Consider the boundary value problem 


d 

dt 


u{t,x) + 

rt 


+ / n(5 



6(s — t, ? 7 , x)u{s, r])dr]ds 


— t)u{s,x)ds + ai{t,x), t>0, 


J — OO 

u{t, 0) = u{t, tt) = 0, t>0, 

u{9,x) = (t>{d,x), 0 < 0, 0 < X < TT, 


dx'^ 


u{t,x) + ao{x)u{t,x), 


0 < X < TT, 


where the functions Uo, a, Ci, bi and <p satisfy appropriate conditions. We can refor- 
mulate this problem in the form of an abstract neutral functional differential equation 
with unbounded delay of the form 

^(x(t) -f G{t, xt)) = Ax{t) -I- F{t, xt), t>a, 

Xa = (f £0., ^ (7.0.1) 

in a Banach space X where is an open subset of some abstract phase space B. The 
history function Xt : (oo,0] X given by Xt{9) := x{t + 9) belongs to B, and the 
nonlinear functions i^, G : [0, a] x — »■ X are continuous and A generates an analytic 
semigroup T{t) of bounded linear operators on X and 0 < cr < a. 

A mild solution for (7.0.1) is a function x : (-oo,cr -h 6) X, 6 > 0, such that 
Xa = (p, the restriction of x(.) on the interval [cr, cr + fe) is continuous and for each 
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a <t<a + bthe function AT{t - s)F(5,a:,), s e [a,t), is integrate and satisfies the 
integral equation 


x{t) = T{t-a)[cp{0) + F{a,ip)]-F{t,Xt) 

- ^T{t - s)F{s, x,)ds + T{t - s)G{s, Xt)ds, a <t. 

We plan to establish similar approximation results for the equation (7.0.1) as we have 
done for the Sobolev type case. Further, we want to obtain error estimates in the 
Faedo-Galerkin procedure for the evolution equations considered in the present study. 

In Chapter 3 and Chapter 5, we have proved the convergence of the approximate 
solution to the second order semi-linear evolution and integro-differential equations. In 
these cases, we get a pair of systems of first order ordinary integro-differential equations 
which can not be easily solved. How to solve these systems will depend on the form of 
nonlinear functions / and 5 . If / and g depend on unknown function only and not on 
the time derivative of the unknown function, then we do not require to solve the second 
system. We may get the solution by solving the first system only which involves af . 
We have planed to develop a method to solve these type of systems. 


In Chapter 6, we have proved the existence of a unique strong solution which con- 
tinuously depends on the initial data with the help of the method of semi-discretization 
in time. Using this method we would Uke to prove the existence result for the sec- 
ond order semi-linear integro-differential equation considered in Chapter 5. We plan to 
consider the following discretization scheme. 

j-i 

+ ASu] + Bu]_^ = /; + h 

i=0 


where 




h? 


j = l,2---n; 
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= Kt" - <?). 1 < ! < J < ni 

We would like to obtain the error estimates also, i.e., the problem of estimating the 
difference between the actual value of the solution u and its approximation Un obtained 
at the point of discretization by Rothe’s method which is important from the numerical 
point of view. 


Second Order Nonlinear Differential Equations 

Let X and Y be two real Hilbert spaces such that Y C X and let the inclusion 
mapping of Y into X be continuous and densely defined. Let Y' denote the dual of Y. 
We identify X with its own dual X'. Consider the problem 

^(t) + A(^)(t) + B«W3/(t,«{t),^(t)). a.e. on (0,r), 

«(0) = HO, ^(0) = Ho, (7.0.2) 

where A is a maximal monotone set in y xY' and A E L{Y, Y') is assumed to satisfy 
the following coercivity condition 

(Au,u) + Q;|up > w||uf, for every Y, 

where a E B} and u; > 0. The nonlinear map / is defined from (0, T) x T xY to H. 

In chapter 3 we have established the convergence results for the approximate so- 
lutions of a second order semi-linear evolution equation. Our aim is to extend these 
results for the second order nonlinear differential inclusion (7.0.2). 
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